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ABSTRACT 


The  Non-Intrusive  Load  Monitor  (NILM)  provides  a  method  of  measuring  component  perfor¬ 
mance  and  source  power  quality  through  a  single  point  of  entry  in  the  power  distribution  system. 
A  study  was  performed  utilizing  the  NILM  onboard  three  different  ships  (along  with  experimen¬ 
tation  in  the  laboratory)  to  determine  its  effectiveness  in  determining  the  state  of  mechanical 
systems  through  analysis  of  electrical  power  data. 

Data  collected  from  the  Auxiliary  Seawater  System  onboard  the  USCGC  SENECA  indicate 
that  the  NILM  is  able  to  predict  several  faults  (clogged  pump  inlet  strainers,  faulty  motor/pump 
coupling,  fouled  heat  exchangers)  as  well  as  provide  a  backup  indication  of  flow  levels  to  heat 
loads.  Data  collected  from  the  Sewage  System  of  the  SENECA  indicate  several  metrics  which  can 
be  applied  to  cycling  systems  in  general  to  differentiate  between  periods  of  heavy  usage  and  fault 
conditions.  Finally,  data  collected  from  the  Steering  System  of  the  SENECA  as  well  as  a  yard 
patrol  boat  operated  by  the  Office  of  Naval  Research  shows  the  potential  of  the  NILM  to  be  used 
as  a  control  system  without  the  need  for  separate  mechanical  transducers. 

Thesis  Supervisor:  Steven  B.  Leeb 

Title:  Associate  Professor  of  Electrical  Engineering  and  Computer  Science 
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Chapter  1 
Introduction 


1.1  Motivation  for  Research 

The  U.S.  Navy  has  a  growing  dependence  on  technology  to  perform  its  mission,  making  ships 
ever  more  complex.  As  a  result  of  the  number  and  complexity  of  shipboard  systems,  it  becomes 
more  important  to  monitor  them  constantly.  The  research  presented  in  this  thesis  is  the  first  step 
towards  developing  an  automated  electrical  load  monitoring  system  for  shipboard  use. 

1.1.1  The  Electric  Ship 

There  is  a  growing  trend  of  “electrification”  in  U.S.  Navy  ships,  most  notably  in  the  development 
of  the  integrated  power  system.  In  the  integrated  power  system,  ships  no  longer  have  separate 
main  engines  to  provide  power  for  propulsion  and  generators  to  provide  electricity.  The  ship  that 
uses  an  integrated  power  system  would  have  only  generators  which  supply  electrical  power  for 
ship  service  power,  combat  system  power,  and  propulsion  power  (via  an  electric  drive  system). 
This  setup  allows  the  generators  to  be  loaded  closer  to  their  rated  power  output  in  order  to 
operate  the  prime  mover  for  these  generators  at  higher  efficiency  (and  therefore  lower  specific  fuel 
consumption). 

Another  way  that  the  electrical  loading  on  naval  ships  is  increasing  is  through  transition 
from  mechanical  systems  to  electrical  systems.  Advances  in  computing  and  power  electronics 
have  made  possible  electrical  equipment  applications  which  have  traditionally  been  mechanical, 
hydraulic,  or  pneumatic.  For  example,  steam  catapults  have  traditionally  been  used  to  launch 
aircraft  from  aircraft  carriers.  The  electromagnetic  aircraft  launch  system  (EMALS),  which  uses 
a  linear  induction  motor,  will  replace  the  traditional  steam  catapults  on  the  next  generation  of 
aircraft  carrier. 

The  final  way  that  the  electrical  loading  on  naval  ships  is  increasing  is  through  the  increased 
computing  power  required  on  ships.  This  is  due  not  only  to  the  increasing  abilities  of  the  combat 
system  suite,  but  also  because  of  the  more  extensive  use  of  information  technology. 

The  importance  of  electrical  monitoring  grows  as  the  number  of  electrical  loads  onboard 
warships  increases.  As  the  shipboard  distribution  systems  get  more  complicated,  it  will  be  more 
difficult  to  keep  track  of  the  loading  and  usage  of  all  electrical  loads;  hence  the  need  for  an 
electrical  load  monitoring  system. 
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1.1.2  A  Growing  Need  for  Information 

Along  with  the  increased  use  of  electrical  equipment  on  ships,  there  is  a  trend  to  make  ships  more 
automated.  This  would  not  only  allow  for  reduced  manning  (and  therefore  reduce  operating  cost) 
but  would  also  make  ship  systems  more  reliable.  Automation  of  equipment  has  two  fundamental 
aspects,  control  and  monitoring.  Since  it  is  impossible  for  an  automated  system  to  know  what 
control  signals  to  send  unless  it  knows  the  current  state  of  the  system  being  controlled,  monitoring 
the  state  of  the  system  becomes  a  crucial  part  of  automation. 

In  addition  to  monitoring  systems  as  part  of  automation,  there  is  also  an  increasing  need  to 
collect  information  for  the  purpose  of  condition-based  maintenance.  The  preventive  maintenance 
system  traditionally  used  by  the  U.S.  Navy  follows  the  idea  that  equipment  should  be  serviced  at 
a  regular  (and  conservative)  interval  in  order  to  prevent  failures  from  occurring.  The  drawback 
of  this  philosophy  is  that  equipment  which  has  a  low  level  of  usage  can  be  serviced  more  often 
than  is  necessary  while  equipment  which  is  heavily  used  may  not  be  serviced  before  failure.  As 
a  cost-saving  measure,  there  is  a  growing  emphasis  on  condition-based  maintenance  in  which 
equipment  is  serviced  based  on  its  actual  condition.  Of  course,  the  condition  of  equipment  can 
only  be  known  if  that  equipment  is  continuously  monitored. 

1.2  Previous  Shipboard  Monitoring 

It  would  be  useful  at  this  point  to  briefly  review  some  of  the  past  efforts  in  monitoring  ship 
systems. 

1.2.1  Mechanical  Monitoring  of  Mechanical  Equipment 

This  kind  of  monitoring  is  the  simplest  and  has  been  going  on  for  the  longest  period  of  time.  With 
the  transition  from  sail  to  steam,  it  became  necessary  to  monitor  the  parameters  of  the  steam 
plant.  This  monitoring  continues  on  mechanical  systems  today,  with  a  watchstander  walking 
around  looking  at  pressure  and  temperature  indications  and  writing  down  the  readings  into  a  log 
once  per  hour.  In  some  cases  these  indications  are  purely  mechanical,  such  as  a  Pitot  tube  (as 
described  in  [1])  used  to  measure  total  pressure  in  a  fluid  system.  In  some  cases  the  indication 
involves  a  transducer  to  convert  from  one  form  of  energy  to  another,  as  in  a  thermocouple  used  to 
measure  temperature.  Whether  a  direct  mechanical  indication  or  a  transducer,  the  information 
recorded  is  always  mechanical  in  nature. 

There  has  been  a  recent  increase  in  automated  log  taking  for  mechanical  systems  in  which  a 
central  processor  collects  information  from  transducers  and  automatically  records  it.  An  example 
of  this  is  a  system  onboard  the  U3CGC  SENECA  that  records  any  changes  in  main  engine 
rpm.  Although  this  alleviates  some  of  the  burden  on  the  watchstanders,  it  does  not  provide  any 
analysis  or  control  functions.  Watchstanders  are  still  required  to  analyze  trends  in  data  and  take 
the  appropriate  actions. 

1.2.2  Electrical  Monitoring  of  Electrical  Equipment 

In  addition  to  monitoring  the  electrical  parameters  of  the  generating  equipment,  many  other 
parameters  are  now  being  monitored  onboard  ships  as  described  below. 

With  the  push  towards  electric  drive  systems,  it  has  become  necessary  to  develop  control  and 
monitoring  systems  for  extremely  large  electric  motors.  The  system  presented  in  [2]  uses  “search 
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coils”  to  monitor  the  amount  of  flux  in  a  large  electric  motor  as  an  indicator  of  motor  condition. 
Flux  information  is  collected  by  the  system  and  compared  to  baseline  flux  levels  in  an  effort  to 
detect  faults  in  the  motor’s  rotor.  For  example,  the  partial  demagnetisation  of  a  single  magnet 
(from  a  total  of  38  pole-piece  pairs  mounted  on  the  rotor)  can  be  detected  as  a  change  in  the  flux 
amplitude  over  a  single  rotation  of  the  rotor. 

■'  The  arc-fault  detection  system  described  in  [3]  is  a  way  of  detecting  the  existence  of  arc  faults 
hi  electrical  switchboards.  As  described  by  the  paper,  arc  faults  were  becoming  an  increasing 
problem  as  the  source  for  fires  in  submarine  switchboards  in  the  1970’s.  The  faults  occur  because 
of  faulty  connections  due  to  corrosion,  faulty  fastenings,  and  vibration  as  well  as  contamination 
and  foreign  objects.  The  circuit  breakers  installed  in  switchboards  do  not  open  for  an  arc  because 
the  arcs  often  draw  much  less  current  than  the  currents  already  present  in  the  switchboard.  As 
a  result,  the  arc  fault  detection  system  that  was  developed  uses  photo  and  pressure  detectors  to 
sense  the  presence  of  an  arc.  At  the  time  that  [3]  was  written,  the  system  has  accumulated  500 
ship-years  of  operation  without  a  single  erroneous  operation. 

1.2.3  The  ICAS  system 

Most  of  the  information  in  this  section  is  taken  directly  from  [4] 

The  Integrated  Condition  Assessment  System  (ICAS)  is  a  Commercial  Off  the  Shelf  (COTS) 
software  product  designed  to  provide  a  computerized  engineering  tool  to  eliminate  paper  logs, 
implement  Condition  Based  Maintenance  (CBM)  via  expert  analysis  and  improve  availability 
and  reliability.  The  ICAS  is  a  Microsoft  Windows  NT  based  system  running  on  state-of-the-art 
ruggedized  commercial  workstations.  It  continuously  assesses  machinery  performance  conditions 
and  transmits  the  information  over  a  fiber  optic  network  for  shipboard  or  shore  side  expert  review 
action  and  planning.  Thus  ICAS  provides  maintenance  personnel  with  real-time  data  display 
at  these  consoles,  continuous  performance  assessment  and  trends,  alarms,  equipment  vibration 
information,  continuous  expert  analysis  and  diagnostics.  The  ICAS  shell  is  loaded  with  a  hull 
specific  Configuration  Data  Set  (CDS)  outlining  the  failure  and  maintenance  modes  of  the  main 
propulsion,  electrical,  auxiliary  and  combat  systems  equipment.  These  failure  and  maintenance 
modes  are  defined  by  Navy  subject  matter  experts  including  ship’s  force,  in-service  engineering 
agents  and  fleet  technical  support  personnel.  Inputs  to  ICAS  consist  of  both  manually  collected 
data  and  on-line  sensors  (including  those  parameters  collected  by  MCS  and  shared  over  the 
common  LAN). 

Although  ICAS  seems  very  capable,  it  actually  has  several  severe  limitations.  The  biggest 
drawback  is  that  ICAS  is  software  for  analyzing  and  trending  data,  not  a  data  collection  system. 
The  Portable  Data  Terminal  (PDT)  and  Portable  Diagnostic  Aid  (PDA)  are  hand  held  units 
utilized  to  collect  manual  data  and  upload  to  the  ICAS  workstation  where  it  can  be  analyzed, 
trended,  and  logged.  The  PDT  is  used  to  enter  log  sheet  data  into  the  ICAS  system.  Sensors 
are  typically  grouped  to  represent  a  paper  log  sheet.  The  operator  punches  data  in  on  the  PDT’s 
alphanumeric  keypad  and  then  downloads  the  data  into  the  ICAS  system  via  the  Data  Transfer 
Cable.  The  PDA  is  used  for  acquiring  and  analyzing  vibration  data  on  rotating  machinery.  It  is 
manufactured  by  Dynamic  Instruments  and  has  an  OEM  model  #  of  DI-306.  The  unit  is  capable 
of  acquiring  both  broadband  and  narrowband  vibration  signatures.  The  data  is  acquired  via  the 
PDA  and  downloaded  into  the  ICAS  system  for  analysis  and  trending.  The  only  data  which  is 
automatically  collected  by  ICAS  is  from  the  Machinery  Control  System  (MCS).  The  data  collected 
by  ICAS  from  MCS  is  not  related  to  the  status  of  the  equipment,  it  is  only  a  record  of  control 
signals  sent  to  the  equipment. 
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1.3  Non-Intrusive  Load  Monitoring 

Non-Intrusive  Load  Monitoring  research  has  been  in  progress  at  the  Massachusetts  Institute  of 
Technology’s  Laboratory  for  Electromagnetic  and  Electronic  Systems  (LEES)  for  the  last  twenty 
years.  The  Non-Intrusive  Load  Monitor  (NILM)  provides  a  method  of  measuring  component 
performance  and  source  power  quality  through  a  single  point  of  entry  in  the  power  distribution 
system  using  only  one  voltage  and  one  current  transducer.  The  NILM  can  be  used  in  single  or 
three-phase  systems  over  a  wide  range  of  both  voltage  and  current.  The  NILM  has  the  ability  to 
measure  multiple  loads  simultaneously  by  employing  a  Pentium  based  personal  computer  running 
a  modified  version  of  the  Linux  operating  system.  The  NILM  is  relatively  easy  to  build,  install, 
and  run  (see  Chapters  2  and  3  for  more  information).  i  ,  ]  .  . 

It  has  been  determined  that  since  the  NILM  has  a  high  data  collection  rate  (8  kHz  for  voltage 
and  current,  120  Hz  for  power  envelopes),  it  can  detect  very  short  duration  transients.  It  has 
also  been  determined  through  experimentation  that  if  the  electrical  load  in  question  is  tied  to  a 
mechanical  system  (e.g.  a  pump,  fan,  etc.),  then  changes  in  the  mechanical  system  can  have  quan¬ 
tifiable  impacts  on  the  electrical  signature  of  that  load.  Best  summarized  in  [5],  the  application 
of  the  NILM  to  date  has  been  primarilay  in  HVAC  systems. 

1.4  Putting  it  all  Together:  The  NILM  on  Ships 

The  following  observations  are  made  from  the  preceding  sections: 

•  Monitoring  of  mechanical  parameters  onboard  ship  requires  crewmembers  to  read  unsophis¬ 
ticated  instruments  which  are  not  accurate  or  prone  to  failure.  In  the  case  where  mechanical 
parameters  are  monitored  electronically,  a  dedicated  transducer  is  required  for  each  param¬ 
eter.  This  transducer  is  subjected  to  the  actual  environment  that  it  is  monitoring. 

'  •  Monitoring  of  electrical  parameters  onboard  ship  has  to  date  been  conducted  at  the  genera¬ 
tion  or  distribution  level  with  the  purpose  of  protecting  the  electrical  system  from  electrical 
faults.  ; 

•  The  ICAS  system  has  the  ability  to  analyze  and  trend  equipment  loading,  but  does  not 
automatically  collect  data. 

Clearly  a  better  system  is  needed  onboard  U.S.  Navy  warships  to  monitor  the  electrical  loads 
as  called  for  in  [6].  Such  a  system  should  have  the  following  attributes: 

•  Fault  Detection.  Data  collected  would  continuously  be  compared  to  a  library  of  normal 
and  abnormal  conditions  to  report  the  existence  of  a  fault  condition  to  an  operator.  The 
operator  “...must  be  presented  with  information  in  a  manner  which  assists  quick  decision 
making  and  error-free  command  inputs.”  [7] 

•  Automatic  Corrective  Action.  Should  a  fault  condition  be  severe,  the  system  would  be  able 
to  shift  over  to  a  redundant  piece  of  equipment  before  catastrophic  failure  of  the  faulted 
equipment  occurred. 

•  Fault  Correction.  Should  a  catastrophic  failure  occur,  data  collected  by  the  system  would 
be  stored  in  a  robust  format  to  be  used  later  to  determine  the  cause  of  the  fault. 


18 


•  Usage  Trending.  Tracking  of  equipment  usage  would  be  used  to  ensure  equal  run  time  for 
systems  with  redundant  electrical  loads. 

The  NILM  would  be  a  perfect  starting  point  for  such  a  system  for  many  reasons: 

:  •  It  can  gather  mechanical  information  about  a  system  from  electrical  information  without 

the  need  for  dedicated  mechanical  transducers.  The  electrical  sensors  used  are  solid-state 
••  and  are  therefore  very  reliable. 

•  Since  the  sensors  used  are  located  inside  an  electrical  panel  or  controller,  they  are  not 
subjected  to  the  harsh  conditions  present  in  most  mechanical  systems. 

•  Electrical  monitoring  is  accomplished  using  existing  power  cabling  so  backfitting  the  system 
on  older  ships  would  not  require  running  an  excessive  amount  of  cabling.  It  would  also  not 
require  creating  penetrations  into  mechanical  systems. 

•  All  NILM  installations  used  to  date  are  COTS  and  so  are  very  cost  effective. 

1.5  Objectives  and  Outline  of  Thesis 

The  research  in  this  thesis  represents  the  first  step  in  developing  the  monitoring  system  mentioned 
above.  Since  electrical  load  monitoring  of  this  kind  has  never  been  attempted  on  a  ship  (to  the 
author’s  knowledge),  the  scope  was  very  large.  The  goal  was  to  find  as  many  ships  as  possible 
with  complex  electromechanical  systems  and  install  NILMs  on  those  systems.  The  data  would 
be  analyzed  to  determine  which  systems  the  NILM  was  best  suited  to  monitor.  A  total  of  12 
systems  were  monitored  on  3  different  ships  with  almost  200  gigabytes  of  data  collected.  As  to  be 
expected  with  research  of  this  kind,  some  systems  are  more  amenable  to  non-intrusive  monitoring 
than  others.  Analysis  of  some  data  has  been  deferred  to  the  future. 

Chapters  2  and  3  discuss  brief  theory  and  construction  of  the  NILM.  Chapters  4,  5,  and 
6  discuss  results  for  data  taken  onboard  the  USCGC  SENECA  for  the  sewage  system,  auxiliary 
seawater  system,  and  steering  system,  respectively.  Chapter  6  also  discusses  results  for  the  steering 
system  onboard  a  yard  patrol  craft  operated  by  the  Office  of  Naval  Research.  Chapter  7  presents 
an  outline  of  future  work  that  will  be  conducted  in  upcoming  years. 
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Chapter  2 
Background 


2.1  Measuring  Real  and  Reactive  Power 

The  NILM  uses  measurements  of  voltage  and  current  to  produce  estimates  of  real  and  reactive 
power  consumption  (henceforth  referred  to  as  P  and  Q,  respectively).  For  a  single  phase  system, 
the  NILM  is  supplied  with  the  voltage  from  line  to  neutral  and  the  current  flowing  to  any  loads 
downstream  of  the  monitoring  point.  Real  power  (P)  is  produced  by  current  which  is  in-phase 
with  voltage  and  reactive  power  (Q)  is  produced  by  current  which  is  90  degrees  out  of  phase  with 
voltage.  It  also  estimates  the  higher  harmonic  content  of  current  to  produce  the  first  three  odd 
harmonics  of  real  and  reactive  power  (referred  to  as  Pk  and  Qk  for  the  kth  harmonic  of  P  and  Q, 
respectively) . 

If  the  same  sinusoidal  voltage  is  applied  to  a  resistor,  capacitor,  and  inductor,  the  NILM 
would  produce  positive  P  with  zero  Q  for  the  resistor,  Positive  Q  with  zero  P  for  the  inductor, 
and  negative  Q  with  zero  P  for  the  capacitor  in  accordance  with  the  convention  proposed  in  [8]. 
The  convention  of  using  a  power  angle  measured  in  the  clockwise  direction  from  the  real  axis  is 
employed,  so  that  the  power  angle  for  an  inductor  is  +90  degrees  while  the  power  angle  for  a 
capacitor  is  -90  degrees.  The  summary  of  phase  relationships  between  voltage  and  current  for  the 
three  components  is  shown  in  Figure  2.1. 

For  a  more  in-depth  description  of  the  NILM  calculations  used  to  determine  P  and  Q,  see  [9]. 
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Figure  2.1:  Voltage  and  Current  Phase  Relationships  for  3  Components 
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2.2  Adaptation  to  a  3-Phase  Ungrounded  Power  Sys¬ 
tem 

For  a  three-phase  ungrounded  power  system,  as  one  would  find  on  a  ship,  things  become  more 
complex.  Since  there  are  three  different  AC  voltages  (each  out  of  phase  by  ±  120  degrees  with 
respect  to  the  other  two)  and  no  neutral  wire,  the  voltage  supplied  to  the  NILM  must  be  line- 
to-line  (between  phases).  Figure  2.2  represents  a  3-phase  generation  scheme  showing  phase  and 
line- to-line  voltages  [10]. 

Im  ‘ 


Re 


Figure  2.2:  3-Phase  Voltage  Generation  Scheme 

For  a  purely  resistive  load  in  a  3-phase  power  system,  Figures  2.3  and  2.4  apply.  The  3-phase 
voltage  source  is  shown  in  a  “wye,J  configuration,  while  the  resistive  load  is  shown  In  a  “delta55 
configuration  as  is  typical  onboard  a  ship1  [11].  In  Figure  2.4,  line-to-line  voltages  are  shown 
(although  differ  from  the  previous  figure  in  that  they  have  been  moved  to  the  origin)  along  with 
line  currents  that  would  be  flowing  from  the  source  [10].  Note  that  line  currents  are  in  phase  with 
phase  voltages  (scaled  by  Ohm5s  Law). 


Figure  2.3:  Wye-Delta  Configuration  of  Resistive  Load 

It  is  noted  from  Figure  2.4,  there  are  now  three  voltage-current  pairs  that  show  a  90  degree 
phase  shift,  (VbcVa)?  (FcaVb)?  and  (Vks,Ic)*  From  the  previous  discussion  of  a  single  phase 
system,  if  any  one  of  these  voltage-current  pairs  was  provided  as  input,  the  NILM  would  output 

1The  difference  between  wye  and  delta  configurations  is  only  in  their  voltage  and  current  amplitudes, 
not  their  phase  relationships  for  an  ungrounded  system. 
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Figure  2.4:  Wye-Delta  Resistive  Load  Phasor  Diagram 


negative  Q  and  zero  P.  Since  the  load  in  question  is  purely  resistive,  power  must  be  positive  P  and 
zero  Q.  Fortunately,  since  the  NILM  outputs  power  data  in  table  form,  the  only  adjustment  that 
needs  to  be  made  is  one  of  interpretation.  The  column  that  represented  P  for  the  single-phase 
system  now  represents  Q  in  the  three-phase  system  and  the  column  that  represented  Q  in  the 
single-phase  system  represents  -P  in  the  three-phase  system.  The  relationship  between  line-to-line 
voltages  and  line  currents  is  exploited  by  other  commercially  available  products  (such  as  the  Fluke 
39  powermeter  which  is  used  to  calibrate  the  NILM  power  measurements  [12]).  The  connection 
of  the  Fluke  powermeter  is  shown  in  Figure  2.5. 


Figure  2.5:  Fluke  39  Powermeter  3-Phase  Power  Measurement 


-  There  is  a  drawback  to  this  method  of  three-phase  data  collection  in  that  it  applies  only  to  a 
balanced  system.  Stating  that  power  measurements  taken  from  these  orthogonal  voltage-current 
pairs  is  proportional  to  the  total  power  consumed  by  a  load  assumes  that  each  phase  is  drawing 
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the  same  amount  of  power.  In  practice  this  is  generally  and  approximately  true,  but  if  a  load  is 
unbalanced,  important  dynamic  information  could  be  lost  due  to  a  poor  choice  of  the  monitored 
voltage-current  pair.  The  current  hardware  configuration  of  the  NILM  does  not  allow  for  summing 
of  individual  phase  loading.  An  example  of  currents  in  a  three-phase  load  onboard  the  research 
vessel  Oceanus  is  shown  in  Figure  2.6.  Note  that  the  three  phases  of  the  motor  are  slightly  out 
of  balance  as  indicated  by  the  oscillation  of  the  current  sum  around  zero. 


Figure  2,6:  Example  of  3-Phase  Induction  Motor  Currents 


Chapter  3 

Hardware  and  Software  for  a  Typical 
Installation 


3.1  Introduction 

A  complete  NILM  setup  is  comprised  of  commercially  available  hardware  that  is  relatively  easy  to 
acquire,  assemble,  and  install.  Commercially  available  hardware  was  chosen  because  of  its  wide 
availability  and  low  cost.  A  setup  consists  of  a  NEMA-type  enclosure  to  house  the  measuring 
transducers  and  their  associated  power  supply,  a  personal  computer  (PC)  with  keyboard  and 
monitor,  and  an  uninterruptible  power  supply.  A  block  diagram  for  a  NILM  setup  is  shown  in 
Figure  3.1  while  a  photograph  of  a  typical  installation  is  shown  in  Figure  3.2 
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From  Ship  Service  Pouer 


Figure  3.1:  Block  Diagram  of  Complete  NILM  Setup 

The  personal  computer  used  in  NILM  setups  runs  the  Linux  operating  system.  Since  the 
source  code  for  Linux  is  freely  available,  custom  software  to  be  used  in  conjunction  with  the 
NILM  can  be  easily  developed  at  low  cost.  Linux  is  also  a  very  reliable  system,  which  is  especially 
important  if  long  periods  of  automatic  data  collection  are  anticipated. 

Installation  of  a  NILM  setup  can  be  accomplished  using  readily  available  hardware  and  tools. 
For  the  shipboard  installations  discussed  in  this  thesis,  the  hardware  was  attached  to  panels, 
bulkheads,  or  mounting  plates  with  nuts  and  bolts.  Cables  were  strung  in  exising  cableways. 
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Figure  3.2:  Typical  NILM  Setup 


3.2  Measuring  Hardware 

Voltage  and  current  inputs  to  the  PC  are  generated  by  two  transducers.  Each  of  the  transducers 
employs  the  Hall-effect  to  produce  output  current  proportional  to  input  voltage  or  current.  In 
each  case,  the  transducers  are  powered  across  two  of  their  terminals  by  ±  1 2 V  in  parallel  with 
capacitors1.  The  specific  transducers  used  for  this  research  were  the  LEM  LV  25-P  to  measure 
voltage  and  either  the  LEM  LA  55-P,  LA  205-S,  or  LA  305-S  to  measure  current.  Specifications 
for  these  transducers  are  listed  in  Appendix  A. 

In  the  case  of  the  voltage  transducer,  pins  are  also  provided  for  input  voltage  (line  and  neutral) 
and  output  current.  The  transducers  are  mounted  on  a  board  in  series  with  resistors  that  are 
selected  based  on  input  voltage  to  prevent  exceeding  the  transducer’s  maximum  current  rating 
of  10  mA.  A  schematic  and  photograph  of  the  voltage  board  is  shown  in  Figures  3.3  and  3.4. 
The  voltage  board  was  originally  designed  to  measure  three  phases  independently,  but  only  one 
line-to-line  voltage  measurement  was  used  for  all  shipboard  installations. 

The  current  transducer  is  placed  directly  in  the  panel  of  the  load(s)  being  measured.  An 
insulated  three-wire  cable  is  routed  to  the  transducer  to  provide  DC  voltage  and  carry  the  output 
signal  back  to  the  NILM.  A  current  transducer  is  shown  in  Figure  3.5  inside  of  a  pump  controller 
(voltage  leads  are  also  shown  connected  to  the  other  two  phases). 

*To  account  for  any  “noise”  in  the  power  supplies  , 
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Figure  3.3:  Voltage  Sensing  Schematic 


Figure  3.4:  Assembled  Voltage  Sensing  Board 


Figure  3.5:  Installed  Current  Transducer 


3.3  Setup/Installation 

The  setup  of  the  NILM  involves  using  the  voltage  board  and  current  transducer  together  with  a 
power  supply  and  wiring  terminal  board  inside  a  NEMA  enclosure.  The  line  and  neutral  portions 
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of  a  grounded  120V  AC  power  cord  are  attached  to  the  power  supply  while  the  ground  is  attached 
to  the  voltage  board.  The  DC  voltages  output  from  the  power  supply  are  connected  to  the 
voltage  board  and  a  terminal  strip  to  supply  the  current  transducer.  The  power  supply  ground  is 
connected  to  the  wiring  terminal  board. 

Outputs  from  the  voltage  board  and  current  transducer  are  connected  to  the  wiring  terminal 
board.  Voltage  drops  across  the  terminal  board  are  produced  by  placing  resistors  across  the 
terminals  for  the  transducer  outputs  and  the  power  supply  ground.  The  resistors  must  be  chosen 
carefully  based  on  the  output  of  the  transducers  to  produce  a  voltage  drop  within  the  range 
required  by  the  data  collection  card  inside  the  computer.  Reference  resistors  must  also  be  placed 
between  the  power  supply  ground  and  pin  27  of  the  wiring  terminal  board.  The  PCI  1710  analog- 
to-digital  data  card  used  converts  an  analog  signal  in  the  range  of  -5  V  to  +5V  to  digital  “counts” 
in  the  range  of  0  to  4096  counts.  Table  3.1  specifies  the  electrical  connections  required  for  a  NILM 
setup. 


Table  3.1:  Electrical  Connection  Scheme  Required  for  NILM 
Setup 


Component 

Terminal 

Power  Supply  +12V 

Voltage  board  +12V  and  junction  strip  terminal  1 

Power  Supply  -12  V 

Voltage  board  -12V  and  junction  strip  terminal  3 

Power  Supply  Ground 

SCSI-II  terminals  33,  34 

Current  Transducer  +12V 

Junction  strip  terminal  1 

Current  Transducer  -12V 

Junction  strip  terminal  3 

Current  Output 

SCSI-II  terminal  67 

Voltage  Output 

SCSI-II  terminal  68 

Measurement 

SCSI-II  terminals  68  and  34 

Resistors 

SCSI-II  terminals  67  and  33 

Reference 

SCSI-II  terminals  34  and  27 

Resistors 

SCSI-II  terminals  33  and  27 

3.3.1  Resistor  Selection 

The  following  example  illustrates  proper  selection  of  resistors  and  a  current  transducer  for  use 
with  a  load  that  draws  6.5  A  in  steady  state  with  450  V  (rms  line-to-line)  applied: 

1.  The  voltage  board  should  contain  three  resistors  in  parallel  at  the  input  of  the  voltage 
transducer.  A  good  rule  of  thumb  for  sizing  each  of  the  resistors  in  the  voltage  board  is 
to  use  resistances  approximately  equal  (in  Ml)  to  the  RMS  voltage  being  measured.  For 
example,  to  measure  450  V  rms  line-to-line  voltages  on  ships,  three  470  kCl  resistors  were 
soldered  into  the  voltage  boards  used.  If  the  three  parallel  resistors  used  are  each  470  fcO 
(equivalent  resistance  =  156.67 kQ),  then  the  rms  current  through  the  transducer  would  be 
2.87  mA  (450V/156.67M2). 

2.  The  rms  output  current  of  7.18  mA  is  determined  from  the  conversion  factor  of  2500:1000 
as  specified  in  Appendix  A. 
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a  CD.  The  menu  screen  provides  the  current  date  and  time  along  with  the  amount  of  data  that 
has  not  yet  been  written  to  a  CD  and  the  number  of  CDs  needed  to  do  so.  Should  the  amount 
of  storage  space  on  the  hard  drive  run  out,  the  application  will  begin  erasing  snapshots  starting 
with  the  oldest  to  write  newer  data.  Once  snapshot  data  has  been  written  to  a  CD,  it  is  erased 
from  the  hard  drive. 


Figure  3.7:  NILM  Menu  Screen 


Manual  collection  of  data  is  also  possible  by  first  exiting  the  menu  system  by  going  to  the 
option  labeled  “Shell  prompt”.  This  places  the  user  in  the  nilmdata  directory.  To  create  a 
new  directory  in  which  to  place  the  test  data,  type  mkdir  directory  name  and  then  change  to 
that  directory  by  typing  cd  directory  name .  To  collect  raw  voltage  and  current  data  in  ASCI2 
format,  type  dd  if=/dev/pcii710  |  pci2asc  >  voltage/current  text  file.  This  will  result 
in  voltage  and  current  being  recorded  as  two  columns  of  text  in  the  file  called  voltage/current 
textfile. 

To  collect  preprocessed  power  data,  type  prep-flush  /dev/pci  1710  >  power  file.  This 
will  result  in  the  fundamental  and  first  three  odd  harmonics  of  real  and  reactive  power  being 
recorded  as  eight  columns  of  text  in  the  file  called  power  file.  ; 

It  is  also  possible  to  perform  both  of  the  preceding  operations  simultaneously  by  using  the  tee 
command.  To  collect  both  voltage  and  current  data  along  with  preprocessed  power  data,  type  dd 
if=/dev/pcil710  [  tee  voltage/current  binaryfile  |  prep-flush  -stdin  -p  >  power  file. 
Unfortunately  this  provides  the  voltage  and  current  data  in  a  binary  format,  so  to  convert  to  a 
more  useful  text  file  type  more  voltage/current  binaryfile  |  pci2asc  >  voltage/current 
textfile. 

In  order  to  be  able  to  plot  the  data  in  a  program  like  MATLAB  or  Octave,  open  voltage/current 
textfile  in  a  text  editor  and  delete  any  non-numerical  information  which  may  occur  at  the  top 
or  bottom  of  the  data. 

To  manually  write  data  to  a  CD,  change  back  to  the  nilmdata  directory  by  typing  cd  . . 
and  then  creating  an  iso  file  of  the  directory  by  typing  mkisofs  -o  directory  name,  iso  -r  -J 
directory  name .  Once  the  iso  file  is  created,  write  to  the  CD  by  placing  a  blank  CD  in  the  drive 
and  typing  cdrecord  -v  -eject  dev=0,0,0  speed=12  directory  name. iso , 
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Chapter  4 

Cycling  System  Results 


4.1  Sewage  System  Onboard  USCGC  SENECA 

The  sewage  system  onboard  the  USCGC  SENECA  was  monitored  because  is  represents  an  op¬ 
portunity  to  perform  trending  analysis  on  a  cycling  system.  This  cycling  system  is  especially 
intriguing  because  its  operation  is  affected  not  only  by  physical  parameters  (e.g.,  the  rate  at 
which  the  system’s  two  vacuum  pumps  raise  vacuum)  but  also  by  human  behavior.  Monitoring 
the  system’s  electrical  loads  and  trending  their  power  usage  could  provide  insight  into  the  system’s 
integrity  or  provide  an  early  indicator  of  changes  in  the  crew’s  health. 

The  system  consists  of  commodes,  urinals,  and  drains  which  discharge  into  a  vacuum  collection 
tank.  Vacuum  is  maintained  in  this  tank  by  two  pumps  which  alternately  cycle  to  maintain 
vacuum  between  approximately  15  and  10  inches  of  mercury  (in.  Hg).  Two  transfer  pumps  cycle 
on  and  off  to  occasionally  move  waste  from  the  vacuum  collection  tank  to  a  larger  sanitary  tank. 
Figure  4.1  shows  the  system’s  vacuum  pumps  with  the  vacuum  collection  tank  in  the  background. 


Figure  4.1:  Sewage  System  Vacuum  Pumps 
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4.2  Method  of  Data  Collection  and  Processing 

There  is  a  complete  NILM  system  (as  discussed  in  Chapter  3)  monitoring  the  two  vacuum  pumps 
and  two  transfer  pumps  that  comprise  the  electrical  loads  of  the  SENECA’s  sewage  sytem.  One  LA 
205-S  current  transducer  is  monitoring  all  four  electrical  loads.  The  resistors  used  for  the  voltage 
measurement,  current  measurement,  and  references  are  130  0,  66.5  Q,  and  56.2  O,  respectively. 
Upon  initial  system  installation,  the  system  was  calibrated  by  collecting  balanced  three-phase 
power  data  with  a  Fluke  39  Powermeter  and  comparing  it  to  the  power  counts  obtained  with  the 
NILM.  The  resultant  scaling  factor  between  the  measured  power  (1.42  kW)  and  counts  (-2287.2) 

"^-0.61946  Si- 

Data  was  collected  in  one  hour  “snapshots”  in  order  to  trend  the  cycling  of  electrical  loads 
(primarily  the  vacuum  pumps),  A  typical  hour  of  power  data  is  shown  in  Figure  4.2.  An  expanded 
view  of  the  first  pump  run  in  this  typical  hour  is  shown  in  Figure  4.3.  Figure  4.2  shows  that  the 
electrical  power  drawn  by  each  pump  is  approximately  equal. 
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Figure  4.2:  Vacuum  Pump  Transients  for  One  Hour  of  Weekend  Use 

Because  of  the  relatively  large  power  transient  associated  with  startup  and  the  subsequent 
variations  in  steady  state  power,  pump  runs  like  those  shown  in  Figure  4.3  are  difficult  to  trend 
with  respect  to  time.  For  this  reason,  a  MATLAB  script  entitled  pumprun  was  written  to  accept 
an  hour  of  sewage  system  “snapshot”  data  as  input  and  provide  the  number  of  pump  runs  and 
average  run  time  per  pump  run  as  outputs.  The  characteristics  of  this  script  are  as  follows  (see 
Appendix  C  for  the  script  itself): 

~  •  Unsealed  power  data  (i.e.  in  terms  of  counts,  not  W  or  kW)  is  processed  into  binary  form 
by  assigning  a  value  of  “1”  to  data  points  of  greater  than  500  counts  and  a  value  of  “0”  to 
data  points  of  less  than  or  equal  to  500  counts.  The  binary  version  of  Figure  4.2  is  shown 
in  Figure  4.4. 

•  The  binary  power  data  is  examined  for  any  changes  from  the  off  (0)  state  to  the  on  (1) 
state.  Anytime  in  the  binary  data  that  a  0  is  followed  by  a  1,  the  run  counter  is  increased 
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Figure  4.3:  Vacuum  Pump  Start 


Figure  4.4:  Binary  Form  of  Vacuum  Pump  Transients 


by  1.  If  the  first  element  of  the  binary  power  data  is  a  1,  then  the  run  counter  starts  with  a 
value  of  1  instead  of  0  since  a  pump  was  already  running  when  data  collection  commenced. 

The  binary  power  is  summed  and  then  divided  by  120  (the  sampling  frequency  in  Hz  for 
the  original  power  data)  and  the  number  of  pump  runs  to  obtain  the  average  run  time  per 
pump  run  in  seconds. 


4.3  Sewage  System  Simulation 

The  fundamental  problem  with  performing  controlled  tests  on  a  system  like  the  sewage  system 
onboard  the  SENECA  is  control  of  the  human  element.  Although  it  is  possible  to  insert  vacuum 
leaks  of  various  quantifiable  sizes  into  the  system  and  then  trend  the  resulting  system  cycles, 
it  is  not  possible  to  control  the  main  factor  which  effects  the  amount  of  system  run  time,  the 
amount  of  usage  by  the  ship’s  crew.  In  order  to  determine  effects  of  vacuum  leaks  and  crew  usage 
on  system  run  times  independently,  the  MATLAB  script  flushsim  was  written  to  simulate  the 
system  behavior.  The  simulation  accepts  as  input  the  average  time  between  system  usage  events 
and  vacuum  leak  rate  (in  inches  of  Hg  per  minute)  and  provides  the  number  of  pump  runs  and 
average  run  time  per  pump  run  as  outputs.  The  characteristics  of  the  simulation,  which  can  be 
found  in  Appendix  C,  are  as  follows: 

•  Although  not  an  explicit  input,  the  number  of  minutes  for  which  the  simulation  runs  along 
with  the  number  of  calculations  used  per  minutes  can  be  varied  by  changing  the  parameters 
“n”  and  “k” ,  respectively. 

•  The  simulation  is  written  from  the  standpoint  of  pressure,  with  pressure  maintained  between 
15  and  10  inches  (Hg)  of  vacuum.  Whenever  vacuum  falls  to  less  that  10  in.,  a  pump  run 
indicator  is  switched  to  the  “on”  state.  When  vacuum  exceeds  15  in.,  the  pump  run  indicator 
switches  back  to  the  “off”  state.  If  vacuum  should  fall  to  less  than  7.5  in.  due  to  a  large  leak 
rate  or  heavy  system  usage,  both  pumps  are  turned  on  (or  in  the  case  of  a  pump  already 
running,  the  other  pump  is  turned  on). 

•  To  increase  realism,  the  two  pumps  in  the  simulation  have  different  vacuum  drawing  rates 
(19  and  17  in.  per  minute). 

•  A  system  usage  event,  henceforth  referred  to  as  a  “flush”,  is  modelled  as  a  discrete  vacuum 
drop  of  1  in.  The  probability  of  a  flush  is  calculated  using  the  triangular  probability  density 
function  (pdf)  shown  in  Figure  4.5.  A  triangular  probability  distribution  was  chosen  because 
it  was  deemed  complex  enough  to  add  realism  to  the  simulation  while  being  simple  enough 
to  allow  for  speed  of  calculation.  This  pdf  is  based  on  the  average  time  between  flushes 
entered  as  input  (tmean)  and  assumes  that  if  the  probability  of  a  flush  occurring  within 
tmean  is  0.5,  then  the  probability  of  a  system  usage  event  occurring  within  2tmean  is  one.  At 
each  point  in  the  simulation,  the  probability  of  a  flush  occurring  is  calculated  (based  on  the 
amount  of  time  which  has  elapsed  since  the  last  flush)  from  the  assumed  pdf  and  compared 
to  a  randomly  generated  number.  If  the  random  number  is  less  than  the  probability  a  flush 
occurs,  vacuum  drops  by  1  in.,  and  the  flush  counter  is  reset.  If  the  random  number  is 
greater  than  or  equal  to  the  probability  no  vacuum  drop  occurs  and  the  probability  of  a 
flush  will  be  slightly  higher. 

•  The  number  of  pump  runs  and  average  run  time  per  run  is  calculated  from  the  pump  run 
indicator  vector  in  the  same  manner  as  for  real  data. 

A  sewage  system  simulation  run  showing  vacuum  Information  and  the  pump  status  indicator 
are  shown  versus  time  in  Figure  4.6.  The  simulation  was  run  for  an  average  of  four  minutes 
between  flushes  with  no  vacuum  leak. 
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Figure  4.5:  Vacuum  System  Usage  Probability  Density  Function 


Figure  4.6:  One  Hour  of  Simulated  Vacuum  and  Vacuum  Pump  Cycling 


4.4  Trending  of  System  Cycling  for  Various  Vacuum 
Leaks 

In  order  to  trend  the  cycling  of  the  vacuum  pumps  with  a  fault  in  the  sewage  system,  vacuum 
leaks  of  varying  sizes  were  inserted  into  the  sewage  system  onboard  the  SENECA.  The  vacuum 
leak  was  controlled  (and  quantified)  by  a  flowmeter  attached  to  a  gauge  line  connection  just  off 
the  vacuum  collection  tank.  The  flowmeter  is  shown  in  Figure  4.7. 

&*'  The  throttle  valve  on  the  flowmeter  was  adjusted  to  achieve  the  desired  flow  rate  and  NILM 
snapshot  data  was  collected.  The  flowmeter  contains  two  different  balls,  one  made  of  stainless 
steel  and  one  made  of  plastic.  For  consistency,  leak  rates  specified  throughout  this  chapter  are 
the  leak  rates  read  off  of  the  plastic  ball  immediately  after  adjusting  the  flowmeter  to  a  new  leak 
rate.  Leak  data  is  contained  in  Appendix  D.  A  snapshot  of  leak  data  is  shown  in  Figure  4.8. 

There  are  two  differences  between  Figure  4.8  and  Figure  4.2  shown  before,  Figure  4.8  shows 
cycling  of  vacuum  pumps  on  a  weekday  with  a  leak  of  62.5  standard  cubic  feet  per  hour  (SCFH) 
present  while  Figure  4.2  shows  cycling  of  vacuum  pumps  on  a  weekend  day  with  no  vacuum  leak 
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Table  4.1:  Comparison  of  Two  Different  Hours  of  Data  from 
SENECA 


Leak  Rate  (SCFH) 

0  (Weekend) 

62.5  (Weekday) 

Vacuum  Pump  Runs 

7 

16 

Average  Run  Time  per  Run  (seconds) 

17.89 

20.77 

Total  Run  Time  (seconds) 

125.23 

332.32 

It  can  be  seen  from  Table  4.1  that  the  total  run  time  for  the  weekday  is  higher  than  for  the 
weekend  day.  But  is  the  higher  total  run  time  due  to  heavier  usage  or  the  presence  of  a  vacuum 
leak?  Secondly,  and  perhaps  of  equal  importance,  does  the  increased  total  run  time  manifest  itself 
as  an  increase  in  the  number  of  pump  runs  per  hour,  the  average  run  time  per  pump  run,  or  a 
combination  of  the  two? 

4.4.1  Dependence  of  Total  Run  Time 

In  order  to  determine  the  variation  of  total  run  time  on  different  factors,  the  sewage  system 
simulation  was  used  to  create  a  large  amount  of  data.  The  simulation  (as  discussed  previously) 
was  run  for  combinations  of  13  different  average  flush  times  and  11  different  leak  rates.  Each  one 
of  these  combinations  was  run  24  times  and  averaged.  The  result  of  these  simulation  runs  was 
three  different  13  by  11  arrays.  Each  element  of  the  arrays  represents  total  run  time,  number  of 
pump  runs,  or  average  run  time  per  run  for  one  hour  of  vacuum  pump  cycling. 

Since  the  ultimate  goal  of  trending  the  sewage  system  data  is  to  be  able  to  detect  a  vacuum 
leak  in  the  presence  of  system  usage  variations,  scattergrams  were  produced  which  plotted  total 
run  time  vs.  either  number  of  runs  or  average  run  time  per  run  for  different  levels  of  usage  with 
a  constant  leak  rate.  The  scattergrams  are  shown  in  Figures  4.9  and  4.10. 


Figure  4.9:  Scattergram  of  Total  Runtime  vs.  Number  of  Runs  for  Three  Simulated  Leak 
Conditions 
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Figure  4,10:  Scattergram  of  Total  Runtime  vs.  Average  Run  Time  per  Run  for  Three 
Simulated  Leak  Conditions 


Figure  4.9  indicates  that  total  run  time  seems  to  be  related  to  number  of  pump  runs  in  a 
more  linear  fashion  than  the  average  run  time  per  run.  For  higher  leak  rates,  the  total  run  time 
appears  to  only  shift  along  the  same  line  to  higher  values.  It  is  for  these  reasons  that  the  number 
of  pump  runs  was  focused  on  in  the  analysis  of  real  leak  data. 

4.4.2  Development  of  the  Leak  Indicator 

The  following  steps  (which  will  be  discussed  in  more  depth  below)  summarize  the  development  of 
the  leak  indicator: 

1.  Snapshot  data  for  several  hours  of  data  is  uncompressed  and  run  through  the  pumprun 
Matlab  script  discussed  previously  to  determine  the  number  of  pump  runs  and  total  run 
time  for  that  hour  of  data.  This  effectively  distills  an  hour  of  power  data  into  only  two 
numbers. 

2.  The  number  of  pump  runs  and  total  run  time  is  tabulated  for  all  hours  of  data  collected. 

3.  The  average  number  of  runs  per  hour  and  standard  deviation  in  number  of  runs  per  hour 
is  calculated.  The  average  used  was  the  arithmetic  mean. 

4.  A  linear  regression  is  performed  on  total  run  time  vs.  number  of  pump  runs  for  the  tabulated 
hourly  data. 

5.  The  leak  “indicator”  is  the  standard  deviation  in  pump  runs  and  the  y-intercept  from  the 
linear  regression  subtracted  from  the  average  number  of  runs  per  hour. 

The  leak  indicator  starts  with  the  average  number  of  pump  runs  per  hour  and  standard 
deviation  in  number  of  pump  runs  as  calculatated  with  equations  4.1  and  4.2.  Runs  is  the 
average  number  of  pump  runs,  n  is  the  number  of  hours  under  consideration,  RunSi  is  the  number 
of  runs  for  the  ith  hour,  and  SRuns  is  the  standard  deviation  in  number  of  pump  runs. 
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Since  it  has  already  been  shown  that  the  amount  of  pump  run  time  is  directly  related  to 
number  of  pump  runs,  it  was  felt  that  the  average  number  of  pump  runs  was  a  logical  starting 
point  for  analysis.  After  calculating  the  statistics  (Runs  and  SRuns )  for  the  number  of  pump  runs, 
it  became  clear  that  the  average  number  of  pump  runs  could  not  be  used  alone  as  an  indicator  of 
a  vacuum  leak.  The  statistics  for  the  number  of  pump  runs  are  summarized  in  Table  4.2. 


Table  4.2:  SENECA  Vacuum  Leak  Statistics 


Leak  Rate  (SCFH) 

Runs 

SRuns 

132.5 

21.57142857 

4.878749615 

102.5 

11.34782609 

2.207604536 

62.5 

14.52173913 

4.429955522 

51 

14.2173913 

5.767909041 

42 

12.40909091 

5.160833182 

0 

9.136363636 

2.623981971 

As  seen  in  Table  4.2,  the  average  number  of  runs  for  the  leak  rate  of  102.5  SCFH  does  not 
follow  the  trend  evident  in  the  rest  of  the  table.  It  should  also  be  noted  however,  that  the  standard 
deviation  in  the  number  of  pump  runs  for  the  0  and  102.5  SCFH  leak  rates  is  much  lower  than  the 
other  leak  rates.  The  previous  observations  are  shown  in  Figure  4.11  with  the  pump  run  statistics 
used  to  plot  normal  distributions  (see  [13])  in  accordance  with  equation  4.3.  It  should  be  noted 
that  plotting  the  pump  run  statistics  as  a  normal  distribution  is  for  illustrative  purposes  only  since 
there  were  no  tests  to  determine  if  the  pump  run  data  actually  followed  a  normal  distribution. 
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The  difference  in  both  the  average  number  of  pump  runs  and  standard  deviation  in  pump 
runs  is  due  to  the  day  of  the  week  that  the  data  was  taken. 

Since  the  ultimate  goal  of  trending  the  leak  data  was  to  develop  a  parameter  which  detected 
the  presence  of  a  vacuum  leak  regardless  of  the  amount  or  variance  in  system  usage  by  the  crew, 
the  standard  deviation  was  subtracted  from  the  average  number  of  pump  runs  to  negate  the 
difference  between  weekend  days  and  weekdays.  Since  a  vacuum  leak  is  a  continuous  drain  on  the 
system  and  therefore  does  not  affect  the  system  differently  from  one  hour  to  the  next,  performing 
this  subtraction  removes  a  value  which  is  only  affected  by  crew  usage.  The  standard  deviation  in 
..number  of  pump  runs  subtracted  from  the  average  number  of  pump  runs  is  dubbed  the  “pump 
run  lower  bound”. 
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Figure  4.11:  Normal  Distributions  for  Number  of  Pump  Runs  at  Different  Leak  Rates 


The  next  step  of  developing  a  leak  indicator  was  to  plot  the  total  run  time  vs.  number  of 
pump  runs  as  shown  in  Figure  4.12. 
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Figure  4.12:  Scattergram  of  Total  Runtime  vs.  Number  of  Runs  for  Three  Real  Leak 
Conditions 


Even  though  the  data  shows  more  scatter  than  in  the  simulated  case1,  it  still  shows  a  generally 
linear  trend.  A  least  squares  linear  regression  was  therefore  performed  in  accordance  with  [13]  for 


"‘  This  could  be  a  result  of  the  simulation  being  too  simplistic  or  unable  to  account  for  all  of  the  com¬ 
plexities  of  the  real  system  \ 


40 


total  run  time  vs.  number  of  runs  to  determine  how  the  slope  and  y-intercept  varied  with  leak 
rate. 

The  following  is  an  outline  of  the  linear  regression  process: 

1.  It  is  desired  to  fit  the  data  to  an  equation  of  the  form 

V  =  0 ix  +  0o  (4.4) 

So  that  the  total  amount  of  deviation  between  the  regression  line  and  the  observed  points 
as  given  in  equation  4.5  is  minimized. 


s  (Po,  Pi)  =^2[yi-  (Po  +  PiXi)Y 

i=  1 


(4.5) 


2.  To  minimize  the  value  of  S  with  respect  to  both  p0  and  Pi,  the  partial  derivative  of  S  with 
respect  to  Po  and  Pi  is  set  to  zero. 
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3.  Rearranging  results  in  two  equations  in  terms  of  /3q  and  /?x. 


n  n 


Pon  + 0i  Y;xi  = 

i= 1  i= 1 


n  n 


00^2  xi  +  0 lJ2Xi  XiVi 

i= 1  i= 1  i=l 

4.  This  results  in  the  following  expressions  for  0q  and  0\ . 


(4.6) 

(4.7) 
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Po  =  y-  0ix 

The  results  of  these  linear  regressions  are  shown  in  Table  4.3. 


(4.10) 


(4.11) 
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Table  4.3  shows  intercepts  that  are  all  non-zero.  Since  the  y-axis  in  the  linear  regression  is 
number  of  pump  runs,  the  interpretation  of  non-zero  intercepts  is  that  for  zero  total  run  time, 
there  will  be  more  than  zero  pump  runs.  There  are  two  possible  explanations  for  this  phenomenon. 
Either  there  is  a  leak  in  the  system  other  than  the  one  placed  in  it  for  testing  or  the  behavior  of 
the  system  for  low  total  run  times  is  non-linear. 

Since  in  reality  there  would  be  zero  pump  runs  for  zero  total  run  time,  the  y-axis  should  be 
shifted  upwards  by  a  correction  factor  equal  to  the  y-intercept  for  each  of  the  data  sets  considered. 
So  the  final  leak  indicator  is  the  pump  run  lower  bound  discussed  above  minus  the  y-intercept 
from  the  linear  regression  of  total  run  time  vs,  number  of  pump  runs  as  shown  in  equation  4.12. 

Leak  Indicator  =  Runs  —  SRuns  —  Pq  (4-12) 

This  final  number  is  an  indicator  of  amount  of  vacuum  pump  cycling  above  the  no-cycling 
level.  The  leak  indicator  is  plotted  vs.  leak  rate  in  Figure  4.13  for  the  leak  data  taken  aboard  the 
SENECA. 


Figure  4.13:  Leak  Rate  vs.  Leak  “Indicator” 
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4.4.3  Validation  of  the  Leak  Indicator 


In  order  to  verify  the  ability  of  the  the  leak  indicator  to  show  the  presence  of  a  leak  regardless 
of  the  amount  of  system  usage,  more  simulation  data  were  created.  Weekend  and  weekday  data 
were  simulated  for  24  hours  for  several  different  leak  rates  and  the  leak  indicator  calculated.  The 
weekday  simulations  had  a  lower  average  time  between  flushes  with  a  higher  variation  in  time 
between  flushes  while  the  weekend  simulations  had  a  higher  average  time  between  flushes  with  a 
lower  variation  in  time  between  flushes  (see  Appendix  D).  The  weekday  and  weekend  results  are 
shown  in  Table  4.4  and  Figure  4.14. 


Table  4.4:  Leak  Indicator  Summary 


Leak  Rate 
(in.  Hg/min.) 

Weekday  Leak 
“Indicator” 

Weekend  Leak 
“Indicator” 

0 

6.296625689 

5.903347459 

0.5 

9.973276855 

11.47976765 

1 

14.21729058 

17.16598438 

1.5 

17.33927728 

18.96290688 

2 

17.48596088 

17.69171294 

2.5 

19.17635935 

18.64716683 

Leak  Rate  (in  Hg/mfn.) 


-Weekday  —■—Weekend  | 


Figure  4.14:  Leak  Rate  vs.  Leak  “Indicator”  for  Simulation  Data 


It  is  encouraging  to  note  from  Figure  4.14  that  as  with  the  leak  indicators  developed  for  the 
actual  shipboard  data,  there  is  a  fairly  linear  relationship  between  leak  indicator  and  leak  rate 
independent  of  system  usage.  It  is  also  important  to  note  that  the  leak  indicators  for  the  weekend 
and  weekday  cases  are  very  similar  for  the  two  no  leak  cases. 
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4.5  Chapter  Summary 

The  results  presented  previously  in  this  chapter  show  the  existence  of  a  leak  indicator  for  cycling 
systems  that  can  be  used  to  detect  a  fault  in  the  system.  This  leak  indicator  uses  the  statistics 
of  pump  runs  to  arrive  at  a  figure  of  merit  which  does  not  depend  on  the  amount  or  variation  in 
system  usage. 

The  leak  indicator  was  developed  from  a  purely  empirical  standpoint  based  on  observations 
made  from  collected  data  and  physical  reasoning.  The  best  way  of  determining  its  ability  to 
detect  faults  is  to  take  data  on  others  system  that  involve  cycling  components 


.Chapter  5 

Fluid  System  Results 


5.1  ASW  System  Onboard  USCGC  SENECA 

The  auxiliary  seawater  (ASW)  system  onboard  the  USCGC  SENECA  provides  cooling  for  all 
heat  loads  onboard  the  ship  other  than  those  associated  with  the  main  engines.  Heat  loads 
cooled  include  the  HVAC  units,  refrigerators  and  freezers,  diesel  generator  air  coolers,  and  diesel 
generator  lube  oil  coolers.  Suction  to  the  two  ASW  pumps  is  taken  through  two  supply  lines 
from  the  bottom  of  the  ship  through  suction  screens,  isolation  globe  valves,  suction  strainers,  and 
butterfly  valves.  The  isolation  globe  valves  determine  which  suction  is  providing  ASW  cooling 
while  the  butterfly  valves  allow  isolation  of  the  strainers  for  cleaning.  A  butterfly  valve  across 
the  inlets  also  allows  cross-connecting  the  system  to  provide  redundancy.  Changing  the  flow  rate 
is  accomplished  by  throttling  two  butterfly  valves  at  the  outlets  of  the  two  pumps:  a  pump  is 
started  with  the  overboard  throttle  valve  completely  open  and  the  heat  load  throttle  valve  shut, 
then  the  heat  load  throttling  valve  is  opened  and  the  overboard  throttle  valve  is  partially  shut. 
The  overboard  valve  is  throttled  until  a  pressure  gauge  in  the  supply  line  for  the  heat  loads  shows 
a  pressure  of  approximately  35  psig.  Figure  5.1  shows  a  one-line  diagram  of  the  system  while 
Figure  5.2  shows  some  of  the  components  at  the  inlets  to  the  pumps. 
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Figure  5.1:  ASW  System  Schematic 


Since  seawater  cooling  systems  are  common  across  a  wide  range  of  ships,  monitoring  the 
seawater  cooling  system  on  a  ship  with  many  heat  loads  was  seen  as  an  opportunity  to  examine 
how  the  power  usage  of  the  pumps  that  maintain  flow  in  the  system  varies  with  changes  in  the  flow 
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Figure  5.2:  ASW  Inlet  Isolation  Valves  and  Strainers 


characteristics  upstream  and  downstream  of  the  pumps.  Specifically,  it  was  desired  to  know  how 
pump  power  would  change  for  differing  amounts  of  flow  that  was  sent  to  heat  loads,  clogging  in  the 
inlet  strainers,  sudden  rapid  failure  of  the  coupling  beween  motor  and  pump,  and  flow  blockage 
that  might  occur  in  any  of  the  heat  loads.  The  last  item  was  especially  important  because  several 
crewmembers  of  the  USCGC  SENECA  had  made  reference  to  mussels  growing  inside  the  endbells 
of  tube-and-shell  heat  exchangers  in  the  diesel  generators  and  the  HVAC  units.  It  was  therefore 
deemed  especially  important  to  focus  attention  on  this  problem  to  provide  a  possible  prediction 
of  heat  exchanger  fouling  well  in  advance  of  a  rise  in  temperature. 

The  NILM  installation  for  the  ASW  system  uses  one  voltage  transducer  and  two  current 
transducers  to  monitor  the  two  seawater  pumps  that  pressurize  the  supply  header.  Since  there  is 
no  one  panel  that  supplies  both  of  the  pumps  (for  reliability),  one  LA  305-S  current  transducer 
had  to  be  placed  in  each  pump  controller  and  the  output  combined  at  the  NILM  box.  The  resistors 
used  for  the  voltage  measurement,  current  measurement,  and  references  are  130  Q,  7.5  O,  and 

56.2  S7,  respectively.  Upon  initial  system  installation,  the  system  was  calibrated  by  collecting 
balanced  three-phase  power  data  with  a  Fluke  39  Powermeter  and  comparing  it  to  the  power 
counts  obtained  with  the  NILM.  The  resultant  scaling  factor  between  real  power  (29.5  kW)  and 
counts  (4150.8)  was  7.11 

5.2  Fluid  Test  System 

The  ASW  system  onboard  the  SENECA  is  very  promising  as  a  system  on  which  to  conduct 
research  into  the  relationships  between  fluid  dynamics  and  electrical  power  usage,  but  there  are 
a  few  drawbacks: 

•  Availability,  The  SENECA  is  only  available  for  short  periods  of  time  while  in  port.  This 
coupled  with  the  very  aggressive  maintenance  and  repair  schedule  while  in  port  means  that 
crewmembers  have  very  little  time  in  their  busy  schedule  to  assist  with  data  collection. 

•  Complexity.  Since  the  ASW  system  is  sized  to  cool  many  different  heat  loads  all  around 
the  ship,  the  system  contains  a  large  number  of  components  as  well  as  a  large  amount  of 
piping.  As  a  result,  there  are  complex  interactions  in  the  system  which  can  be  difficult  to 
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model  or  predict.  Also,  the  effect  of  a  small  change  in  one  of  the  many  components  may 
not  be  significant  enough  to  be  noticed  in  the  power  signature  of  the  pumps 

•  Controllability.  In  order  to  predict  the  response  of  the  ASW  system  to  specific  faults, 
it  would  be  necessary  to  take  test  data  with  those  faults  inserted  into  the  system.  For 
example,  to  see  how  mussels  growing  inside  of  a  diesel  generator  lube  oil  cooler  effects  the 
y.  power  signature  of  the  two  pumps,  it  would  be  necessary  to  clog  tubes  in  the  heat  exchanger, 

f  "  Since  this  is  impractical,  the  clog  can  only  be  simulated  by  throttling  the  inlet  valve  to  the 

heat  exchanger.  Although  it  might  be  possible  to  create  the  same  pressure  drop  across  the 
heat  exchanger  by  doing  this,  it  is  does  not  necessarily  reproduce  the  more  complicated 
hydrodynamic  phenomena  (or  the  same  phenomena)  that  occur  when  tubes  are  actually 
blocked. 

Because  of  the  limitations  stated  above,  a  fluid  test  system  was  constructed  to  collect  data  in 
the  laboratory.  Such  a  system  has  unlimited  availability,  a  relatively  low  number  of  components, 
and  a  much  higher  level  of  controllability  (in  the  sense  that  it  is  easier  to  change  the  fundamental 
characteristics  of  the  components  that  comprise  it) . 

The  fluid  test  system  that  was  constructed  consists  of  a  three-phase  0.5  hp  pump  that  takes 
a  suction  on  a  large  water  reservoir,  pumps  the  water  through  some  combination  of  three  parallel 
flow  paths,  then  returns  the  water  to  the  water  reservoir.  A  one-line  diagram  of  the  fluid  test 
system  is  shown  in  Figure  5.3  while  Figure  5.4  provides  a  front  view. 


Figure  5.3:  Fluid  Test  System  Schematic 

The  pump  suction  run  consists  of  a  0.75  in.  globe  valve  attached  to  1  in.  PVC  piping.  A 
small  ball  valve  and  plastic  tubing  are  provided  to  easily  drain  the  system  when  not  in  use.  The 
pump  outlet  piping  consists  of  three  different  branches,  each  consisting  of  .75  in.  copper  tubing 
with  a  .75  in.  globe  to  provide  isolation  and  throttling. 

•'  Two  of  the  branches  contain  “heat  loads”  while  the  third  branch  returns  directly  to  the  water 
reservoir.  One  of  the  two  “heat  load”  branches  is  a  section  of  pipe  which  makes  seven  180  degree 
turns  to  mimic  the  construction  of  a  water-to-air  heat  exchanger  as  one  might  find  in  the  generator 
section  of  a  diesel  generator. 

The  other  “heat  load”  branch  is  constructed  as  a  tube-and-shell  heat  exchanger,  as  one  might 
use  to  cool  lubricating  oil  in  a  diesel  engine  or  refrigerant  in  an  HVAC  plant.  The  heat  exchanger 
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Figure  5.4:  Fluid  Test  System  Front  View 


(henceforth  abbreviated  HX)  consists  of  seven  0,5  in.  PVC  tubes  attached  to  tubesheets  at  both 
ends,  mounted  inside  of  a  4  in.  PVC  tube.  Figure  5.5  shows  the  HX  prior  to  assembly. 


Figure  5.5:  Fluid  Test  System  HX  Prior  to  Assembly 


Water  in  the  test  system  flows  through  .75  in.  copper  piping,  into  the  4  in.  inlet  HX  endbell, 
through  the  inside  of  seeven  0.5  in.  PVC  tubes,  into  the  4  inch  outlet  HX  endbell,  and  finally 
through  a  run  of  .75  in.  copper  piping  back  to  the  water  reservoir.  This  flowpath  represents  the 
flow  of  the  cooling  medium  (such  as  seawater)  through  the  HX.  Since  the  phenomena  of  interest 
are  hydrodynamic  and  not  thermodynamic,  the  HX  was  not  designed  to  actually  transfer  any 
heat.  If  it  had  been  designed  to  do  so,  there  would  be  a  second  flowpath  in  which  the  fluid  being 
cooled  would  flow  around  the  outside  of  the  small  tubes,  which  would  be  made  of  copper  (or  some 
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other  metal  to  provide  a  high  coefficient  of  heat  transfer)  and  not  PVC. 

Electrically,  the  pump  power  and  NILM  monitoring  are  combined  in  a  single  NEMA  enclosure. 
The  enclosure  contains  the  components  for  a  NILM  box  discussed  in  Chapter  3  as  well  as  a  three- 
phase  switch  and  an  LA  55-P  current  transducer.  The  resistors  used  for  the  voltage  measurement, 
current  measurement,  and  references  are  130  Q,  66.5  fi,  and  56.2  C,  respectively.  Upon  initial 
system  construction,  the  system  was  calibrated  by  collecting  balanced  3-phase  power  data  with 
Jr  Fluke  39  Powermeter  and  comparing  it  to  the  power  counts  obtained  with  the  NILM.  The 
'resultant  scaling  factor  between  real  power  (152  W)  and  counts  (-440.58)  was  -0.345  X~lnt  • 

5.3  Building  a  Centrifugal  Pump  Simulation 

A  simulation  of  a  centrifugal  pump  startup  was  developed  for  the  following  reasons: 

•  Simulations  based  on  physical  principles  offer  a  means  of  validating  understanding  of  the 
system  upon  which  they  are  based. 

•  They  offer  the  ability  to  quickly  collect  large  amounts  of  simulated  data  while  varying  many 
different  parameters,  often  in  an  automated  fashion. 

•  They  provide  an  approximation  of  system  behavior  for  situations  which  are  not  possible  or 
feasible  to  reproduce  in  the  laboratory  (especially  catastrophic  failures). 

The  simulation  that  was  constructed  takes  known  parameters  for  an  induction  motor  and 
couples  them  to  a  fictitious  pump  impeller  which  is  geometrically  similar  to  the  pumps  onboard 
the  SENECA.  The  motor  model  from  [15]  is  a  fifth-order  model  of  an  induction  motor  with  a 
rated  power  of  3  hp.  Since  pump  curves  for  the  ASW  pumps  onboard  the  SENECA  were  available 
in  [16],  it  was  possible  to  scale  the  pump  impeller  to  the  much  smaller  motor.  The  pumps  onboard 
the  SENECA  are  powered  by  40  hp  induction  motors,  but  the  operating  point  specified  in  [16] 
corresponds  to  a  brake  horsepower  (BHP)  of  only  31  hp.  Using  the  same  ratio  of  motor  power  to 
pump  BHP,  it  was  decided  that  the  scaled  pump  head  would  therefore  have  a  brake  horsepower 
of  2.25  hp.  The  ASW  pump  motors  also  have  two  poles  versus  the  simulation  motor’s  four,  so  the 
simulation  pump  was  assigned  an  operating  speed  of  1775  revolutions  per  minute  (rpm)  versus 
the  ASW  pump’s  3550  rpm. 


5.3.1  Inertial  Considerations 


The  first  consideration  when  adding  a  pump  impeller  to  an  electric  motor  is  the  additional  moment 
of  inertia  added  by  the  pump  impeller  itself  as  well  as  the  moment  of  inertia  added  by  the  water 
being  pumped.  The  first  step  in  estimating  these  two  moments  of  inertia  was  to  determine  the 
pump  impeller  diameter  through  geometrical  scaling.  The  diameter  of  the  pump  impeller  was 
scaled  using  the  power  coefficient  given  in  [1]  by  equation  5.1. 


Cp  — 


P 

pu>3D 5 


(5.1) 


Maintaining  this  coefficient  constant  for  the  same  pumped  fluid  leads  to  an  expression  for 
-the  diameter  of  the  scaled  pump  impeller  as  shown  in  equation  5.2.  The  calculated  value  of  the 
simulated  pump  diameter  was  5.30  in. 
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I>2  =  D\ 


(5.2) 


Given  the  diameter  of  the  pump  impeller,  a  three-dimensional  model  was  created  in  Rhino 
3D©  software1.  This  program  was  used  because  its  built  in  analysis  functions  allowed  easy 
calculation  of  the  impeller’s  moment  of  inertia.  The  pump  impeller  from  the  fluid  test  system 
described  earlier  was  used  as  a  model  for  the  geometry  of  the  simulated  impeller.  The  only 
difference  was  that  the  simulated  pump  impeller  is  designed  with  an  additional  plate  on  the  inlet 
side  of  the  rotor  to  make  it  a  “closed  rotor”  design  instead  of  the  “open  rotor”  design  used  in 
the  fluid  test  system.  It  was  felt  that  this  was  much  closer  to  the  style  of  impeller  used  on  larger 
pumps  like  those  found  in  the  ASW  system  of  the  SENECA.  In  addition  to  having  the  scaled 
diameter  determined  previously,  the  front  and  back  plates  are  .25  in.  thick  with  vanes  .5  in.  thick. 
There  is  a  .5  in.  diameter  hole  in  the  back  plate  for  the  motor  shaft  and  a  1  in.  diameter  hole 
in  the  front  plate  as  an  inlet  into  the  Center  of  the  impeller.  The  fluid  test  system  impeller  and 
simulation  impeller  are  shown  in  Figures  5.6  and  5.7,  respectively.  The  simulated  pump  impeller 
is  shown  with  the  inlet-side  plate  removed. 


Figure  5.6:  Fluid  Test  System  Pump  Impeller 


As  stated  before,  moments  of  inertia  were  calculated  from  the  rendering  program.  A  solid 
disk  encompassing  the  entire  pump  rotor  was  first  created  and  its  second  moment  of  volume 
calculated.  The  second  moment  of  volume  for  the  pump  impeller  was  separately  calculated  and 
subtracted  from  the  total  disk  moment  to  arrive  at  a  second  moment  of  volume  for  the  entrained 
water.  These  second  moments  of  volume  are  purely  geometric  quantities,  so  had  to  be  multiplied 
by  material  densities  to  arrive  at  moments  of  Inertia.  Densities  of  1  and  7.85  -2-*  were  used  to 
arrive  at  moments  of  inertia  of  0.000366  kg  —  m2  and  0.00364  kg  —  mr  for  the  water  and  impeller, 
respectively.  The  value  of  density  for  the  pump  impeller  assumes  that  it  is  made  of  stainless  steel. 

An  alternate  method  of  determining  the  moment  of  inertia  for  the  pump  impeller  and  entrained 
water  from  [17]  is  showin  in  equation  5.3.  P  in  the  equation  is  the  shaft  power  in  kW  and  N  is 
the  rotational  speed  in  rpm. 

(p  \  0.9556 

w)  (5'3) 

1Rhino-CAD  is  a  product  of  Robert  McNeel  and  associates 
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Figure  5.7:  Simulated  Pump  Impeller 


According  to  [17],  the  formula  was  developed  from  a  linear  regression  of  “...248  data  points 
from  5  different  pump  manufacturers” .  Using  a  shaft  power  or  1.679  kW  (2.25  hp  as  stated  earlier) 
and  a  rotational  speed  of  1775  rpm  yields  a  moment  of  inertia  of  0.0119  kg  —  m2.  This  value  is 
over  2.5  times  larger  than  the  value  calculated  before.  This  formula  was  not  used  since  there 
is  not  enough  known  about  the  type  of  pumps  considered  in  the  development  of  the  empirical 
relationship.  Specifically,  there  is  not  description  of  the  type  of  pump  impellers  used  and  it  was 
not  known  if  pumps  of  low  power  (less  than  5  hp)  were  included. 


5.3.2  Hydrodynamic  Considerations 

To  determine  the  amount  of  torque  on  the  pump  impeller  as  a  function  of  speed,  the  power 
coefficient  of  equation  5.1  was  again  employed.  In  this  case,  it  is  not  maintained  constant  for 
two  pump  impellers  which  are  geometrically  similar  but  for  a  pump  impeller  rotating  at  different 
speeds.  Since  the  density  of  the  fluid  and  pump  impeller  diameter  do  not  change  with  speed, 
equation  5.1  reduces  to  equation  5.4. 


U!3  U>$ 


(5.4) 


Since  mechanical  power  is  just  the  torque  placed  on  the  rotor  times  its  angular  speed,  this 
further  reduces  to  equation  5.5. 


(jJ 

Ta;  =  P0-3  (5.5) 

-  •  ...  u0 

Dividing  through  by  the  angular  speed  gives  an  expression  for  torque  in  terms  of  angular 
speed  as  in  equation  5.6. 


T  =  P0 


Wn 


(5.6) 
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Finally,  since  it  was  decided  that  the  pump  would  draw  2.25  hp  at  1775  rpm  (185.88  r^p),  the 
final  expression  for  torque  is  developed  in  equation  5.7.  It  should  be  noted  that  this  development 
of  torque  vs.  angular  speed  is  inherently  quasi-steady.  Applying  this  relationship  to  determine 
the  hydrodynamic  torque  assumes  that  the  stator  is  stepping  through  a  series  of  steady  states  as 
it  speeds  up  without  any  consideration  of  more  complicated  hydrodynamics  such  as  transitions 
from  laminar  to  turbulent  flow.  This  is  adequate  however,  because  the  simulation  was  used  to 

explore  changes  in  steady  state  power  levels,  not  changes  in  the  start  up  transient  for  the  pump. 

!  '  • 

T  =  0.0002614a;2  (5.7) 

The  scripts  used  to  implement  the  centrifugal  pump  simulation  are  contained  in  Appendix  E. 

5.4  Steady  State  Power  Level  as  an  Indicator  of  Pump 
Performance 

There  are  fundamentally  three  “bands”  of  pump  power  discussed  in  this  section.  The  first  is  the 
band  that  corresponds  to  operation  with  the  overboard  valve  throttled  to  various  levels.  This 
normal  operation  band  can  be  used  to  predict  the  total  amount  of  flow  from  the  pump.  The 
second  band  corresponds  to  a  problem  on  the  inlet  side  of  the  pump,  such  as  a  clogged  inlet 
strainer,  resulting  in  a  lower  steady  state  power.  The  last  band  is  below  the  second  power  level 
and  corresponds  to  a  problem  with  the  coupling  between  the  pump  motor  and  the  pump  impeller. 
For  this  last  case,  power  levels  during  an  uncoupling  even  drop  below  zero. 

5.4.1  The  NILM  as  an  Indicator  of  Flow 

Test  Data 

One  of  the  first  goals  of  placing  the  NILM  on  the  ASW  system  of  the  SENECA  was  to  see  if 
steady  state  values  of  power  could  be  used  as  an  indicator  of  flow  through  the  system.  To  that 
end,  the  system  was  run  with  the  overboard  valve  throttled  to  various  degrees  and  power  data 
collected.  An  HVAC  unit  and  two  diesel  generators  were  aligned  to  the  system  as  heat  loads.  As 
previously  discussed,  the  overboard  valve  is  throttled  to  control  the  amount  of  flow  to  the  heat 
loads  in  the  system.  If  the  overboard  valve  is  throttled  down,  there  is  less  water  flowing  directly 
over  the  side  and  therefore  more  water  flowing  to  the  heat  loads  (see  Figure  5.1). 

The  overboard  throttling  valve  is  a  butterfly  valve  that  has  discrete  settings  from  notches  cut 
in  the  operating  handle.  A  valve  of  this  type  (but  not  the  actual  valve  in  question)  is  shown  in 
Figure  5.8. 

For  each  setting  of  the  overboard  throttling  valve,  the  pump  discharge  pressure,  heat  load 
pressure  (pressure  downstream  of  the  heat  load  throttling  valve),  and  pump  power  were  collected. 
The  values  of  the  pressure  gauge  zeroes  were  3,0  psig,  1.5  psig,  and  6.0  psig  for  pump  one  discharge 
pressure,  pump  two  discharge  pressure,  and  heat  load  pressure,  respectively.  These  zeroes  were 
subtracted  from  the  readings  taken  from  the  gauges  to  remove  hydrostatic  pressures  from  the 
readings  (as  well  as  account  for  errors  in  the  gauge  zeroes).  Table  5.1  summarizes  data  recorded 
for  ASW  pump  one. 
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Figure  5.8:  ASW  System  Throttling  Valve 


Table  5.1:  Pressure  and  Power  Data  for  Throttling  of  Over¬ 
board  Valve 


Notches 

Throttled 

Pump  Outlet 
Pressure  (psig) 

Heat  Load 
Pressure  (psig) 

Pump 

Power  (kW) 

1 

32 

7.75 

28.43 

2 

32.5 

8 

28.21 

3 

32.5 

8.2 

28.18 

4 

33 

9.5 

28.01 

5 

34 

12.5 

27.87 

6 

35.5 

17.5 

27.17 

7 

41.5 

30.5 

26.68 

Use  of  Pump  Curves  to  Determine  Flow 

Pump  data  was  available  for  the  ASW  pumps  in  [16],  so  the  pressure  information  in  Table  5.1  was 
converted  into  flow  rates.  Among  the  “pump  curves”  plotted  in  [16]  is  the  total  discharge  pressure 
of  .the  pumps  plotted  vs.  flow  rate.  Points  on  this  curve  were  fit  to  a  third-order  polynomial  in 
MATLAB  to  arrive  at  equation  5.8.  In  the  equation,  TDH  is  total  discharge  head  (the  terms 
“head”  and  pressure  will  henceforth  be  used  interchangeably)  in  feet  of  water  while  Q  is  volumetric 
flow  rate  in  gpm. 

TDH  =  -3.3508  •  10~8Q3  +  9.9359  •  10~6Q2  -  4.7494  •  10_3Q  +  110.16  (5.8) 
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For  each  value  of  discharge  head  collected  onboard  the  SENECA,  the  value  of  flow  rate  was 
calculated  by  solving  for  the  flow  rate  that  would  lead  to  that  TDH  in  equation  5.8.  This  flow 
rate  however,  represents  the  total  flow  rate  coming  out  of  the  pump,  not  the  amount  of  flow  going 
to  the  heat  loads. 

Unfortunately,  there  are  too  many  variables  to  calculate  the  amount  of  flow  going  to  the  heat 
load  branch  from  only  two  pressure  readings.  In  order  to  calculate  the  amount  of  flow  going  to 
that  branch,  it  would  be  necessary  to  know  the  loss  -coefficients  for  each  of  the  three  sections  of 
pipe  involved  (pump  outlet  to  split,  split  to  overboard,  and  split  to  heat  loads). 

The  flow  rates  calculated  using  equation  5.8  from  the  collected  data  is  shown  in  Table  5.2  and 
Figure  5.9. 

Table  5.2:  Power  and  Flow  Data  for  Throttling  of  Overboard 
Valve 


Notches  Pump  Total  Flow 
Throttled  Power  (kW)  Rate  (gpm) 


26.6  26.8  27  27.2  27.4  27.6  27.8  28  28.2  28  A  28.6 

Power  (kW) 


Figure  5.9:  Flow  vs.  Power  for  ASW  Pumps 
All  of  the  previous  data  and  discussion  were  from  tests  run  on  pump  one.  The  testing  of  pump 
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one  followed  the  expected  trend  of  pressure  (at  the  pump  discharge  and  in  the  heat  load  line) 
increasing  as  the  overboard  valve  was  throttled  shut.  The  data  also  showed  an  expected  decrease 
in  required  power  for  a  decrease  in  flow  rate.  The  data  collected  for  pump  two  however,  did  not 
follow  these  trends.  An  error  in  testing  may  be  the  culprit,  although  crewmembers  onboard  the 
SENECA  also  expressed  very  little  confidence  in  the  accuracy  of  pressure  gauges  in  the  system. 
Performing  these  tests  again  should  be  a  high  priority  in  future  ASW  system  tests  to  validate  the 
data  collected  for  pump  one. 

Verification  of  Results  using  the  Fluid  Test  System 

A  similar  set  of  tests  was  conducted  for  the  fluid  test  system.  The  system  was  operated  with  water 
flowing  through  the  “aircooler”  and  straight  branches  while  power  data  was  collected.  As  on  the 
SENECA,  the  straight  branch  that  returns  water  directly  to  the  water  reservoir  was  throttled 
to  various  levels.  Unlike  the  ASW  system  onboard  the  SENECA,  the  valves  involved  are  globe 
valves,  so  a  restriction  of  fluid  flow  only  occurs  when  the  valves  are  nearly  shut  (there  are  9  turns 
from  fully  open  to  fully  shut).  Unfortunately,  there  are  no  pressure  gauges  in  the  system,  so  a 
determination  of  flow  could  not  be  made.  Table  5.3  shows  the  results  of  these  tests. 


Table  5.3:  Power  Data  for  Throttling  of  Straight  Return 
Valve  in  Fluid  Test  System 


Turns 

Throttled 

Pump 
Power  (W) 

0 

785.5 

7 

778.9 

7.5 

773.5 

8 

769.3 

8.55 

755.5 

9 

698.5 

The  table  indicates  a  trend  consistent  with  data  collected  from  the  SENECA.  There  is  clearly 
a  decrease  in  power  required  to  run  the  pump  as  more  flow  is  sent  through  the  “aircooler”  branch. 

5.4.2  Prediction  of  Inlet  Strainer  Clogging 

Simulated  Clogging  in  ASW  System 

In  addition  to  a  change  in  head  loss  at  the  outlet  of  the  ASW  pumps  being  reflected  in  pump 
power,  it  was  expected  that  a  change  in  head  loss  at  the  inlet  of  the  ASW  pumps  would  also  be 
seen  in  pump  power.  This  change  could  be  used  as  an  indicator  for  a  blocked  suction  screen  or 
clogged  inlet  strainer. 

To  simulate  clogs  in  the  inlet  strainers  of  the  ASW  pumps  onboard  the  SENECA,  the  strainer 
isolation  valves  were  throttled  to  the  3/4  shut  position  with  the  pumps  off.  The  pumps  were 
each  started  in  this  condition  and  the  power  data  collected.  In  each  case,  there  was  a  marked 
decrease  in  the  power  level  drawn  by  each  of  the  pumps.  An  example  of  a  pump  start  with  the 
inlet  throttled  (and  a  normal  start)  is  shown  in  Figure  5.10. 
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Figure  5.10:  Throttled  Inlet  Starts  for  ASW  System 


It  is  noted  from  Figure  5.10  that  the  transients  themselves  are  very  similar,  it  is  only  the 
steady  state  power  level  which  is  significantly  different.  For  the  pump  start  on  pump  one  with  the 
strainer  isolation  valve  throttled  to  3/4  closed,  the  average  of  1000  points  of  steady  state  power 
data  is  19.96  kW.  Referring  back  to  Table  5.2  shows  that  this  is  nearly  7  kW  lower  than  the  lowest 
normal  steady  state  power  level  for  that  pump.  This  shows  that  it  is  possible  to  distinguish  a 
fault  condition  like  a  clogged  inlet  strainer  from  the  normal  range  of  steady  state  power  levels. 

Simulated  Clogging  in  Fluid  Test  System 

As  with  the  ASW  system  onboard  the  SENECA,  the  fluid  test  system  in  the  laboratory  was  also 
started  with  the  pump  inlet  valve  throttled.  The  pump  was  started  with  the  pump  inlet  valve 
opened  only  0.5  turns  and  the  power  data  collected.  An  example  of  a  pump  start  with  the  inlet 
throttled  is  shown  in  Figure  5,11 

As  with  the  pump  starts  from  the  ASW  system,  the  transients  are  similar  while  the  steady 
state  power  level  is  different.  For  the  pump  start  with  the  pump  inlet  valve  throttled  to  0.5  turns 
open,  the  average  of  2500  points  of  steady  state  power  data  is  611.3  W.  Referring  to  table  5.3 
shows  that  this  is  over  80  W  lower  than  the  lowest  normal  steady  state  power  level  for  the  pump. 
This  corroborates  the  idea  that  a  clogged  inlet  strainer  or  other  obstruction  in  the  inlet  line  is 
distinguishable  from  the  normal  range  of  steady  state  power  levels. 

Head  Loss  at  the  Inlet  of  a  Pump  vs.  the  Outlet 

The  natural  question  now  is:  “Why  should  a  change  to  the  system  at  the  inlet  side  of  the  pump 
have  a  more  pronounced  effect  than  a  change  at  the  outlet  side?”  After  all,  both  are  changes  to  the 
amount  of  flow  resistance  in  the  system.  Although  as  stated  in  [18]  the  total  amount  of  head  loss 
in  a  fluid  system  is  the  sum  of  the  head  loss  upstream  of  the  pump  and  the  head  loss  downstream 
of  the  pump,  the  pump  must  have  a  minimum  pressure  at  its  inlet  to  function  properly. 
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Figure  5.11:  Throttled  Inlet  Starts  for  Fluid  Test  System 


On  the  suction  side  of  a  pump,  low  pressures  are  created.  If  pressure  falls  below  the  vapor 
pressure  of  the  fluid,  bubbles  of  water  vapor  are  formed  in  a  process  called  cavitation.  For 
centrifugal  pumps,  the  minimum  pressure  required  at  the  pump  inlet  to  avoid  cavitation  is  termed 
the  net  positive  suction  head  required  (NPSHR).  If  a  large  pressure  drop  occurs  somewhere 
upstream  of  the  pump  (like  that  caused  by  a  clogged  inlet  strainer),  the  available  net  positive 
suction  head  (NPSHA)  may  fall  below  the  NPSHR  and  cavitation  occurs.  If  cavitation  occurs 
in  the  pump  inlet,  the  hydrodynamic  resistance  may  fall  sharply  since  the  average  density  of  the 
fluid  is  greatly  reduced  by  the  presence  of  vapor  bubbles.  This  will  lead  to  a  sharp  decrease  in 
the  amount  of  power  required  to  run  the  pump. 

5.4.3  Prediction  of  Pump  Uncoupling 

ASW  Pump  Uncoupling 

The  ASW  pumps  onboard  the  SENECA  do  not  have  a  rigid  connection  between  the  pump  motor 
and  pump  impeller.  Instead,  a  flexible  coupling  is  provided  between  the  two.  Couplings  of  this 
kind  are  usually  used  to  reduce  the  vibration  transmitted  down  a  shaft  or  to  prevent  an  over¬ 
torque  condition.  The  coupling  is  a  thick  rubber  ring  with  teeth  on  the  inside  and  outside  edges. 
These  teeth  sit  in  grooves  in  cast-iron  flanges  that  are  connected  to  the  pump  motor  and  pump 
impeller  shafts.  An  example  of  the  coupling  used  to  connect  the  pump  motor  and  pump  impeller 
is  shown  in  Figure  5.12. 

Since  the  coupling  is  actually  composed  of  two  halves,  a  metal  ring  is  used  to  hold  the  two 
halves  together.  This  ring  can  slip  off  during  operation,  causing  the  pump  motor  and  impeller  to 
become  uncoupled.  According  to  a  ship’s  crewmember  “...when  it  breaks  it  is  not  a  clean  break  so 
when  two  high  spots  meet  it  grabs  until  they  slip  by.”  [19]  During  the  performance  of  tests  onboard 
the  SENECA,  a  pump  start  in  which  the  motor  and  impeller  came  uncoupled  was  captured  and  is 
presented  in  Figure  5.13.  The  second  trace  appears  to  be  the  coupling  still  making  some  contact 
with  the  flange  after  the  initial  uncoupling,  exactly  as  described  by  the  crewmember. 
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Figure  5.12:  ASW  Pump  Motor-to-Impeller  Coupling 


Figure  5,13  also  includes  a  start  immediately  after  the  uncoupling  occurred  in  which  the  motor 
is  completely  uncoupled  from  the  pump  impeller.  According  to  the  same  crewmember,  this  is  also 
a  possibility  if  the  set  screw  holding  the  metal  flange  to  the  motor  shaft  comes  loose:  “I  have  not 
had  any  problems  with  the  set  screws  but  they  have  in  the  past  and  the  shaft  will  just  spin  in  the 
coupling.”  [19] 


Figure  5.13:  ASW  Pump  Starts  for  Various  Levels  of  Coupling 
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Simulated  Pump  Uncoupling 

This  uncoupling  event  was  reproduced  by  the  pump  start  simulation  developed  earlier.  To  repro¬ 
duce  the  pump  uncoupling,  the  hydrodynamic  counter-torque  and  moment  of  inertia  were  both 
removed  at  0.8  seconds  into  the  simulation.  At  0.75  seconds  into  the  simulation,  they  were  par¬ 
tially  re-inserted.  This  means  that  a  random  number  between  1  and  0.9  was  used  as  a  multiplier 
of  the  moment  of  inertia  and  the  hydrodynamic  counter-torque.  A  set  of  simulated  starts  similar 
ftp  Figure  5.13  is  shown  in  Figure  5.14. 


Figure  5.14:  Simulated  Pump  Starts  for  Various  Levels  of  Coupling 


The  simulated  pump  starts  produced  in  Figure  5.14  exhibit  the  same  behavior  as  the  ASW 
pumps  shown  earlier  in  Figure  5.13.  It  is  to  note  the  large  decrease  in  pump  power  when  an 
uncoupling  occurs.  This  decrease  is  even  larger  than  would  occur  for  a  clogged  inlet  strainer  at 
the  inlet  side  of  the  pump,  indicating  that  an  automatic  data  collection  system  should  be  able  to 
distinguish  between  a  clogged  inlet  strainer  and  the  pump  impeller  becoming  uncoupled. 

5.5  Power  Spectrums  as  an  Indicator  of  Fluid  System 
Blockage 

Although  the  previous  section  showed  that  steady-state  power  level  can  be  used  as  a  way  to 
determine  the  status  of  the  fluid  system,  they  were  all  based  on  large  changes  to  the  head  loss  in 
the  system.  Small  changes  to  an  individual  component  in  the  fluid  system  however,  can  represent 
much  more  subtle  differences. 

After  the  success  of  the  NILM  to  detect  changes  in  the  fluid  system  at  the  inlet  or  outlet 
of  the  pumps,  it  was  decided  to  see  if  the  NILM  could  detect  changes  at  the  heat  loads  of  the 
system.  To  this  end,  one  of  the  air  conditioning  (HVAC)  units  along  with  two  ship  service  diesel 
generators  (SSDGs)  were  aligned  to  the  system.  The  amount  of  ASW  flow  to  the  heat  loads  was 
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then  throttled  in  various  combinations  and  the  steady-state  power  level  recorded.  The  results  are 
summarized  in  Table  5.4,  The  data  is  for  two  minutes  of  steady-state  operation. 

Table  5.4:  Steady-State  Power  Data  for  Throttling  to  Heat 
Loads 


Condition 

Mean  Power 
(kW) 

Standard  Deviation 
of  Power  (kW) 

No  Throttling 

24.43 

0.75 

SSDGs  Throttled 

24.37 

0.96 

HVAC  Throttled 

24.71 

0.86 

SSDGs  &  HVAC  Throttled 

24.28 

0.75 

As  Table  5.4  indicates,  there  is  a  difference  in  the  amount  of  steady  state  power  level  drawn  for 
throttled  flow  to  the  various  heat  loads,  but  it  is  not  significant  compared  the  standard  deviation 
in  the  power  data.  The  interpretation  of  this  table  is  that  it  is  not  necessarily  possible  to  see  a 
change  in  the  steady-state  power  level  outside  of  that  power  level’s  noise  for  a  change  in  the  flow 
level  to  a  heat  load. 

5.5.1  Fan  Data 

The  first  indication  that  the  frequency  content  of  a  power  signature  could  contain  vital  information 
was  from  a  fluid  system  that  pumped  air  and  not  one  that  pumped  water  (considering  both  air 
and  water  as  fluids;  one  compressible,  the  other  not).  A  NILM  setup  was  installed  at  a  large  fan 
that  is  part  of  the  SENECA’s  ventilation  system.  Although  the  NILM  was  originally  installed 
to  determine  if  clogged  inlet  filters  could  be  detected,  the  shape  of  the  fan’s  startup  transients 
Indicated  that  there  were  some  oscillations  in  the  power  drawn  by  the  fan.  An  illustration  of  the 
fan  is  shown  in  Figure  5.15. 

One  LA  205-S  current  transducer  is  monitoring  the  current  to  the  fan.  The  resistors  used  for 
the  voltage  measurement,  current  measurement,  and  references  are  130  fi,  66.5  O,  and  56.2  fi, 
respectively.  Upon  initial  system  installation,  the  system  was  calibrated  by  collecting  balanced 
three-phase  power  data  with  a  Fluke  39  Powermeter  and  comparing  it  to  the  power  counts  obtained 
with  the  NILM.  The  resultant  scaling  factor  between  real  power  (1.99  kW)  and  counts  (3099.5) 
was  0.642  S- 

The  fan  was  run  with  no  filters  installed,  four  filters  installed,  and  with  two  of  the  filters  covered 
with  cardboard  (to  simulate  blockage).  The  frequency  spectrum  of  power  data  for  the  steady-state 
operation  and  start  transients  were  then  created  using  a  discrete  fourier  transformation  (see  [20] 
for  a  description  of  the  discrete  fourier  transformation).  These  power  spectrums  are  shown  In 
Figures  5.16  and  5.17,  respectively. 

For  the  steady  state  fan  runs,  the  most  important  difference  is  between  the  four  filter  and 
blocked  filter  case.  They  both  show  large  peaks  at  17.5  Hz,  but  the  case  in  which  filters  are 
blocked  shows  a  peak  approximately  25%  higher  than  the  unblocked  case.  For  the  fan  starts, 
there  is  an  oscillation  band  in  the  power  spectrum  which  decreases  in  frequency  as  more  blockage 
is  placed  at  the  inlet  of  the  fan.  The  band  can  be  seen  between  11  and  17  Hz  for  the  no  filter 
case,  between  7  and  10  Hz  for  the  four  filter  case  and  between  3  and  6  Hz  for  the  case  with  two 
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Figure  5.15:  SENECA  HVAC  Fan 


Figure  5.16:  Steady  State  Power  Spectrum  for  Various  Filter  Conditions 

liters  blocked.  These  data  were  the  first  indication  that  blockage  in  a  fluid  system  could  be  seen 
in  a  power  spectrum. 

5.5.2  SENECA  Testing 

Since  is  was  now  known  that  blockage  could  affect  the  power  spectrum,  data  from  the  SENECA 
was  taken  and  analyzed  to  determine  if  blockage  could  be  seen  in  the  power  spectrums  of  the 
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Frequency  (Hz) 

Figure  5.17:  Transient  Power  Spectrum  for  Various  Filter  Conditions 


ASW  pumps.  Blockage  in  heat  exchangers  is  an  actual  problem  that  is  known  to  occur  in  the 
ASW  system  so  was  deemed  a  high  priority.  Many  crewmembers  reported  that  mussels  would 
grow  in  the  endbells  of  heat  exchangers  while  the  ship  was  in  cold  water,  then  detach  themselves 
in  warmer  waters,  causing  a  reduction  of  flow  through  the  heat  exchangers. 

The  first  step  was  to  determine  the  difference  between  the  ASW  system’s  power  spectrum 
with  and  without  mussels  present.  Underway  data  was  found  for  number  one  SSDG  with  mussels 
present  in  the  endbell  of  the  lube  oil  cooler.  A  power  spectrum  was  developed  for  this  data  and 
compared  to  a  power  spectrum  for  the  time  period  two  hours  later  (after  the  mussels  had  been 
removed).  These  two  power  spectrums  are  shown  in  Figure  5.18. 

There  are  two  key  observations  made  from  Figure  5.18: 

1.  There  is  10  Hz  peak  in  power  that  appears  to  be  characteristic  of  the  ASW  system.  This 
peak  may  be  due  to  resonant  effects  anywhere  in  the  system,  including  close  to  the  pump 
itself. 

2,  This  peak  changes  dramatically  for  blockage  in  the  number  one  SSDG  lube  oil  cooler.  The 
difference  in  the  two  peaks  is  nearly  a  factor  of  three.  The  ship’s  logs  only  state  that 
temperature  was  rising  in  number  one  SSDG  and  that  it  was  subsequently  shut  down  at 
0828.  There  is  no  indication  of  rising  temperatures  anywhere  else  in  the  ASW  system,  so  it 
assumed  that  the  change  is  due  only  to  blockage  in  the  number  one  SSDG  lube  oil  cooler. 

In  order  to  verify  the  validity  of  the  power  spectrum  with  blockage  present,  tests  were  run 
with  flow  into  number  one  SSDG  lube  oil  cooler  throttled.  The  throttled  and  unthrottled  results 
are  shown  in  Figure  5.19. 

It  should  be  noted  that  the  peak  at  exactly  10  Hz  shows  the  same  reduction  in  amplitude 
as  in  the  blocked  case.  However,  there  are  peaks  to  the  sides  of  10  Hz  that  show  a  much  higher 
amplitude  for  the  throttled  case  than  for  the  unthrottled  case.  This  phenomenon  was  not  seen 
when  the  lube  oil  cooler  was  actually  clogged  with  mussels.  Clearly,  throttling  an  inlet  valve  to 
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Figure  5.18:  Power  Spectrum  for  Blocked  SSDG  Lube  Oil  Cooler 


Figure  5.19:  Power  Spectrum  for  Throttled  SSDG  Lube  Oil  Cooler 


a  heat  exchanger  is  not  the  same  as  blocking  some  of  the  individual  heat  exchanger  tubes.  The 
difference  is  that  in  the  former  case,  there  is  lessened  flow  to  all  tubes  while  in  the  latter  case, 
the  blocked  tubes  have  reduction  in  flow  while  the  unblocked  tubes  have  an  increase  in  flow  (to 
maintain  the  same  volumetric  flow  rate). 

Ideally,  testing  would  occur  with  obstructions  inserted  into  the  lube  oil  heat  exchanger.  Since 
•testing  of  this  nature  is  not  feasible  on  an  operational  Coast  Guard  Cutter,  the  fluid  test  system 
was  again  utilized  for  more  in-depth  experimentation. 
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5.5.3  Power  Spectrum  and  Vibration  Experiments  using  the  Fluid 
Test  System 

Unlike  tests  conducted  onboard  the  SENECA,  tests  conducted  on  the  fluid  test  system,  also 
involved  the  collection  of  vibration  data.  It  was  felt  that  the  oscillations  in  the  pump  power  might 
be  caused  by  pressure  waves  being  transmitted  back  to  the  pump  impeller.  If  these  pressure  waves 
transmitted  even  some  of  their  energy  into  the  piping,  then  it  should  be  detected  externally. 

To  this  end,  vibration  transducers  were  connected  at  the  inlet  side  of  the  fluid  test  system 
HX.  Each  vibration  transducer  consists  of  a  magnet  which  is  able  to  vibrate  inside  a  loop  of  wire, 
causing  a  voltage  to  be  ouput.  The  vibration  transducers  were  able  to  detect  vibration  along  two 
perpendicular  axes  by  being  attached  to  the  inlet  piping  of  the  HX  as  shown  in  Figure  5.20.  The 
convention  used  was  that  of  cylindrical  coordinates  in  which  the  “Z”  direction  is  along  the  axis 
of  the  cylinder  while  the  “R”  direction  is  outward  from  the  axis  of  the  cylinder. 

* 

|  "R"  Motion 


Figure  5.20:  Heat  Exchanger  Vibrometer  Degrees  of  Freedom 


Throttling  Data 

The  first  experiment  conducted  on  the  fluid  test  system  was  with  the  inlet  valve  to  the  HX 
throttled  to  only  one  turn  open.  This  data  was  taken  to  validate  the  throttling  data  taken 
onboard  the  SENECA.  The  spectrums  for  power,  R  direction  vibration  and  Z  direction  vibration 
are  shown  in  Figures  5.21,  5.22,  and  5.23,  respectively. 

Just  as  on  the  SENECA,  the  power  spectrum  for  throttled  flow  into  the  fluid  test  system  HX 
shows  higher  amplitudes  for  the  throttled  condition.  The  vibration  shows  lower  amplitudes  in 
both  the  R  and  Z  directions,  due  to  the  decrease  in  total  flow  through  the  heat  exchanger. 

Clogging  Data 

Unlike  throttling  the  total  flow  through  the  HX  with  the  use  of  a  single  valve,  there  are  many 
different  ways  that  blocking  of  tubes  can  occur.  Along  with  the  unblocked  condition,  data  were 
taken  for  three  different  types  of  the  blocked-tube  condition.  These  different  modes  of  blocking 
are  illustrated  in  Figure  5.24.  . 

For  each  type  of  clog,  one  minute  of  steady  state  power  and  vibration  data  were  taken  along 
with  a  measurement  of  the  pump’s  impeller  speed.  The  pump  impeller  speed  was  measured  by 
painting  a  small  white  dot  on  the  rotor  of  the  pump  and  then  using  a  laser  tachometer  to  measure 
rpm,  This  process  was  repeated  five  times  for  each  of  the  clogged  conditions,  the  Frequency 
spectrum  for  that  run  of  data  was  created  (for  both  power  and  vibration),  and  the  spectrums  for 
each  of  the  five  rims  were  averaged.  The  average  impeller  speed  was  3520  rpm  with  the  variation 
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Figure  5.21:  Power  Spectrum  for  Throttled  Flow  to  Fluid  Test  System  Heat  Exchanger 


Figure  5.22:  R-Motion  Spectrum  for  Throttled  Flow  to  Fluid  Test  System  Heat  Exchanger 


fri  impeller  speed  within  the  error  of  the  tachometer  (±  1  rpm).  The  frequency  spectrum  for 
pump  power  is  shown  in  Figure  5.25. 

Although  not  necessarily  obvious  from  Figure  5.25,  the  important  frequency  for  the  power 
spectrum  is  at  2.5  Hz.  It  is  important  because  every  test  case  has  a  large  spike  at  that  frequency. 
The  frequency  spectrum  for  vibration  in  the  R  direction  is  shown  in  Figure  5.26. 

Figure  5.26  shows  two  large  peaks  in  frequency,  one  at  58.73  Hz  and  the  other  in  the  vicinity 
of  8.5  Hz.  Since  the  pump  rotational  speed  of  3520  rpm  corresponds  to  a  frequency  of  58.67 
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Frequency  (Hz) 

Figure  5.23:  Z-Motion  Spectrum  for  Throttled  Flow  to  Fluid  Test  System  Heat  Exchanger 


3  Tubes,  6  Tubes 

Asymmetric 


Figure  5.24:  Fluid  Test  System  HX  Clog  Types 


Hz,  the  peak  at  58.73  Hz  is  just  mechanical  vibration  transmitted  through  the  piping  due  to  the 
rotation  of  the  pump  impeller  and  not  due  to  fluid  moving  through  the  pipe.  The  peaks  at  58.73 
Hz  and  8.5  Hz  are  shown  in  Figures  5.27  and  5.28,  respectively. 

There  are  two  interesting  phenomena  to  note  in  Figure  5.28.  First,  for  the  cases  where  there 


Figure  5.25:  Power  Spectrum  for  Clogged  Flow  to  Fluid  Test  System  Heat  Exchanger 
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Figure  5.26:  R-Motion  Spectrum  for  Clogged  Flow  to  Fluid  Test  System  Heat  Exchanger 


Frequency  (Hz) 


-  unclogged 

—  3  tubes  clogged  (symm) 

. .  3  tubes  clogged  (asymm) 

- 6  tubes  clogged  (symm) 


is  symmetric  flow  through  the  HX,  frequency  peaks  do  not  get  larger,  but  they  do  shift  upwards 
in  frequency  as  the  amount  of  clogging  increases.  Secondly,  the  only  frequency  peak  which  is 
larger  than  the  others  is  for  the  case  where  three  tubes  are  clogged  asymmetrically.  Compare 
these  observations  with  Figure  5.25.  Note  that  in  Figure  5.25,  the  frequency  peak  in  the  vicinity 
of  2.5  Hz  for  the  unclogged  case  and  the  case  with  three  tubes  clogged  asymmetrically  are  at  the 
same  frequency.  As  with  vibration  in  the  R  direction,  the  frequency  peak  is  higher  when  there  is 
asymmetric  blockage.  Looking  farther  back  to  Figure  5.16,  note  that  when  two  filters  are  blocked, 
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Figure  5.27:  58.73  Hz  Peaks  in  R- Direction  for  Clogged  Flow  to  Fluid  Test  System  Heat 
Exchanger 
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Figure  5.28:  8.5  Hz  Peaks  in  R-Direetion  for  Clogged  Flow  to  Fluid  Test  System  Heat 
Exchanger 


there  is  a  higher  frequency  peak  (at  the  same  frequency)  than  when  the  fan  filters  are  unblocked. 
The  blocked  filters  in  the  SENECA  HVAC  fan  case  also  represent  asymmetric  blockage  because 
the  top  two  filters  were  blocked  while  the  bottom  two  filters  were  left  unblocked. 

For  vibration  of  the  fluid  test  system  in  the  Z  direction,  Figures  5.29  and  5.30  show  the  fre¬ 
quency  spectrum  for  various  blockage  conditions.  The  largest  frequency  peak  is  at  approximatelev 


74  Hz  with  two  side  bands.  The  peaks  to  either  side  of  74  Hz  are  identified  as  sidebands  because 
they  are  spaced  evenly  around  74  Hz. 


0.035 

0.03 

0.025 

0.02 

0.015 

0.01 

0.005 

0 

50  60  70  80  90  100 

Frequency  (Hz) 

Figure  5.29:  Z-Motion  Spectrum  for  Clogged  Flow  to  Fluid  Test  System  Heat  Exchanger 
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Figure  5.30:  73.88  Hz  Peaks  in  Z-Direction  for  Clogged  Flow  to  Fluid  Test  System  Heat 
Exchanger 
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All  three  of  the  peaks  in  Figure  5.29  show  the  same  trend,  so  the  peak  presented  in  Figure 
5.30  at  73.88  Hz  is  discussed.  The  Z  direction  vibration  does  not  follow  the  same  trends  as  the 
R  direction  vibration  or  the  power  oscillations.  The  Z  direction  vibration  shows  the  highest  peak 
for  the  unclogged  condition,  which  is  not  true  for  the  other  two  spectrums.  Vibration  in  the  Z 
direction  should  probably  not  be  used  as  a  predictor  of  the  power  spectrums  because  it  was  noted 
during  testing  that  the  thrust  of  water  against  the  side  of  the  water  reservoir  in  the  fluid  test 
system  was  causing  vibration  in  the  Z  direction  of  the  HX. 

To  summarize  the  trends  in  the  the  various  frequeny  spectrums,  Table  5.5  provides  rankings 
of  the  frequency  peaks  for  the  power,  R  direction  vibration,  and  Z  direction  vibration  spectrums. 
The  highest  peak  is  given  a  rank  of  one  while  the  lowest  peak  is  given  a  rank  of  four.  Shifts  in 
frequency  are  also  indicated  where  applicable. 


Table  5.5:  Frequency  Peak  Rankings  for  Various  Fluid  Test 
System  HX  Clogs 


Clog 

Power 

R 

R  (58.73  Hz) 

Z  (73.88  Hz) 

None 

3  (2.7  Hz) 

2  (8.1  Hz) 

2 

1 

3  (symm) 

2  (2.63  Hz) 

2  (8.5  Hz) 

4 

3 

3  (asymm) 

1  (2.7  Hz) 

1  (8.3  Hz) 

3 

2  (tie) 

6 

4  (2,63  Hz) 

2  (8.7  Hz) 

1 

2  (tie) 

It  is  noted  from  Table  5,5  that  there  is  no  clear  trend  in  the  spectrum  of  frequencies  that 
corresponds  to  the  clog  severity.  It  is  noted  however,  that  the  situation  which  is  most  severe 
shows  a  large  decrease  in  the  frequency  peak  in  the  vicinity  of  2.5  Hz.  If  2.5  Hz  is  the  frequency 
that  is  an  indicator  of  HX  clogging,  a  decrease  in  the  peak  at  that  frequency  for  a  highly  clogged 
HX  is  consistent  with  the  results  from  the  SENECA.  It  is  also  noted  from  Table  5.5  that  the 
upward  shift  in  frequency  for  vibration  in  the  R  direction  (in  the  vicinity  of  8.5  Hz)  may  be  the 
best  indicator  of  clogging. 

5,6  Chapter  Summary 

Large  changes  in  head  loss  in  a  fluid  system  are  readily  apparent  in  the  steady  state  power  levels 
of  the  pumps  which  provide  flow  through  the  system.  It  is  possible  to  distinguish  between  large 
head  losses  at  the  pump  inlet  and  the  pump  outlet. 

Subtle  changes  in  the  fluid  system  at  individual  components  are  much  harder  to  detect.  It 
is  possible  that  each  component  in  the  fluid  system  has  a  unique  signature  in  the  spectrum  of 
frequencies  over  which  power  is  transmitted.  It  remains  to  be  seen  how  changes  to  individual  fluid 
system  components  affect  these  frequency  signatures.  Specifically,  a  more  detailed  understanding 
of  the  fluid  mechanics  of  clogged  heat  exchangers  must  be  developed. 


Chapter  6 

Hydraulic  System  Results 


6.1  Introduction 

Like  seawater  cooling  systems,  hydraulic  systems  are  common  across  a  wide  range  of  ship  types 
and  sizes.  One  of  the  most  important  type  of  hydraulic  system  found  onboard  ships  is  the  steering 
system.  Since  the  steering  system  (in  conjunction  with  the  main  engines)  provides  the  ship  its 
maneuverability,  the  status  of  its  functionality  must  be  known  at  all  times.  Monitoring  the 
steering  systems  onboard  the  USCGC  SENECA  and  the  ONR  YP  provided  on  opportunity  to 
determine  if  the  NILM  could  detect  any  degradation  of  the  system’s  performance.  Additionally, 
it  was  hoped  that  the  NILM  could  provide  a  backup  indication  of  system  parameters  (e.g.  rudder 
angle)  with  the  hope  of  developing  a  NILM-based  autopilot  in  the  future. 

6.2  SENECA  Steering  System 

The  steering  system  onboard  the  SENECA  consists  of  two  30  hp  pumps  drawing  hydraulic  fluid 
from  a  common  reservoir  to  pressure  a  common  header.  Rudder  actuation  is  accomplished  through 
two  4  in.  hydraulic  rams  (one  for  each  rudder  stock).  Additionally,  parallel  rudder  motion  is 
accomplished  through  the  use  of  a  large  tie  rod  which  connects  the  two  rudder  stocks.  Figure  6.1 
shows  a  schematic  representation  of  the  SENECA’s  steering  system  while  Figure  6.2  shows  one 
of  the  rudder  hydraulic  rams  and  the  connecting  tie  rod. 

Different  system  configurations  are  achieved  by  repositioning  valves.  The  valve  positions 
requird  for  these  different  configurations  is  shown  in  Figure  6.3.  The  diagram  is  a  copy  of  the 
actual  plate  on  the  side  of  the  steering  system  itself. 

The  NILM  installation  for  the  SENECA  system  uses  one  voltage  transducer  and  two  current 
transducers  to  monitor  the  two  hydraulic  pumps  that  pressurize  the  system.  Since  there  is  no  one 
panel  that  supplies  both  of  the  pumps  (for  reliability),  one  LA-205  current  transducer  had  to  be 
ffiaced  in  each  pump  controller  and  the  output  combined  at  the  NILM  box.  The  resistors  used 
for  the  voltage  measurement,  current  measurement,  and  references  are  130  fi,  15  fi,  and  56.2  Q, 
respectively.  Upon  initial  system  installation,  the  system  was  calibrated  by  collecting  balanced 
three-phase  power  data  for  pump  one  with  a  Fluke  39  Powermeter  and  comparing  it  to  the  power 
counts  obtained  with  the  NILM.  The  resultant  scaling  factor  between  real  power  (3.6  kW)  and 
counts  (570.11)  was  6.31  . 

It  should  be  noted  that  the  steering  system  onboard  the  SENECA  is  the  only  NILM  installation 
discussed  in  this  thesis  that  did  not  have  its  own  dedicated  computer.  Test  data  were  taken  with 
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Figure  6.1:  SENECA  Steering  System  Schematic 


Figure  6.2:  SENECA  Steering  System  Rudder  Ram  and  Tie  Rod 


a  computer  that  was  connected  to  the  NILM  box  only  at  the  time  of  testing.  Two  different  test 
sessions  used  two  different  computers.  It  was  discovered  that  the  number  of  counts  produced  by 
the  second  computer  for  the  same  pump  power  was  different  than  for  the  first  computer.  The 
resultant  scaling  factor  between  real  power  (3.6  kW)  and  counts  (1033.3)  was  3.48  for  the 
second  computer. 

The  first  scaling  factor  was  used  for  data  taken  in  October,  2003  while  the  second  scaling 
factor  was  used  for  data  taken  in  2004, 
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Figure  6.3:  SENECA  Steering  System  Valve  Lineups 


6.2.1  In-Port  Rudder  Transients  Before  and  After  a  Repaired 
Cylinder 

The  rudder  transients  that  were  conducted  consisted  of  “fishtailing”  the  rudder  back  and  forth 
through  various  angles.  The  rudder  was  first  moved  from  0  degrees  to  left  5  degrees.  After  a  brief 
pause,  the  rudder  was  moved  back  to  0  degrees.  After  another  brief  pause,  the  rudder  was  moved 
to  right  5  degrees.  The  rudder  was  then  moved  back  to  0  degrees  following  another  brief  pause. 
This  process  was  continued  for  10  degrees,  15  degrees,  etc.  up  to  25  degrees  of  rudder  deflection. 
The  power  data  for  these  transients  are  shown  in  Figure  6.4.  Notice  that  the  transients  are  longer 
as  the  rudder  is  deflected  to  higher  angles. 

The  data  contained  in  Figure  6.4  was  one  of  the  first  sets  of  data  taken  for  the  steering  system. 
After  collecting  this  data,  it  was  desired  to  disable  one  of  the  rudder  ramss  to  determine  how  the 
transients  would  look  for  rudder  motion  in  which  only  one  of  the  rams  was  providing  power.  Due 
to  an  incorrect  valve  lineup,  the  return  path  for  the  hydraulic  ram  that  had  been  disabled  (number 
one)  was  closed,  so  the  cylinder  began  leaking  hydraulic  fluid.  It  was  later  discovered  that  an 
O-ring  had  become  dislogded  from  the  excessive  pressure  in  the  cylinder.  The  cylinder  and  ram 
were  removed  for  repair  and  no  further  testing  was  accomplished  before  the  ship  went  out  to  sea. 

The  power  data  in  Figure  6.4  show  a  variation  in  the  amount  of  power  drawn  when  the  rudders 
are  moving  left  or  right.  This  is  due  to  a  difference  in  the  two  cylinders.  When  these  data  were 
collected,  it  was  noted  that  cylinder  number  one  would  shake  when  the  ram  was  retracting  into  its 
cylinder,  which  corresponds  to  a  left  turn.  The  transients  that  correspond  to  the  rudders  deflecting 
left  show  a  lower  power  value  than  the  transients  that  correspond  to  the  rudders  deflecting  right. 

Following  the  SENECA’s  return  to  port,  a  new  round  of  testing  was  begun  with  cylinder 
number  one  repaired  and  reinstalled  in  the  system.  The  first  test  was  the  same  “fishtailing”  of 
the  rudders  that  had  been  done  before.  The  power  data  for  these  transients  with  the  repaired 
cylinder  are  shown  in  Figure  6.4 

The  power  data  in  Figure  6.5  do  not  show  a  variation  in  the  amount  of  power  drawn  when 
the  rudders  are  moving  left  or  right.  There  was  clearly  a  fault  present  in  number  one  cylinder 
even  before  the  overpressure  condition  caused  the  O-ring  to  become  dislodged.  Whatever  fault 
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Figure  6.4:  SENECA  Rudder  Transients,  Both  Rams,  Pump  1  Running 


Figure  6.5:  SENECA  Rudder  Transients,  Both  Rams,  Pump  1  Running,  After  Repairs 


was  present  in  the  eyinder  might  have  actually  been  the  cause  of  the  dislodged  O-ring.  This  is 
corroborated  by  the  failure  of  the  system’s  relief  valve  (which  several  people  heard  lift)  to  protect 
the  cylinder.  ,  : 


6.2.2  Rudder  Transients  with  Disabled  Cylinders 

The  same  rudder  “fishtailing”  that  was  conducted  before  was  again  conducted  to  determine  if 
there  was  a  difference  in  power  transients  when  only  one  ram  was  functioning  instead  of  both  of 
them.  The  collected  power  data  are  shown  in  Figures  6.6  and  6.7. 
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Figure  6.6:  SENECA  Rudder  Transients,  Ram  1  only 


Figure  6.7:  SENECA  Rudder  Transients,  Ram  2  only 


Figures  6.6  and  6.7  show  that  when  there  is  only  one  powered  ram,  there  is  a  difference  between 
the  power  drawn  by  the  hydraulic  pumps  for  left  rudder  movement  and  right  rudder  movement. 
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If  ram  one  is  powered,  more  power  is  drawn  for  a  left  turn,  whereas  more  power  is  drawn  for  a 
right  turn  if  ram  two  is  powered.  This  difference  can  be  explained  by  the  difference  in  head  loss 
experienced  at  one  side  of  the  ram  vs.  the  other. 

According  to  [1],  the  amount  of  head  loss  experienced  by  a  fluid  moving  at  some  velocity  is 
given  by  equation  6.1  where  hi,  is  the  head  loss,  Kl  is  the  head  loss  coefficient,  V  is  the  fluid 
velocity,  and  g  is  acceleration  due  to  gravity.  The  head  loss  coefficient  in  this  equation  is  usually 
only  a  function  of  the  piping  geometry  under  consideration. 

V2 

hL  =  Kl—  (6.1) 

For  the  purposes  of  this  analysis,  an  expansion  is  defined  as  when  a  fluid  moves  from  a  pipe  of 
small  cross-sectional  area  to  a  pipe  of  larger  cross-sectional  area.  Likewise,  a  contraction  is  when 
a  fluid  moves  from  a  pipe  of  large  cross-sectional  area  to  a  pipe  of  smaller  cross-sectional  area. 
For  a  cross-sectional  area  ratio  of  ^  where  Ai  is  the  cross-sectional  area  of  the  large  pipe  and  A 2 
is  the  cross-sectional  area  of  the  small  pipe,  the  head  loss  coefficients  for  the  cases  of  contraction 
and  expansion  are  given  in  Figure  6.8. 


Figure  6.8:  Head  Loss  Coefficients  for  Expansion  and  Contraction 

The  steering  system  on  the  SENECA  is  composed  of  4  in.  cylinders  and  1  in.  hydraulic 
lines.  Therefore,  Operation  of  the  steering  system  consists  of  fluid  expansions  and  contractions  to 
and  from  the  hydraulic  cylinders.  When  both  of  the  rams  are  powered,  the  head  losses  due  to 
expansion  and  contraction  are  equal  for  left  or  right  rudder  deflection.  If  only  one  ram  is  powered, 
however,  the  head  losses  are  different  for  left  and  right  rudder  deflections. 

It  is  clear  that  the  side  of  the  hydraulic  cylinder  which  does  not  contain  the  ram  (e.g.  the 
left  side  of  cylinder  one  in  Figure  6.1)  has  a  larger  cross-sectional  area  through  which  fluid  can 
move.  For  the  steering  system  onboard  the  SENECA,  the  ratio  for  the  side  of  the  cylinder 
which  does  not  contain  the  ram  and  the  hydraulic  piping  is  0.25.  The  estimated  ^  ratio  for  the 
side  of  the  cylinder  which  contains  the  ram  and  the  hydraulic  piping  is  0.5.  Using  these  area 
ratios,  it  is  possible  to  generate  a  table  for  the  head  loss  coefficients  in  the  two  hydraulic  cylinders 
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for  left  and  right  rudder  deflection  for  only  one  powered  cylinder.  Table  6.1  gives  the  head  loss 
coefficients  from  Figure  6.8  for  left  and  right  rudder  deflection  with  only  cylinder  one  powered. 
Only  the  portion  of  each  cylinder  that  acts  as  a  load  on  the  system  is  considered. 


Table  6.1:  Hydraulic  Head  Loss  Coefficients  for  Left  and 
Right  Rudder  Motion  w/  Ram  1  Powered 


Left  Deflection 

Right  Deflection 

Cyl.  1  Contraction 

0.425 

0.375 

Cyl.  2  Expansion 

0.5625 

0.25 

Cyl.  2  Contraction 

0.425 

Total 

1.3625 

1.05 

Table  6.1  indicates  that  the  total  head  loss  coefficient  for  left  rudder  deflection  is  higher  than 
for  right  rudder  deflection  when  only  cylinder  one  is  powered.  This  means  that  there  is  a  larger 
pressure  drop  for  left  rudder  deflections,  so  the  hydraulic  pumps  must  consume  more  power  to 
maintain  pressure  in  the  system,  exactly  as  is  shown  in  Figure  6.6. 

For  the  situation  in  which  cylinder  two  is  powered,  the  entries  in  Table  6.1  are  reversed,  and 
right  rudder  deflection  requires  a  higher  amount  of  hydraulic  pump  power.  This  is  the  situation 
shown  in  Figure  6.7. 

6.3  ONR  YP  Steering  System 

In  addition  to  the  extensive  testing  that  was  performed  on  the  USCGC  SENECA,  tests  were  also 
performed  on  a  yard  patrol  craft  operated  by  the  Office  of  Naval  Research  (referred  to  as  the 
“ONR  YP”).  The  steering  system  onboard  the  ONR  YP  consists  of  two  3  hp  pumps  drawing 
hydraulic  fluid  from  a  common  reservoir  to  pressure  a  common  header.  Unlike  the  steering  system 
on  the  SENECA,  rudder  actuation  is  accomplished  through  only  one  hydraulic  ram.  The  single 
hydraulic  ram  moves  in  two  different  directions  to  accomplish  left  and  right  rudder  deflections. 
The  hydraulic  ram  moves  the  starboard  rudder  directly.  Parallel  motion  of  the  port  rudder  is 
achieved  through  the  use  of  a  tie  rod  which  connects  the  port  an  starboard  rudder  stocks.  Figures 
6.9  shows  a  schematic  of  the  system  while  Figures  6.10  and  6.11  shows  one  of  the  hydraulic  pumps 
that  pressurizes  the  system  and  the  system  tie  rod,  respectively. 

The  NILM  installation  for  the  ONR  YP  steering  system  uses  one  voltage  transducer  and  two 
current  transducers  to  monitor  the  two  hydraulic  pumps  that  pressurize  the  system.  Since  there 
is  no  one  panel  that  supplies  both  of  the  pumps  (for  reliability),  one  LA-55  current  transducer 
bad  to  be  placed  in  each  pump  controller  and  the  output  combined  at  the  NILM  box.  The 
resistors  used  for  the  voltage  measurement,  current  measurement,  and  references  are  47  Q,  95  O, 
and  47  Q,  respectively.  Upon  initial  system  installation,  the  system  was  calibrated  by  collecting 
balanced  three-phase  power  data  with  a  Fluke  39  Powermeter  and  comparing  it  to  the  power 
counts  obtained  with  the  NILM.  The  resultant  scaling  factor  between  real  power  (1.19  kW)  and 
counts  (3938.7)  was  0.3021  SgSa. 
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Figure  6.9:  ONR  YP  Steering  System  Schematic 


Figure  6.10:  ONR  YP  Steering  System  Hydraulic  Pumps 

6.3.1  Normal  In-Port  Rudder  Transients 

The  rudder  transients  that  were  conducted  consisted  of  “fishtailing”  the  rudder  back  and  forth 
through  various  angles,  just  as  on  the  SENECA.  The  power  data  for  these  transients  are  shown 
in  Figure  6.12. 

It  is  important  to  note  in  Figure  6.12  that  the  power  data  is  not  nearly  as  clear  as  for  the 
rudder  transients  collected  onboard  the  SENECA.  Not  only  is  there  a  larger  variation  in  the 
steady  state  power  level  between  power  levels,  but  the  rudder  transients  themselves  do  not  show 
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Figure  6.11:  ONR  YP  Steering  System  Tie  Rod 


Figure  6.12:  ONR  YP  Rudder  Transients 


the  flat  plateau  that  the  rudder  transients  for  the  SENECA  did.  This  is  not  necessarily  bad 
news  however,  the  differences  in  power  transients  can  actually  be  used  to  determine  the  nature  of 
rudder  motion.  Figure  6.13  shows  the  last  three  rudder  transients  from  Figure  6.12. 

In  Figure  6.13,  the  rudder  is  moving  from  left  25  degrees  to  0  degrees,  from  0  degrees  to  right 
25  degrees,  and  finally  from  right  25  degrees  back  to  0  degrees.  Figure  6.13  clearly  shows  that  the 
rudder  transient  for  deflection  towards  0  degrees  has  a  much  larger  initial  power  peak  than  rudder 
motion  away  from  0  degrees.  Ultimately  this  information  could  be  used  by  an  autopilot  (along 
w/  recognition  in  the  difference  between  left  and  right  rudder  deflection)  to  determine  rudder 
position. 

6.3.2  Faulted  In-Port  Rudder  Transients 

In  addition  to  collecting  rudder  transients  with  both  rudders  being  moved  by  the  single  hydraulic 
ram,  transients  were  also  collected  with  the  tie  rod  disconnected  so  that  the  hydraulic  ram  was 
moving  only  the  starboard  rudder.  The  same  “fishtailing”  transients  were  collected  as  for  when 
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Figure  6.13:  Larger  View  of  ONR  YP  Rudder  Transients 


the  tie  rod  was  connected.  The  power  data  for  these  transients  are  shown  in  Figure  6.14. 
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Figure  6.14:  ONR  YP  Rudder  Transients,  Starboard  Rudder  Only 


It  should  noted  from  a  comparison  of  Figures  6.12  and  6.14  that  the  power  levels  for  the 
tie  rod  disconnected  are  higher  than  for  the  tie  rod  connected.  It  might  at  first  seem  that  this 
higher  power  level  might  be  a  good  indicator  of  the  tie  rod  being  disconnected.  The  discrepancy 
is  that  the  power  drawn  by  the  pump  when  there  is  no  rudder  transient  present  is  also  higher. 
Since  disconnecting  the  tie  rod  between  rudder  stocks  does  not  in  any  way  alter  the  lineup  of  the 
hydraulic  piping,  there  is  no  reason  it  should  be  drawing  more  power.  It  is  better,  therefore  to 


look  at  the  difference  between  the  power  drawn  during  a  rudder  transient  and  the  power  drawn 
with  no  rudder  transient.  Figure  6.15  shows  the  last  three  rudder  transients  from  Figure  6.14. 


Figure  6.15:  ONR  YP  Rudder  Transients,  Starboard  Rudder  Only 


As  for  the  transients  with  the  tie  rod  connected,  Figure  6.15  shows  the  rudder  moving  from  left 
25  degrees  to  0  degrees,  from  0  degrees  to  right  25  degrees,  and  finally  from  right  25  degrees  back 
to  0  degrees.  The  characteristic  difference  between  transients  for  deflection  towards  0  degrees 
and  transients  for  deflection  away  from  0  degrees  are  still  present.  The  power  level  drawn  for 
transients  is  also  about  200  W  higher  than  the  baseline  power  level,  just  as  for  the  case  with  the 
tie  rod  connected. 

The  only  difference  in  transients  with  the  tie  rod  disconnected  is  the  shape  of  the  power  drawn 
during  the  transients.  In  Figure  6.13  the  power  transients  with  the  tie  rod  connected  show  more 
oscillations  during  the  rudder  deflections.  The  power  transients  in  Figure  6.15  however,  show 
much  smaller  oscillations  in  power. 

6.4  Chapter  Summary 

Based  on  the  results  presented  previously,  it  is  possible  to  detect  faults  associated  with  hydraulic 
cylinders  in  steering  sytems  onboard  ships.  For  the  two-cylinder  system  onboard  the  SENECA,  a 
disabled  cylinder  is  indicated  by  a  difference  in  the  steady  state  power  level  drawn  for  left  and  right 
rudder  deflection.  For  the  single-cylinder  system  onboard  the  ONR  YP,  the  difference  between 
one-rudder  deflection  and  two-rudder  deflection  is  indicated  by  the  shape  of  power  transient. 

For  the  purposes  of  providing  an  indication  of  rudder  angle,  the  two-cylinder  system  onboard 
the  SENECA  does  not  currently  show  an  obvious  difference  in  left  and  right  rudder  deflection. 
The  single-cylinder  system  onboard  the  ONR  YP  only  shows  a  difference  in  deflection  away  from 
zero  versus  deflection  towards  zero.  If  a  method  is  developed  to  distinguish  between  left  and 
right  rudder  deflections,  rudder  angle  can  be  easily  determined  from  the  duration  of  the  rudder 
deflection  transient. 
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Chapter  7 
Future  Work 


7.1  Current  Systems 

7.1.1  NILM  Hardware 

The  NILM  Hardware  used  in  the  installations  discussed  in  this  thesis  was  able  to  collect  data  for 
analysis  satisfactorily,  but  there  is  room  for  improvement. 

Although  the  analog-to-digital  conversion  card  used  in  the  PCs  performed  well,  they  are 
extremely  expensive. .  A  universal  serial  bus  (USB)  interface  developed  at  MIT  by  fellow  LEES 
students  will  be  used  in  the  next  generation  of  NILM  installations.  In  addition  to  being  more 
cost-effective,  the  USB  NILM  promises  a  much  higher  data  collection  rate  than  the  current  SCSI 
system. 

A  more  advanced  interface  between  the  voltage  and  current  transducers  and  the  PC  should  also 
be  investigated.  If  wiring  terminal  board  resistors  are  sized  to  prevent  an  over-voltage  condition 
in  the  analog-to-digital  card  for  a  load  start,  fidelity  is  lost  in  the  steady  state  signal.  An  interface 
needs  to  be  developed  that  adjusts  the  sensitivity  of  the  signal  to  the  computer  based  on  the 
amount  of  current  being  drawn  by  the  load.  This  hardware  would  ideally  prevent  damage  during 
a  large  current  peak  but  provide  a  high  amount  of  fidelity  for  steady  state  signals. 

In  some  instances  (especially  where  power  level  and  not  frequency  content  is  important)  high 
or  low  pass  filters  should  also  be  considered  at  the  the  interface  between  the  transducers  and  the 
PC.  This  might  allow  recognition  of  much  smaller  changes  in  steady  state  power  level. 

7.1.2  NILM  Software 

The  current  state  of  the  art  in  NILM  software  consists  of  the  preprocesser  that  converts  the  voltage 
and  current  signals  into  a  power  envelope,  the  automatic  data  collection  discussed  in  this  thesis, 
and  “demonstrator”  software  not  discussed  in  this  thesis  that  has  been  developed  to  track  loads 
turning  on  and  off.  With  the  demonstrator  software  in  “train”  mode,  load  on  and  off  transients 
are  shown  to  the  system.  When  the  software  is  switched  to  “run”,  it  is  able  to  recognize  on  and 
off  transients  for  any  load  that  it  has  been  trained  to  see. 

The  next  step  in  NILM  software  should  combine  automatic  data  acquisition  and  storage,  load 
on/off  recognition,  and  recognition  of  large  changes  in  steady  state  power  level.  Although  not 
expected  to  have  the  advanced  capabilities  of  the  ideal  system  listed  in  Chapter  1,  it  should  be 
able  to  trend  equipment  performance  over  time,  recognize  faults  and  create  the  corresponding 
error  message,  and  store  data  for  later  analysis. 


I 


83 


7.1.3  Cycling  Systems 

Although  the  “leak  indicator”  developed  earlier  is  a  good  start  for  detecting  a  constant  drain  on 
a  cycling  system,  it  needs  to  be  refined  through  further  testing  and  simulation  . 

Improvements  to  the  Sewage  System  Simulation 

First,  the  simulation  is  based  on  pressure,  and  not  leak  flow  rate.  Leaks  in  a  vacuum  systems  do 
not  actually  develop  as  a  constant  loss  of  vacuum,  but  rather  as  an  intake  of  air  into  the  system 
due  to  a  loss  of  system  integrity  (e.g.  loose  fittings  or  holes  in  piping).  This  air  intake  is  not 
constant  over  time  for  a  loss  of  system  integrity,  but  actually  varies  with  the  amount  of  vacuum 
in  the  system. 

Second,  the  probably  density  function  for  flushes  used  was  probably  not  very  realistic.  It 
is  probably  more  realistic  to  use  a  probability  density  function  whose  cumulative  distribution 
function  approaches  one  as  time  goes  to  infinity.  The  binomial  or  Poisson  distributions  are  good 
candidates  (see  [13]),  they  were  not  originally  used  because  of  the  need  to  generate  large  amounts 
of  data  quickly. 

Cross-System  Validation 

Another  way  that  the  leak  indicator  can  be  improved  upon  is  through  cross-system  validation. 
Experiments  should  be  run  on  other  types  of  cycling  systems  to  see  if  they  are  affected  in  the 
same  way  as  the  sewage  system  onboard  the  SENECA.  The  HVAC  system  onboard  the  ONR  YP 
is  a  good  candidate  for  such  a  system  because  it  already  contains  a  NILM  installation. 

The  NILM  installation  for  the  ONR  YP  is  using  one  voltage  transducer  and  one  LA  205-S 
current  transducer  to  monitor  two  compressors,  two  heater  banks,  one  seawater  pump,  and  one 
chillwater  pump.  A  schematic  of  the  system  is  shown  in  Figure  7.1. 

Chillwater  is  pumped  through  various  thermostatically-controlled  spot  coolers  around  the 
ship,  Chillwater  flow  through  these  spot  coolers  is  constant  while  airflow  through  them  cycles  on 
and  off  to  maintain  temperature  in  each  zone.  The  compressors  for'  the  HVAC  system  in  turn 
cycle  on  and  off  to  maintain  chillwater  temperature.  An  example  of  compressor  cycling  is  shown 
in  Figure  7.2. 

The  ONR  YP  HVAC  system  would  be  an  excellent  system  to  trend  because  it  is  affected  by 
many  variables.  Not  only  are  there  many  different  combinations  of  thermostat  settings  which 
could  be  tried,  but  trends  could  also  be  analyzed  for  changes  in  air  and  seawater  temperature. 

7.1.4  Fluid  Systems 

Although  large  changes  to  the  head  loss  of  a  fluid  system  (like  those  is  section  5.4)  are  easy  to 
see  in  the  steady  state  power  level  of  a  pump,  it  has  already  been  shown  that  small  changes  in  a 
fluid  system  are  not  as  easily  detectable. 

The  equipment  used  thusfar  to  detect  pressure  waves  which  may  be  influencing  the  frequency 
spectrum  of  pump  power  measures  vibration  in  piping,  not  fluid  pressures.  There  are  too  many 
things  which  affect  the  vibration  of  the  piping  other  than  fluid  flow,  primarily  mechanical  vibration 
of  the  pump  itself. 

In  the  future,  pressure  transducers  should  be  installed  in  the  fluid  test  system  to  directly 
measure  any  pressure  waves  that  are  formed  in  the  piping.  This  would  better  allow  the  frequency 
of  interest  in  the  power  spectrum  to  be  pinpointed  to  predict  clogging  failures. 
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Figure  7.2:  Compressor  Cycling  in  ONR  YP  HVAC  System 


A  good  candidate  is  the  K-series  pressure  transducer  manufactured  by  Ashcroft.  Available 
from  Grainger  Tools,  this  pressure  transducer  comes  in  a  wide  variety  of  pressure  ranges,  and  has 
a  stainless  steel  membrane  to  ensure  longevity. 
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7.1.5  Hydraulic  Systems 

The  results  for  the  steering  systems  onboard  the  SENECA  and  ONR  YP  are  all  for  inport  data. 
There  are  two  drawbacks  to  collecting  data  for  a  steering  system  while  moored.  First,  there  are  no 
hydrodynamic  forces  on  the  rudders  as  a  result  of  the  ship’s  motion  through  the  water.  Second, 
the  system  is  powered  by  shore  power  instead  of  the  ship’s  onboard  generating  equipment.  To 
predict  performance  while  underway,  tests  must  be  conducted  underway. 

Some  data  was  collected  with  the'  ONR  YP  underway.  Turns  Were  made  at  different  speeds 
with  different  rudder  angles,  but  there  was  a  constant  need  to  avoid  traffic  in  the  Chesapeake 
Bay.  A  plot  of  the  track  taken  for  underway  testing  on  the  ONR  YP  is  shown  in  Figure  7.3. 
Unfortunately,  the  constant  maneuvering  that  was  required  makes  the  data  nearly  impossible  to 
interpret. 


Figure  7,3:  Track  of  YP  Motion  for  Underway  Testing 

Because  the  steering  system  NILM  installation  on  the  SENECA  does  not  have  a  dedicated 
computer,  no  underway  data  has  been  taken  for  that  system.  If  data  were  to  be  collected  automat¬ 
ically,  it  would  have  to  be  accompanied  by  a  log  of  corresponding  rudder  motion.  Although  this 
log  should  be  maintained  anyway  by  the  ship’s  quartermasters,  it  might  be  difficult  to  reconcile 
it  with  power  data  from  the  NILM  if  the  two  clocks  used  are  out  of  agreement.  Therefore,  the 
best  way  to  collect  data  on  the  SENECA  is  to  put  someone  on  the  ship  for  controlled  underway 
testing  whose  only  duties  are  to  collect  NILM  data  while  recording  rudder  motions. 

7.2  Future  Systems 

Although  there  is  more  research  that  can  be  done  on  the  systems  featured  in  this  thesis,  data  has 
been  collected  for  other  systems  which  has  yet  to  be  analyzed.  There  are  also  some  promising 
systems  for  which  very  little  data  has  been  taken. 
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7.2.1  Anchor  Windlass 


The  anchor  windlass  consists  of  a  20  hp  variable  speed  drive  connected  through  a  reduction  gear 
to  a  capstan  and  “wildcat” .  A  wildcat  is  a  geared  drum  which  is  used  to  pull  up  the  anchor  while 
a  capstan  is  a  taller,  smooth  drum  that  is  located  on  the  same  axis  as  the  wildcat  and  is  used  in 
conjunction  with  lines  to  moor  the  ship.  There  is  a  clutch  on  the  wildcat  to  disengage  it  from  the 
capstan  so  that  the  two  can  move  independently. 

tc  It  was  reported  by  the  Engineer  Officer  of  the  SENECA  that  the  brake  for  the  starboard  anchor 
windlass  would  sometimes  rub  during  operation.  The  anchor  windlass  was  therefore  tested  for 
•  various  combinations  of  capstan  and  wildcat  engagement.  A  Tektronix  3054  oscilloscope  was  used 
to  measure  current  to  the  variable  speed  drive.  A  plot  of  currents  for  fast  speed  operation  is  shown 
in  Figure  7.4.  For  the  case  when  the  anchor  has  come  off  the  bottom  and  is  being  pulled  up  by 
the  wildcat,  there  is  noticeable  increase  in  current  drawn  by  the  motor. 


Figure  7.4:  Anchor  Windlass  Currents  for  Different  Loading  Conditions 

Installing  a  NILM  on  the  anchor  windlass  could  provide  valuable  information  about  the  power 
signatures  of  motors  tasked  with  moving  large  weights. 

7.2.2  Roll  Stabilization  Fin 

The  SENECA  has  fins  at  the  forward  end  of  the  hull  on  the  port  and  starboard  sides  that  are 
tasked  with  minimizing  the  effects  of  sea  state  on  the  roll  motion  of  the  ship.  Actuation  of  the 
fins  is  hydraulic  with  30  hp  motors  powering  pumps  that  maintain  pressure  in  the  system.  With 
the  fins  in  automatic,  pressure  transducers  that  sense  wave  motion  are  used  as  the  control  input 
to  position  the  fins.  Actuation  speed  of  the  fins  is  variable  and  fin  motion  can  be  controlled 
manually. 

It  is  possible  that  the  power  consumed  by  the  hydraulic  pumps  that  maintain  pressure  could 
be  used  as  an  indicator  of  fin  position.  It  might  even  be  possible  to  use  this  power  data  as  an 
indication  of  sea  state.  An  example  of  power  transients  for  fast  and  slow  speed  motion  with 
deflections  of  10  degrees  is  shown  in  Figure  7.5. 


87 


16 


.  “Speed" 10 


1000  2000  3000  4000  5000  6000  7000 


Figure  7.5:  Roll  Stabilizing  Fin  Power  Transients  for  Two  Speeds 


7.2.3  Propeller  Pitch  Control 

Yet  another  hydraulic  system  on  the  SENECA  is  the  system  that  maintains  pressure  for  the 
propeller  pitch  control  The  angle  of  attack  of  the  propeller  blades  is  adjusted  by  this  system  so 
that  the  main  engines  can  be  operated  at  more  efficient  rpm. 

There  is  a  NILM  installation  that  has  collected  underway  data  for  the  pumps  that  maintain 
system  pressure.  It  is  hoped  that  future  analysis  of  this  data  (along  with  controlled  testing)  may 
provide  indicators  of  shaft  speed  based  on  loading.  The  NILM  data  might  also  be  useful  as  a 
backup  indication  of  propeller  pitch  angle 

7.3  Conclusions 

There  are  too  many  shipboard  systems  which  can  be  monitored  by  the  NILM  to  list  them  all  in 
this  thesis.  The  results  presented  previously  are  all  for  systems  that  are  ultimately  mechanical 
in  nature.  The  application  of  the  NILM  to  more  “purely”  electrical  systems  (e.g.  radar  or 
communications)  has  yet  to  be  explored. 

The  results  presented  in  previous  chapters  show  the  suitability  of  the  NILM  for  monitoring 
electrical  loads  on  ships.  It  has  been  shown  that  it  is  capable  of  providing  backup  indications  of 
system  performance,  trending  equipment  performance,  and  recognizing  different  fault  conditions. 

Further  development  of  hardware  and  software  along  with  continued  research  into  the  behavior 
of  shipboard  systems  will  make  a  highly  capable  NILM-based  automated  monitoring  system  a 
reality  in  the  near  future. 
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Appendix  A 

Transducer  Specifications 


The  information  provided  on  the  following  pages  provides  technical  specifications  for  the  voltage 
and  current  transducers  typically  used  for  shipboard  NILM  installations. 
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Voltage  Transducer  LV  25-P 

For  the  electronic  measurement  of  voltages :  DC,  AC,  pulsed..., 
with  a  galvanic  isolation  between  the  primary  circuit  (high  voltage) 
and  the  secondary  circuit  (electronic  circuit). 


Electrical  data 


Primary  nominal  r.m.s.  current 

10 

mA 

u 

Primary  current,  measuring  range 

0  ..  ±  14 

mA 

K 

Measuring  resistance 

£> 

Minin 

p 

with  ±  12  V  @  ±  10  mA^ 

30 

190 

Q 

@±14mA,„„ 

30 

100 

a 

With  ±  15  V  <§>  ±  10  rrsAMS 

100 

350 

a 

100 

190 

a 

Secondary  nominal  r.m.s.  current 

25 

mA 

Conversion  ratio 

2500: 

1000 

Vc 

Supply  voltage  (±  5  %) 

±12.. 

15 

V 

'c 

Current  consumption 

10(@±15V)+I 

mA 

va 

R.m.s.  voltage  for  AC  isolation  test J),  50  Hz,  1  mn 

2.5 

kV 

Accuracy  -  Dynamic  performance  data 

X0 

Overall  Accuracy  @  lPK ,  TA  ~  25*C  @±12..  15V 

±0.9 

% 

@  ±  1 5  V  (±  5  %) 

±0.8 

% 

Linearity 

<0.2 

% 

Typ 

Max 

■o 

Offset  current  @  lp  «  0,  TA  ~  25*0 

±0.15 

mA 

'or 

Thermal  drift ofl0  0*C  ..  +  25*C 

±0.06 

±0.25 

mA 

+  25SC..+;70*C 

±0.10 

±0.35 

mA 

tr 

Response  time*5  @  90  %  of  VP{y)^ 

40 

MS 

Genera!  data 

T* 

Ambient  operating  temperature 

0 +  70 

°C 

Ts 

Ambient  storage  temperature 

-25 

+  85 

*C 

R, 

Primary  coil  resistance  @  TA  ~  70° C 

250 

a 

Rs 

Secondary  coil  resistance  @  TA  »  70*C 

110 

a 

m 

Mass 

22 

Q 

Standards  35  EN  50178 


Notes  :  Between  primary  and  secondary 

*  R «  «  25  kd  (UR  constant,  produced  by  the  resistance  and  Inductance 
of  the  primary  circuit) 

A  list  of  corresponding  tests  is  available 


LEM  Components 


lPN  =  10  mA 
VPN  =  10  ..  500  V 


Features 

*  Closed  loop  (compensated)  voltage 
transducer  using  the  Hal!  effect 

*  Insulated  plastic  case  recognized 
according  to  UL  94-V0. 

Principle  of  use 

*  For  voltage  measurements,  a  current 
proportional  to  the  measured  voltage 
must  be  passed  through  an  external 
resistor  R1  which  is  selected  by  the 
user  and  installed  in  senes  with  the 
primary  circuit  of  the  transducer. 

Advantages 

*  Excellent  accuracy 

*  Very  good  linearity 

*  Low  thermal  drift 

*  Low  response  time 

*  High  bandwidth 

*  High  immunity  to  external 
interference 

*  Low  disturbance  in  common  mode. 

Applications 

*  AC  variable  speed  drives  and  servo 
motor  drives 

*  Static  converters  for  DC  motor  drives 

*  Battery  supplied  applications 

*  Uninterruptible  Power  Supplies 
(UPS) 

*  Power  supplies  for  welding 
applications. 
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Mechanical  characteristics  Remarks 

*  General  tolerance  ±  0.2  mm  *  ls  is  positive  when  Vp  is  applied  oh  terminaf  +HT. 

*  Fastening  &  connection  of  primary  2  pins  ♦  This  is  a  standard  model.  For  different  versions  (supply 

0.835  x  0.635  mm  voltages*  turns  ratios,  unidirectional  measurements...}* 

*  Fastening  &  connection  of  secondary  3  pins  0  1  mm  please  contact  us. 

*  Recommended  FOB  hole  1.2  mm  :  .  ... 


Instructions  for  use  of  the  voltage  transducer  model  LV  25-P 

Primary  resistor  R 1 :  the  transducer's  optimum  accuracy  is  obtained  at  the  nominal  primary  current.  As  far  as  possible,  Ui  should  be 
calculated  so  that  the  nominal  voltage  to  be  measured  corresponds  to  a  primary  current  of  10  mA. 

Example:  Voltage  to  be  measured  -250  V  a)  R 1 =  25  kQ/2.5  W,  lp  =  10mA  Accuracy  =  ±0.8%  of  {%  JA  *'+ 25*C) 

b)Rt  =  50  mi  1.25  W*  ip  *  5  mA  Accuracy  *  ±1.8  %  of  (@  TA  =  +25*0) 

Operating  range  (recommended) :  taking  Into  account  the  resistance  of  the  primary  windings  {which  must  remain  low  compared  to  H,  in  order 
to  keep  thermal  deviation  as  low  as  possible)  and  the  isolation*  this  transducer  is  suitable  for  measuring  nominal  voltages  from  10  to  500  V. 


LEM  reserves  the  right  to  carry  out  modifications  on  its  transducers*  in  order  to  improve  them*  without  previous  notice. 


Current  Transducer  LA  55-P 


For  the  electronic  measurement  of  currents  :  DC,  AC,  pulsed..., 
with  a  galvanic  isolation  between  the  primary  circuit  (high  power) 
and  the  secondary  circuit  (electronic  circuit). 


Electrical  data 

w 

Primary  nominal  r.m.s.  current 

50 

A 

'p 

Primary  current,  measuring  range 

0 

±70 

A 

R„ 

Measuring  resistance  @ 

Ta  =  70°C 

Ta  *  85°C 

with  ±  12  V  @±50Atmk 

10  100 

60  95 

n 

@±70Am 

10  50 

601*  60=5 

a 

with  ±  15  V  @  ±  50  A 

50  160 

135  155 

a 

@±70A_ 

50  90 

135?>135^ 

a 

Secondary  nominal  r.m.s.  current 

50 

mA 

Conversion  ratio 

1: 

1000 

i 

vc 

Supply  voltage  (±  5  %) 

±12  ..  15 

V 

'o 

Current  consumption 

10(@±15V}+1. 

mA 

v. 

R.m.s,  voltage  for  AC  isolation  test,  50  Hz, 

1  mn  2.5 

kV 

Accuracy  -  Dynamic  performance  data 


X 

Accuracy  @  lps, ,  TA  “  25*0  @±15V(±5%) 

±0.65 

% 

@±12  ..15V  (±5  %) 

±  0.90 

% 

e. 

linearity 

<0.15 

% 

Typ- 

Max 

«o 

Offset  current  @  ~  0,  TA  ~  25°C 

±0.2 

mA 

Residual  current 33  @  lp  »  0,  after  an  overload  of  3  x  lpN 

±0.3 

mA 

»OT 

Thermal  drift  of  l0  0°C  70°C 

±0.1 

±0.5 

mA 

-  25°C  ..  +  85°C  : 

±0.1 

±0.6 

mA 

Reaction  time  @  10  %  of  lP{r>35f 

<500 

ns 

Response  time  @  90  %  of  fF 

<1 

ps 

dl/dt  di/dt  accurately  followed 

>200 

A/ps 

f 

Frequency  bandwidth  (-  1  dB) 

DC  >.  200 

kHz 

General  data 

Ta 

Ambient  operating  temperature 

-25 .. 

+  85 

*C 

Ts 

Ambient  storage  temperature 

-40.. 

+  90 

Secondary  coil  resistance  @  TA~706C 

80 

a 

*  tA*85’C 

85 

a 

m 

Mass 

18 

9 

Standards4' 

EN  50178 

Notes  :  Measuring  range  limited  to  ±  60  Amaj( 

Measuring  range  limited  to  ±  55 
Result  of  the  coercive  field  of  the  magnetic  circuit 
^  A  list  of  corresponding  tests  is  avaiiabie 


LEM  Components 


U  =  50  A 


Features 

*  Closed  loop  (compensated)  current 
transducer  using  the  Hall  effect 

♦  Printed  circuit  board  mounting 

♦  Insulated  plastic  case  recognized 
according  to  UL  94-VO. 

Advantages 

♦  Excellent  accuracy 

•  Very  good  linearity 

•  Low  temperature  drift 

»  Optimized  response  time 

•  Wide  frequency  bandwidth 

*  No  insertion  losses 

•  High  immunity  to  external 
interference 

♦  Current  overload  capability. 

Applications 

•  AC  variable  speed  drives  and  servo 
motor  drives 

*  Static  converters  for  DC  motor  drives 

*  Battery  supplied  applications 

*  Uninterruptible  Power  Supplies 

/.  (UPS)  '"/"i  ■ 

•  Switched  Mode  Power  Supplies 

(SMPS)  :  .. 

•  Power  Supplies  for  welding 
applications. 


980706/8 


www.lem.com 
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Dimensions  LA  55-P  (in  mm.  1  mm  =  0.0394  inch) 


Bottom  view 


{ 
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i 
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-  5,08  22.86 
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Front  view 


\  V 

Year  Week 


Left  view 
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. 

■ 

l . .  •  •  L. . Jj 

- . . .  . :  . J 

■ 

ON 

:  K 

_ :■* 

in 

m 

to 

27.6  , 

4.5  +/-0.3 

Secondary  terminals 

T armina!  +  :  supply  voltage  +  12..  15  V 
Terminal  - :  supply  voltage  *  12 15  V 
Terminal  M  :  measure 


Connection 


LA  55-P 


■  i  Rh 

. -4Z3 - oOV 


Mechanical  characteristics 


Remarks 


*  General  Iterance 

*  Primary  through-hole 

*  Fastening  &  connection  of  secondary 

Recommended  PCB  hole 


±  0.2  mm  *  ls  is  positive  when  lp  flows  in  the  direction  of  the  arrow. 

12.7  x  7  mm  *  Temperature  of  the  primary  conductor  should  not  exceed 

3  pins  ;  90SC. 

0.63  x  0.56mm  *  Dynamic  performances  (di/dt  and  response  time)  are  best 

0.9  mm  With  a  single  bar  completely  filling  the  primary  hole. 

•  In  order  to  achieve  the  best  magnetic  coupling,  tile  primary 
windings  have  to  be  wound  over  the  top  edge  of  the  device, 

•  This  is  a  standard  model  For  different  versions  (supply 
voltages,  turns  ratios,  unidirectional  measurements...}, 
please  contact  us. 


LEM  reserves  the  right  to  carry  out  modifications  on  its  transducers,  in  order  to  improve  them,  without  previous  notice. 
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Current  Transducer  LA  205-S 

For  the  electronic  measurement  of  currents  ;  DC,  AC,  pulsed.,, 
with  a  galvanic  isolation  between  the  primary  circuit  (high  power) 
and  the  secondary  circuit  (electronic  circuit). 


PN 


200  A 


Electrical  data 


Ip* 

Primary  nominal  r.m.s;  current 

200  : 

A 

Features 

ip 

Primary  current,  measuring  range 

0.. 

±300 

A 

*P  max 

Measuring  overload 

600 

A 

♦  Closed  loop  (compensated)  current 

R* 

Measuring  resistance  @ 

ta= 

70°C 

1  ta  = 

85*  C 

transducer  using  the  Hall  effect 

p 

"  M 

i  max 

rrtsrt  rtiaic 

*  Insulated  plastic  case  recognized 

with  ±  12  V 

#±20GA,m 

0 

68 

0 

66 

0 

according  to  UL  94-V0 

#±300 

0 

33 

0 

30 

0 

•  Patent  pending. 

with  ±  15  V 

@±200Afr- 

5 

95 

5 

93 

0 

'«±300'A^ 

5 

50 

5 

49 

0 

Advantages 

fSN  Secondary  nominal  r.m.s.  current 

Conversion  ratio 
Vc  Supply  voltage  (±  5  %) 

lc  Current  consumption 

V&  R.m.s  rated  voltage 2),  safe  separation 

basic  isolation 


100  mA 

1:2000 

±12  ..15  V 

2dtf±l5V)  +  ls  mA 
162$  V 

3250  V 


Accuracy  -  Dynamic  performance  data 


XG 

.  Overall  accuracy  @  1PN  ,  TA  -  25*C 

±0.8 

% 

e. 

Linearity 

<0.1 

% 

Typ 

Max 

lo 

Offset  current  @  L  -  0,  T.  =  25°C 

■ !  ■. 

±0.15 

rriA 

'OM 

Residual  current 3>@  1  «  0,  after  an  overload  of  3  x  L.j 

±0.50 

mA 

•or 

Thermal  drift  of  lG 

-  10aC  ..  +  85°C 

±0.15 

±0.30 

mA 

L 

"Reaction  time  @  10  %  of  lPmax 

<500 

ns 

tf 

Response  time4)  @  90  %  of  lp?rai< 

<1 

Ms 

di/dt 

di/dt  accurately  followed 

>  100 

A/ps 

f 

Frequency  bandwidth  {-  3  dB) 

DC  ..100 

kHz 

General  data 

Ta 

Ambierit  operating  temperature 

-id.: 

+  85 

•C 

Ambient  storage  temperature 

-40.. 

+  90 

Ac 

Rs 

•  Secondary  coil  resistance  @ 

TA  =  70°C 

35 

a 

:  Ta  =  85»C 

37 

0 

m 

Mass 

110 

g 

Standards  5' 

h  :  ■  ,  i  -7  -  ' 

EN  50178 

•  Excellent  accuracy 

♦  Very  good  linearity 

•  Low  temperature  drift 

•  Optimized  response  time 

•  Wide  frequency  bandwidth 

*  No  insertion  losses 

•  High  immunity  to  external 
interference 

*  Current  overload  capability. 

Applications 

♦  AC  variable  speed  drives  and  servo 
motor  drives 

•  Static  converters  for  DC  motor  drives 

♦  Battery  supplied  applications 

♦  Uninterruptible  Power  Supplies 

;  {UPS)  . 

♦  Switched  Mode  Power  Supplies 
(SMPS) 

•  Power  supplies  for  welding 
applications. 


Notes  :  ^  3  mn/hour  @  Vc  =  ±  1 5  V,  R  5  0 

2>  Pollution  class  nr  2.  With  a  non  insulated  primary  bar  v/hich  fills  the 
through-hole 

3)  The  result  of  the  coercive  field  of  the  magnetic  circuit 

4)  With  a  di/dt  of  100  Af\i$ 

A  list  of  corresponding  tests  is  available 


980716/5 


LEM  Components 


www.lem.com 
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Dimensions  LA  205-S  {in  mm.  1  mm  =  0.0394  inch) 
Front  view 


HGLEX  5046-04/AG 


Top  view 

Mechanical  characteristics 

*  General  tolerance 

*  Fastening 

*  Primary  through-hole 

*  Connection  of  secondary 


±  0.5  mm 
2  holes  0  5.5  mm 
23  x  18  mm 
Molex  5048-04/ AG 


Secondary  terminals 

Pin  1  :  supply  Voltage  +12  „  15  V 
Pin  2  ;  measure 
Pin  3  ;  supply  voltage  -12 15  V 
Pin  4  :  NC 


Connection 


LA  205-S 


Remarks 

*  fs  is  positive  When  lp  flows  In  the  direction  of  the  arrow. 

*  Temperature  of  the  primary  conductor  should  not  exceed 

100°C.  /  >-*;  ; V ;/ I  -  • 

*  Dynamic  performances  (dl/dt  and  response  time)  are  best 
with  a  single  bar  completely  filling  the  primary  hole. 

*  This  is  a  standard  model.  For  different  versions  {supply 

voltages,  turns  ratios,  unidirectional  measurements...), 
please  contact  us.  \4"  •  '  =  :'  :  .  . 


LEM  reserves  the  right  to  carry  out  modifications  on  its  transducers,  in  order  to  Improve  them,  without  previous  notice. 
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Current  Transducer  LA  305-S 

For  the  electronic  measurement  of  currents  :  DC,  AC,  pulsed. i, 
with  a  galvanic  isolation  between  the  primary  circuit  (high  power) 
and  the  secondary  circuit  (electronic  circuit). 


Electrical  data 

Primary  nominal  r.m.s.  current 

800 

A 

K 

Primary  current,  measuring  range 

0.. 

±500 

A 

K 

Measuring  resistance  @ 

Ta  = 

70°C 

Ta  =  85“C 

D 

n,Mmsx 

with  ±  12  V  @  ±  300  A  inax 

0 

52 

0  50 

a 

@±500ArMX 

0 

17 

0  15 

a 

with  ±  15  V  @  ±  300  A  m 

0 

75 

6  73 

a 

@  ±  500  A 

0 

31 

5  29 

a 

Secondary  nominal  r.m.s.  current 

120 

mA 

Conversion  ratio 

1  : 

2500 

vc 

Supply  voltage  (±  5  %) 

±12  ..  15 

V 

«c 

Current  consumption 

20{g±15V)  +  ls 

mA 

vb 

R.m.s.  rated  voltage  1),  safe  separation 

1750 

V 

basic  isolation 

3500 

V 

Accuracy  -  Dynamic  performance  data 


Overall  accuracy  @  lPN  ,  Ta  ~25X 

±0.8 

% 

e. 

Linearity 

<0.1 

% 

Typ 

Max 

»o 

Offset  current  @  L  =  0,  T.  ~  25X 

±0.20 

mA 

Residual  current 75  @  ip-  0,  after  an  overload  of  3  x  lPN 

±  0.40 

mA 

'OT 

Thermal  drift  of  !c 

-10X..  +  85X 

±  0.12 

±0.30 

mA 

Reaction  time  @10%  of  lPfr3X 

<500 

ns 

t 

Response  time3}  @  90  %  of  lPrrax 

<  1 

Ms 

di/dt  di/dt  accurately  followed 

>100 

A/ps 

f 

Frequency  bandwidth  (-  3  dB) 

DC..  100 

kHz 

General  data 

Ta 

Ambient  operating  temperature 

-10.. 

±85 

X 

T5 

Ambient  storage  temperature 

-40 

+  90 

X 

Secondary  coil  resistance  @ 

Ta  =  70“C 

35 

a 

Ta  -  85X 

37 

a 

m 

Mass 

200 

g 

Standards 

EN  50178 

Notes  :  Pollution  class  2.  With  a  non  insulated  primary  bar  which  fills  the 
through-hole 

&  The  result  of  the  coercive  field  of  the  magnetic  circuit 
3)  With  a  di/dt  of  100  A/p$ 

4'  A  list  of  corresponding  tests  is  available 


LEM  Components 


Ip*  =  300  A 


Features 

*  Closed  loop  (compensated)  current 
transducer  using  the  Hall  effect 

•  Insulated  plastic  case  recognized 
according  to  UL.  94-VG 

*  Copyright  protected. 

Advantages 

•  Excellent  accuracy 

•  Very  good  linearity 

♦  Low  temperature  drift 

*  Optimized  response  time 

♦  Wide  frequency  bandwidth 

♦  No  Insertion  losses 

*  High  immunity  to  external 
interference 

♦  Current  overload  capability. 

Applications 

•  AC  variable  speed  drives  and  servo 
motor  drives 

•  Static  converters  for  DC  motor  drives 

•  Battery  supplied  applications 

*  Uninterruptible  Power  Supplies 
(UPS) 

•  Switched  Mode  Power  Supplies 

.  .  ,{SMPS)  . 

*  Power  supplies  for  welding 
applications. 


980707/6 
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Dimensions  LA  305-S  (in  mm.  I  mm  =  0.0394  inch) 


Front  view 


4~f~ 

— 


Left  view 


MOLEX  S046-S4/AS 


Secondary  terminals 


Pin  1  :  supply  voltage  *  12  ,.15  V 
Pin  2  :  measure 
Pin  3  :  supply  voltage  - 12  ..  15  V 
Pin  4  :  NC 


Connection 


2l%- 


Top  view 


Mechanical  characteristics 


Remarks 


*  General  tolerance 

*  Fastening 

*  Primary  through-hole 

*  Connection  of  secondary 


±  0.5  mm 
2  holes  0  5,5  mm 
25.5  x  25,5  mm 
Molex  5046-04/AG 


■  fs  is  positive  when  lF  flows  in  the  direction  of  the  arrow 
•  Temperature  of  the  primary  conductor  should  not  exceed 
100*C.  ;  ■  ;•  '■  -  '  '  -  : 

‘  Dynamic  performances  (dt/dl  and  response  time)  are  best 
with  a  single  bar  completely  filling  the  primary  hole. 

;  This  Is  a  standard  model.  For  different  versions  (supply 
voltages,  turns  ratios,  unidirectional  measurements,..), 
please  contact  us. 


LEM  reserves  the  right  to  carry  out  modifications  on  its  transducers,  in  order  to  improve  them,  without  previous  notice. 
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Appendix  B 

NILM  Hardware  Listing 


The  purpose  of  this  appendix  is  to  provide  construction  information  for  individuals  performing 
future  work  with  the  NILM. 


Table  B.l:  Electrical  Components  Required  for  Voltage  Sens¬ 
ing  Board 


Part 

Value 

Qty. 

Supplier 

Part  No. 

Capacitor 

0.1  \xF  ceramic 

2 

Various 

N/A 

l/j,F  ceramic 

2 

10 fiF  electrolytic 

2 

Resistor 

ioo  m 

9 

Digikey 

BC100KW-3JCT 

470 kCl 

9 

Digikey 

BC470KW-3JCT 

Transducer 

N/A 

3 

LEM 

LV  25-P 

Fuseholder 

N/A 

1 

Littlefuse 

L60060M-3PQ 

Fuse 

100  mA/600  V 

3 

Littlefuse 

KLKR1/10 

Wire  and  Crimp  Lug 

12  AWG 

3 

Various 

N/A 

Straight  Term.  Block 

2-pin 

1 

Digikey 

ED2354-ND 

Diag.  Term.  Block 

2-pin 

4 

Digikey 

281-1421-ND 

Table  B.2:  Additional  Electrical  Components  Required  for 
NILM  Setup 


Part 

Value 

Qty. 

Supplier 

Part  No. 

Power  Supply 

±12V 

1 

Jameco 

148689 

Wiring  Terminal 

1 

Advantech 

ADAM-3968 

A/D  Converter 

PCI 

1 

Advantech 

PCI-1710 

SCSI-II  Cable 

68-pin 

1 

Various 

N/A 

Power  Cable 

120V 

1 

Various 

N/A 
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Table  B.3:  Hardware  Required  for  NILM  Box 


Part 

Value 

Qty. 

Purpose 

NEMA  Box 

12”xl0”x5” 

House  voltage  board, 
terminal  board,  and 
power  supply 

Backplane 

12”xl0”x5” 

1 

Mounting  surface  for 
voltage  board 

Phillips  Machine  Screw 

■Mil 

6 

Mount  voltge  board 
and  junction  strip  to 
backplane 

Phillips  Machine  Screw 

10x32,  0.5”  length 

2 

Mount  DIN  rail  to  box 

Spacer 

10x32,  smooth  inte¬ 
rior,  0.25”  length,  alu¬ 
minum 

4 

Between  voltage  board 
and  backplane 

Lock  Washer  and  Nut 

10x32 

H 

For  voltage  board  and 
junction  strip 

Spacer 

6x32,  threaded,  0.5” 
length,  aluminum 

H 

Mount  power  supply  to 
side  of  NEMA  box 

Lock  Washer  and  Nut 

6x32 

4 

Secure  spacers  to 
NEMA  box 

Phillips  Machine  Screw 

6x32,  0.375”  length 

4 

Attaching  power  sup¬ 
ply  to  spacers 

Junction  Strip 

3-terminal 

■ 

Provide  power  to  cur¬ 
rent  transducer 

Molex  Connector 

3-pin 

■ 

Power  cord  connection 
to  power  supply 

Molex  Connector 

4-pin 

1 

Power  supply  connec¬ 
tion  to  voltage  board 
and  junction  strip 

Wire  Clamps 

0.375”  “Halex” 

3 

Strain  relief  for  volt¬ 
age  input,  current  in¬ 
put,  and  power  cord 

DIN  Rail 

10”  length 

1 

Mounting  for  SCSI-II 
wiring  terminal  board 
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Appendix  C 

Sewage  System  MATLAB  Scripts 


The  following  script  was  written  to  trend  vacuum  pump  cycling  in  the  sewage  system  of  the 
USCGC  SENECA.  The  input  vector  A  is  power  data  sampled  at  120  Hz  for  the  sewage  system. 
The  output  is  the  number  of  times  that  one  or  both  of  the  pumps  was  running  along  with  the 
average  time  of  each  pump  run. 

j=0; 

z=ones (size (y) ) ; 
for  i=l : length (y) ; 
if  y(i)>500 
z(i)=l; 
else  z(i)=0; 
end 

end 

z; 

if  z(l)==l 

j=j+i; 

end 

for  i=2 : length (z) ; 
if  z(i)==l 

if  z(i-l)==0 

j=j+i; 

end 

end 

end 

runs=j 

average_run_time=sum(z)/(j*120) 
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The  following  script  simulates  operation  of  the  sewage  system  vacuum  pumps  onboard  the 
USCGC  SENECA.  The  average  time  between  system  usages  is  and  amount  of  vacuum  leak  is 
input  with  number  pump  runs  and  average  run  time  per  run  output. 

f unct ion  [runs , runt ime] =f lushs im (t ime_f lush , Pdot _leak) 

n=60;  %  simulation  time  (min);  ' 

k=200;  %  number  of  points  per  minute  in  simulation 

P_f lush=l ;  %  pressure  decrease  caused  by  toilet  flush  (in  Hg) 

Pdot_pumpl=19;  %  vacuum  rate  gained  by  pumpl  (in  Hg/min) 

Pdot_pump2=17;  %  vacuum  rate  gained  by  pump2  (in  Hg/min) 

Plower_setpt=7 . 5 ;  %  vacuum  at  which  two  pumps  come  on  (in  Hg) 

Plow_setpt=10;  %  vacuum  at  which  pump  comes  on  (in  Hg) 

Phigh_setpt=15;  %  vacuum  at  which  pump  shuts  off  (in  Hg) 

P(l)=15; 

Pump_status=0; 

Twopump_status=0 ; 
f lush_count=0 ; 
j=0; 

for  i-2:n*k; 

if  P(i-l)  <  Plow_setpt 
Pump_status=l ; 
elseif  P(i-l)  >  Phigh^setpt 
Pump_status=0 ; 

Tw0pump_status=O ; 

end 

if  P(i-l)  <  Plower_setpt 
Twopump_status=l ; 

end 

if  Twopump_status==l 
Pump_status=0; 

end 

Pump_Counter(i)=Pump_status  +  Twopump_status ; 
if  Pump_Counter (i)==l 

if  Pump _ Counter ( i - 1 ) ==0 

v  \ 

end 

end 

if  flush_count/k  <  time_flush 
f lushprob=(f lush_eount/k) "2/ (2*time_f lush~2) ; 
else  , 

f lushprob=l- (  (2*t ime_f lush-  (f  lush_count/k)  )  "2/  (2*t ime_f  lush~2)  )  ; 
end 

f lushchance=rand; 

if  flushchance  <  flushprob 
dP_f lush=P_f lush ; 
f lush_count=0 ; 

else 

f lush_count=f lush_count+i ; 
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dP_f lush=0; 


end 

if  rem(j,2)==0 

Pdot_pump=Pdot_pump2 ; 

else 

Pdot_pump=Pdot_pumpl; 

end 

P(i)=P(i-l) -dP_f lush- (Pdot_leak) * (1/k)  +  Pump_status*Pdot_pump* ( 1/k) 
+  Twopump_status*(Pdot_pumpl  +  Pdot_pump2)*(l/k) ; 

end 

runs=j ; 

runtime=sum(Pump_Counter)*60/(j*k)  ; 
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Appendix  D 

Sewage  System  Leak  Data 

The  following  data  was  taken  onboard  the  USCGC  SENECA  between  February  20  and  March  6, 
2004: 


Leak  Data  for  Leak  Rate  of  0  SCFH 


Hour 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

1800-1900 

15 

23.78 

356.7 

1900-2000 

io 

26.14 

261.4 

2000-2100 

10 

19.56 

195.6 

2100-2200 

6 

17.05 

102.3 

2200-2300 

7 

18.47 

129.29 

2300-0000 

8 

21.62 

172.96 

0000-0100 

15 

23.99 

359.85 

0100-0200 

7 

18.78 

131.46  | 

0200-0300 

6 

16.49 

98.94 

0300-0400 

7 

17.12 

119.84  •"] 

0400-0500 

8 

19.16 

153.28 

0500-0600 

6 

17.24 

103.44 

0600-0700 

7 

17.11 

119.77 

0700-0800 

13 

21.39 

278.07 

0800-0900 

io 

19.57 

195.7 

0900-1000 

8 

19.63 

157.04 

1000-1100 

10 

18.62 

186.2 

1100-1200 

10 

22.28 

222.8 

1200-1300 

11 

18.41 

202.51 

1300-1400 

9 

21.79 

196.11 

1400-1500 

10 

17.47 

174.7 

1500-1600 

8 

17.18 

i37.44 
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Leak  Data  for  Leak  Rate  of  42  SCPH 


Hour 

Runs 

Time  per  run  (seconds) 

12 

17.01 

204.12 

1900-2000 

12 

22.32 

267.84 

2000-2100 

9 

19.89 

179.01 

2100-2200 

9 

19.88 

178.92 

2200-2300 

8 

21.28 

170.24 

2300-0000 

8 

24.88 

199.04 

0000-0100 

7 

15.9 

111.3 

0100-0200 

7 

13.18 

92.26 

0200-0300 

7 

15.41 

107.87 

0300-0400 

10 

15.71 

157.1 

0400-0500 

mnm 

15.45 

108.15 

0500-0600 

mm 

18.85 

207.35 

0600-0700 

mm 

20.33 

345.61 

0700-0800 

24 

20.46 

491.04 

0800-0900 

15 

20.56 

308.4 

0900-1000 

24 

21.94 

526.56 

1000-1100 

14 

24.7 

345.8 

1100-1200 

17 

27.1 

460.7 

1200-1300 

19 

24.75 

470.25 

1300-1400 

14 

27.13 

379.82 

IKH 

12 

23.02 

276.24 

1500-1600 

10 

18.98 

189.8 

Leak  Data  for  Leak  Rate  of  51  SCFH 


Hour 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

1700-1800 

15 

26.59 

398.85 

16 

25.5 

408 

mu 

10 

22.05 

220.5 

2000-2100 

11 

21.15 

232.65 

2100-2200 

9 

15.55 

139.95 

2200-2300 

9 

16.82 

151.38  ! 

jtkirilri'iiil 

20.86 

208.6  : 

8 

17.05 

136.4 

8 

15.27 

122.16 

8 

17.7 

141.6 

9 

18.14 

163.26 

9 

15.78 

142.02 

0500-0600 

12 

20.51 

246.12 

iiMiliHilil 

22 

21.66 

476.52 

21 

21.33 

447.93 

18 

26.62 

479.16 

25 

24.62 

615.5 

21 

30.51 

640.71 

16 

22.27 

356.32 

1200-1300 

23 

22.15 

509.45 

1300-1400 

23 

23.65 

543.95 

1400-1500 

11 

25.23 

277.53 

1500-1600 

13 

20.11 

261.43 
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Leak  Data  for  Leak  Rate  of  62.5  SOPH 


Hour 

Time  per  run  (seconds) 

1700-1800 

13 

24.93 

324.09 

1800-1900 

16 

20.77 

332.32 

1900-2000 

12 

20.46 

245.52 

2000-2100 

10 

17.76 

177.6 

2100-2200 

13 

17.85 

232.05 

2200-2300 

14 

19.79 

277.06 

mmm 

10 

18.41 

184.1 

0000-0100 

10 

16.96 

169.6 

0100-0200 

9 

15.29 

137.61 

0200-0300 

10 

15.4 

154 

0300-0400 

11 

15.61 

171.71 

0400-0500 

9 

18.49 

166.41 

0500-0600 

14 

18.51 

259.14 

0600-0700 

23 

20.82 

478.86 

0700-0800 

13 

20.63 

268.19 

0800-0900 

13 

30.5 

396.5 

0900-1000 

20 

22.39 

447.8 

111^121111 

22 

22.94 

504.68 

19 

26.99 

512.81 

1200-1300 

22 

22.95 

504.9 

1300-1400 

18 

28.09 

505.62 

1400-1500 

18 

20.74 

373.32 

1500-1600 

15 

19.11 

286.65 
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Leak  Data  for  Leak  Rate  of  102.5  SCFH 


Hour 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

1700-1800 

16 

22.11 

353.76 

1800-1900 

11 

19.34 

212.74 

1900-2000 

12 

22.6 

271.2 

11 

17.35 

190.85 

iMTOITll 

10 

21.1 

211 

2200-2300 

11 

18.48 

203.28 

11 

18.28 

201.08 

10 

15.39 

153.9 

9 

18.13 

163.17 

9" 

16.88 

151.92  ; 

imsuTii 

11 

26.55 

292.05  ; 

9 

16.56 

149.04 

9 

17.65 

158.85 

9 

17.3 

155.7 

14 

20.3 

284.2 

14 

19.71 

275.94 

12 

22.32 

267.84 

11 

20.17 

221.87 

15 

20.03 

16 

17.56 

280.96 

10 

23.13 

231.3 

10 

20.27 

202.7 

1500-1600 

11 

21.61 

237.71 

Leak  Data  for  Leak  Rate  of  132.5  SCFH 


Hour 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

19 

25.32 

481.08 

16 

22.43 

358.88 

17 

22.4 

380.8 

16 

18.26 

292.16 

2100-2200 

15 

18.38 

275.7 

19 

16.78 

318.82 

23 

16.12 

370.76 

0000-0100 

21 

16.76 

351.96 

0100-0200 

22 

15.8 

347.6 

0200-0300 

18.5 

462.5 

0300-0400 

25.79 

618.96  ! 

0400-0500 

26 

15.37 

399.62 

lllrilllBlliillll 

28 

16.7 

467.6 

31 

20.14 

624.34 
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The  following  data  was  produced  by  the  MATLAB  simulation  script  f  lushsim: 


Simulated  Weekday  Data  for  Leak  Rate  of  0  in./minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

5 

m 

20.625 

247.5 

5 

bb 

19.96363636 

219.6 

5 

mm 

20.075 

240.9 

5 

12 

20.075 

240.9 

2 

21 

20.95714286 

440.1 

1 

29 

23.4 

678.6 

5 

12 

20.075 

240.9 

2 

22 

21.40909091 

471 

5 

12 

20.075 

240.9 

1 

32 

22.33125 

714.6 

2 

21 

20.92857143 

439.5 

5 

12 

20.075 

240.9 

5 

11 

19.96363636 

219.6 

5 

12 

20,075 

240.9 

7 

9 

19.96666667 

179.7 

7 

10 

20.07 

200.7 

7 

10 

18.3 

183 

5 

12 

20.075 

240.9 

5 

11 

19.96363636 

219.6 

7 

20.07 

200.7 

5 

20.075 

240.9 

5 

11 

20.26363636 

222.9 

2 

22.245 

444.9 

5 

11 

19.96363636 

219.6 
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Simulated  Weekday  Data  for  Leak  Rate  of  0.5  in./minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

5 

18 

18.28333333 

329.1 

5 

19 

18.53684211 

352.2 

5 

19 

18.37894737 

349.2 

5 

16 

19.40625 

310.5 

2 

26 

20.63076923 

536.4 

1 

36 

22.23333333 

5 

19 

18.93157895 

359.7 

2 

26 

21.20769231 

55114 

5 

18 

18.75 

337.5 

1 

37 

22.23243243 

822.6 

2 

26 

20.37692308 

529.8 

5 

18 

19.26666667 

346.8 

5 

17 

19.39411765 

329.7 

5 

17 

18.74117647 

318.6 

7 

16 

17.925 

286.8 

7 

15 

19.38 

290.7  | 

7 

16 

18.76875 

300.3 

5 

18 

19.06666667 

343.2 

5 

18 

17.78333333 

320.1 

7 

15 

19.5 

292.5 

5 

17 

19.05882353 

324 

5 

18 

18.76666667 

337.8 

2 

26 

20.81538462 

541.2 

5 

17 

18.9 

321.3 
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Simulated  Weekday  Data  for  Leak  Rate  of  1.0  in./minute 


Mean  Time 
Between 


Runs 


Time  per  run  (seconds) 


Total  runtime  (seconds) 


Flushes 

1 

5 

5 

5 

2 

1 

5 

2 

5 

1 

2 

5 

5 

5 

7 

7 

7 

5 

5 

7 

5 

5 

2 

5 


22 

23 

21 

23 

32 

42 

22 

29 

23 

42 

32 

22 

23 

23 
20 
20 

19 

24 
22 

20 
22 
22 
32 
22 


19.15909091 

19.46086957 

19.95714286 

18.61304348 

20.371875 

21.95 

19.2 

20.87586207 

18,88695652 

22.00714286 

20.08125 

19.54090909 

19.16086957 

19.6173913 

18.765 

18.735 

20.1 

18.9875 

19.77272727 

19.59 

19.30909091 

19.59545455 

19.903125 

19.62272727 


421.5 

447.6 

419.1 

428.1 

651.9 

921.9 

422.4 

605.4 

434.4 
924.3 

642.6 

429.9 

440.7 

451.2 

375.3 

374.7 

381.9 

455.7 
435 

391.8 

424.8 
431.1 

636.9 
431.7 
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Simulated  Weekday  Data  for  Leak  Rate  of  1.5  in./minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

5 

26 

20.27307692 

527.1  j 

5 

27 

20.1 

542.7 

5 

27 

19.67777778 

531.3 

5 

26 

20.50384615 

533.1 

2 

35 

20.76857143 

726.9 

1 

45 

23.27333333 

1047.3 

5 

26 

20.66538462 

537.3 

2 

34 

20.69117647 

703.5 

5 

26 

20.49230769 

532.8 

1 

45 

22.42666667 

1009.2 

2 

36 

21.2 

763.2 

5 

27 

20.48888889 

553.2 

5 

26 

20.86153846 

542.4 

5 

26 

19.91538462 

517.8 

7 

25 

19.62 

490.5 

7 

19.188 

479.7 

7 

19.644 

491.1 

5 

26 

20.43461538 

531.3 

5 

26 

20.43461538 

531.3 

7 

25 

19.728 

493.2 

5 

25 

20.772 

519.3 

5 

26 

20.79230769 

540.6 

2 

35 

20.72571429 

725.4 

5 

27 

20.61111111 

556.5 
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Simulated  Weekday  Data  for  Leak  Rate  of  2.0  in./minute 


Mean  Time 

Between 

Flushes 

Runs 

Total  runtime  (seconds) 

5 

31 

HHH 

638.1 

5 

31 

20.24516129 

627.6 

5 

31 

20.43870968 

633.6 

5 

33 

20.00909091 

660.3 

2 

38 

21.77368421 

827.4 

1 

45 

24.87333333 

1119.3 

5 

31 

20.53548387 

636.6 

2 

39 

22.12307692 

862.8 

5 

na 

20.55483871 

637.2 

1 

n 

24.28695652 

1117.2 

2 

41 

21.15365854 

867.3 

5 

31 

,  20.64193548 

639.9 

5 

32 

20.4375 

654 

5 

31 

20.48709677 

635.1  ! 

7 

29 

20.37931034 

591 

7 

29 

20.07931034 

582.3 

7 

30 

19.87 

596.1 

5 

32 

20.596875 

659.1 

5 

HB 

19.88709677 

616.5  j 

7 

KM 

20.35862069 

590.4 

5 

20.990625 

671.7 

5 

19.82903226 

614.7 

2 

22.29473684 

847.2 

5 

20.79677419 

644,7 

Simulated  Weekday  Data  for  Leak  Rate  of  2.5  in./minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

5 

35 

21.20571429 

742.2  j 

5 

35 

21.23142857 

743.1 

5 

35 

21.06 

737.1 

5 

34 

21.21176471 

721.2 

2 

41 

22.70487805 

930.9 

1 

49 

24.97346939 

1223.7 

‘5 

33 

21.61818182 

713.4 

2 

42 

22.60714286 

949.5 

5 

21.30857143 

745.8 

1 

47 

25.85744681 

1215.3 

2 

41 

23.2902439 

954.9 

5 

35 

20.82 

728.7 

5 

35 

21.58285714 

755.4 

5 

35 

21.36 

747.6  1 

7 

32 

21.95625 

702.6 

7 

21.35625 

683.4 

7 

20.9 

689.7 

5 

35 

21.24 

743.4 

5 

21.07714286 

737.7 

7 

20.98181818 

692.4  | 

5 

34 

21.46764706 

729.9 

5 

34 

21.75882353 

739.8 

2 

22.91707317 

939.6 

5 

21.60882353 

734.7 
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Simulated  Weekend  Data  for  Leak  Rate  of  0  in. /minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

10 

8 

20.0625 

160.5 

10 

7 

19.92857143 

139.5 

10 

7 

19.92857143 

139.5 

io 

8 

20.0625 

160.5 

10 

7 

19.92857143 

139.5 

10 

8 

20.0625 

160.5 

10 

7 

19.92857143 

139.5 

io 

8 

19.8375 

158.7 

7 

19.92857143 

139.5 

ii 

6 

20.1 

120.6 

li 

7 

19.92857143 

139.5 

ii 

7 

19.92857143 

139.5 

13 

6 

20.1 

120.6 

13 

6 

20.1 

120.6 

13 

5 

19.86 

99.3 

13 

5 

19.86 

99.3 

13 

6 

20.1 

120.6 

13 

6 

20.1 

120.6 

10 

8 

20.0625 

160.5 

10 

7 

19.92857143 

139.5 

10 

8 

20.0625 

160.5 

10 

7 

19.92857143 

139.5 

10 

6 

20.1 

120.6 

10 

7 

19.92857143 

139.5 
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Simulated  Weekend  Data  for  Leak  Rate  of  0.5  in. /minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

!0 

13 

18.27692308 

237.6 

13 

18.55384615 

241.2  ] 

14 

18.34285714 

256.8 

| 

13 

18.99230769 

246.9 

10 

13 

19.08461538 

248.1 

10 

14 

18.38571429 

257.4 

10 

18.75 

262.5 

10 

sa 

18.975 

227.7 

10 

14 

18.34285714 

256.8  j 

11 

13 

18.9 

245.7  ! 

11 

12 

19 

228 

11 

18.64615385 

242.4 

13 

12 

18.35 

220.2 

13 

12 

18.9 

226.8 

13 

12 

18.375 

220.5 

13 

13 

18.41538462 

239.4 

13 

12 

18.1 

217.2 

13 

12 

18.375 

220.5 

10 

13 

19.26923077 

250.5 

10 

13. 

18.39230769 

239.1 

10 

18.46153846 

240 

mmmm 

13 

18.73846154 

243.6 

10 

12 

18.575 

222.9 

lb 

14 

18.55714286 

259.8 
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Simulated  Weekend  Data  for  Leak  Rate  of  1.0  in./minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

19 

353.7 

10 

18 

18.65 

335.7 

10 

18 

18.93333333 

340.8 

10 

18 

19.65 

353.7 

10 

18 

18.68333333 

336.3 

10 

18 

19.31666667 

347.7 

10 

18 

18.9 

340.2 

10 

18 

18.93333333 

340.8 

18 

18.95 

341.1 

11 

17 

19.27058824 

327.6 

11 

17 

18.77647059 

319.2 

11 

18 

18.98333333 

341.7 

13 

18 

18.71666667 

336.9 

13 

17 

18.98823529 

322.8 

13 

16 

19.33125 

309.3 

13 

18 

18.46666667 

332.4 

13 

17 

18.88235294 

321 

13 

18 

19.25 

346.5 

mmmm 

18 

18.93333333 

340.8 

18 

18,85 

339.3 

10 

19 

18.41052632 

349.8 

10 

18 

19.U666667 

344.1 

10 

18 

18.9 

340.2 

10 

18 

18.76666667 

337.8 

Simulated  Weekend  Data  for  Leak  Rate  of  1.5  in. /minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

10 

23 

19.33043478 

444.6 

10 

23 

19.27826087 

443.4 

10 

22 

20.1 

442.2 

10 

23 

19.77391304 

454.8 

10 

23 

19.56521739 

450 

10 

23 

19.06956522 

438.6 

10 

23 

19.39565217 

446.1 

10 

23 

19.04347826 

438 

10 

22 

19.99090909 

439.8 

11 

23 

19.21304348 

441.9 

11 

23 

19.09565217 

439.2 

11 

23 

19.27826087 

443.4 

13 

22 

19.62272727 

431.7 

13 

21 

19.2 

403.2 

13 

21 

19.74285714 

414.6 

13 

21 

19.34285714 

406.2 

13 

22 

18.83181818 

414.3 

13 

22 

19.2 

422.4 

10 

22 

20.07272727 

441.6 

10 

23 

19.72173913 

453.6 

10 

23 

19.47391304 

447.9 

10 

22 

19.74545455 

434.4 

10 

24 

18.6625 

447.9 

10 

23 

19.56521739 

450 
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Simulated  Weekend  Data  for  Leak  Rate  of  2.0  in./minute 


Mean  Time 

Between 

Flushes 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

10 

mi 

20.15555556 

544.2 

10 

27 

19.8 

534.6 

27 

19.97777778 

539.4 

27 

20.02222222 

540.6 

10 

19.85357143 

555.9 

10 

19.71428571 

552 

10 

28 

19.97142857 

559.2 

10 

26 

20.97692308 

545.4 

26 

20.08846154 

522.3 

11 

26 

20.58461538 

535.2 

11 

27 

19.96666667 

539.1 

11 

20.11111111 

543 

13 

20.66538462 

537.3 

13 

26 

19.95 

518.7 

13 

26 

19.83461538 

515.7 

13 

27 

20 

540 

13 

26 

19.76538462 

513.9 

13 

27 

18.96666667 

512.1 

10 

27 

537 

10 

27 

20.54444444 

554.7 

10 

27 

20.14444444 

543.9 

10 

27 

20.1 

542.7 

10 

26 

20.22692308 

525.9 

10 

27 

20.4 

550.8 

Simulated  Weekend  Data  for  Leak  Rate  of  2.5  in./minute 


Mean  Time 

Between 

Flushes 

Runs 

Time  per  run  (seconds) 

Total  runtime  (seconds) 

10 

31 

20.66129032 

640.5 

10 

31 

21.09677419 

654 

10 

31 

20.83548387 

645.9 

10 

31 

21 

651 

10 

31 

21.07741935 

653.4  1 

10 

31 

20.81612903 

645.3 

10 

31 

21.0483871 

652.5 

10 

32 

20.71875 

663 

10 

Bl 

20.39032258 

632.1 

11 

KB 

20.50645161 

635.7  i 

11 

mm 

20.24516129 

627.6 

11 

mm 

20.49677419 

635.4 

13 

mm 

20.27419355 

628.5 

13 

K2H 

20.98 

629.4 

13 

30 

20.61 

618.3 

13 

mm 

20.63225806 

639.6 

13 

— 

20.81 

624.3 

13 

■PM 

20.52 

615.6 

10 

B 

20.79677419 

644.7 

EQHHHHI 

32 

20.04375 

641.4 

31 

20.81612903 

645.3 

10 

31 

20.87419355  . 

647.1 

10 

31 

20.7 

641.7 

10 

31 

20.92258065 

648.6 
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Appendix  E 

Fluid  System  MATLAB  Scripts 


This  appendix  contains  scripts  used  to  simulate  operation  of  a  centrifugal  pump.  The  first  three 
scripts  are  used  to  calculate  the  currents,  voltages  and  torque  used  by  the  motor.  The  next  three 
scripts  perform  a  MATLAB  version  of  the  “prep”  software  used  for  real  current  and  voltage  data 
to  develop  a  power  envelope  for  the  simulated  data.  The  final  script  puts  all  of  the  scripts  together 
to  run  the  simulation  and  plot  the  real  power  used  by  the  electric  motor  vs.  time.  The  scripts 
indparam.m,  ind.m,  and  convind.m  were  written  by  Steven  B.  Leeb  as  part  of  [15].  The  scripts 
prep_fft.m,  prep  .rotate  .m,  and  matlabprep.m  were  written  by  Robert  Cox,  a  fellow  student  at 
MIT’s  Laboratory  for  Electromagnetic  and  Electronic  Systems  (LEES).  The  script  indsim.m  was 
written  by  the  author. 


E.l  Indparam 


%  This  script  loads  the  machine  parameters  for  a  3hp,  180  V  (L-N,  peak)  AC 
%  induction  machine  into  the  global  environment . 

%  Run  this  code  to  load  sample  machine  parameters  before  simulating  with 
%  ind . 

7. 

7,  This  software  is  distributed  in  the  hope  that  it  will  be  useful,  but 
7.  WITHOUT  ANY  WARRANTY.  It  is  for  educational  use  only.  Please  do  not 
7o  distribute  or  sell  this  software,  or  remove  the  copyright  notice. 

1 

7.  Copyright,  1995,  1997,  1998,  2000  Steven  B.  Leeb 

global  P  rs  rr  Xm  Xls  Xlr  we  J  B1  vds  vqs  vqr  vdr  T1 

7.  These  are  the  machine  parameters  for  a  3  Hp,  180V  (L-N,  peak)  AC 

7.  induction  machine . 


P  =  4; 
rs  =  0.435; 
rr  =  0.816; 
Xm  =  26.13; 
Xls  =  0.754; 
Xlr  =  0.754; 
we  =  377 . ; 

J  =  0.089; 


7.  Number  of  poles  (*not*  pole  pairs) 

7.  Stator  resistance 
7.  Rotor  resistance 

7.  Magnetizing  Impedance,  in  Ohms  on  a  60  Hz  base 
7.  Stator  Side  Leakage 1  Impedance,  in  Ohms  on  a  60  Hz  base 
7.  Rotor  Side  Leakage  Impedance,  in  Ohms  on  a  60  Hz  base 
7.  Base  electrical  frequency,  rads  per  second  (60  Hz) 

7.  Rotor  Inertia 
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B1  =  0;  %  Load  Damping  Coefficient 

vds  =0*0;  %  D  axis  stator  voltage 

vqs  =  180*0;  %  Q  axis  stator  voltage 

vqr  =0*0;  %  D  axis  rotor  voltage 

vdr  =0*0;  %  0  axis  rotor  voltage 

T1  =  0.0;  %  Load  torque 


E.2  Ind 

function  [slopes]  =  ind(t,statev) 

%  [slopes]  =  ind(t,statev) 

%  .  1  .  .  '  '  ;  .  .  .  ; 

%  TMs  script  computes  the  state  variable  derivatives  for  a  fifth 

%  order  model  of  a  balanced,  three  phase  AC  induction  machine. 

%  The  state  variables  are  the  D  and  Q  stator  and  rotor  fluxes,  and 
%  the  rotor  speed  (wr) . 

%  Use  this  script  with  0DE45  to  simulate  the  performance  of  the 
%  induction  machine.  '■ 

% 

%  Run  indparam.m  to  load  sample  machine  parameters  before  simulating. 

% 

%  This  software  is  distributed  in  the  hope  that  it  will  be  useful,  but 
%  WITHOUT  ANY  WARRANTY.  It  is  for  educational  use  only.  Please  do  not 
%  distribute  or  sell  this  software,  or  remove  the  copyright  notice. 

% 

%  Copyright,  1995,  1998,  2002  Steven  B.  Leeb 

%  The  variable  w  determines  the  reference  frame  in  which  the  simulation 
%  will  be  conducted.  With  w  =  377,  the  simulation  variables  will  be  in 
%  a  synchronously  rotating  reference  frame, 
global  P  rs  rr  Xm  Xls  Xlr  we  J  R1  vds  vqs  vqr  vdr  T1 
w  =  377.; 

Lm  =  Xm/we;  %  This  is  the  magnetizing  inductance 
Lis  =  Xls/we;  %  Stator  leakage  * 

Llr  =  Xir/we;  %  Rotor  leakage 
Las  =  Lis  +  Lm; 

Lax  =  Llr  +  Lm; 

lamqs  =  statev(l) ; 
lamds  =  statev(2); 
lamqr  =  statev(3) ; 
lamdr  =  statev<4); 
wr  =  statev(5) ; 

%th  =  statev(6) ; 

D  =  Lm*Lm  -  Las*Lar ; 

idr  =  (Lm*lamds  -  Las*lamdr)/D; 

iqr  =  (Lm*lamqs  -  Las*lamqr)/D; 
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iqs  =  (Lm*lamqr  -  Lar*lamqs)/D; 
ids  =  (Lm*lamdr  -  Lar*lamds)/D; 

si  =  (vqs  -  w*lamds  -  rs*iqs) ; 
s2  =  (vds  +  w*lamqs  -  rs*ids) ;  - 

s3  =  (vqr  -  (w  -  wr)*lamdr  -  rr*iqr) ; 
s4  =  (vdr  +  (w  -  wr)*lamqr  -  rr*idr) ; 

%  P  is  the  number  of  poles,  *not*  pole  pairs 
T  =  (3/2)*(P/2)*(lamqr*idr  -  lamdr*iqr) ; 
s5  =  (P/2) * (T-0 . 0002614*wr~2)/ ( J  +  0.00401); 

%s6  =  377; 

slopes  =  [si  s2  s3  s4  s5]’;%  s6] ’ ; 

E.3  Convind 

function  [m,m2]  =  convind (t  ,y)  ; 

1  [m,m2]  =  convind (t ,statev) 

°/o 

7o  This  script  transforms  the  DQ  stator  and  rotor  fluxes  computed  using 
%  ind.m  and  ode45  back  into  laboratory  frame  stator  currents  and  voltages, 
%  e.g.,  ias  and  vas  for  phase  a. 

% 

7o  The  variable  w  determines  the  reference  frame  in  which  the  simulation 
7a  will  be  conducted.  With  w  =  377,  the  simulation  variables  will  be  in 
7o  a  synchronously  rotating  reference  frame. 

% 

7o  On  return,  the  output  matrices  m  and  m2  contain: 

7o  m  =  [ids  iqs  idr  iqr]  ; 

7®  m2  =  [T  ias  ibs  ics  vas  vbs  vcs]  ; 

7o  This  script  also  plots  the  simulated  rotor  torque  versus  speed  on  return 

7o 

7a  This  software  is  distributed  in  the  hope  that  it  will  be  useful,  but 
7a  WITHOUT  ANY  WARRANTY.  It  is  for  educational  use  only.  Please  do  not 
7a  distribute  or  sell  this  software,  or  remove  the  copyright  notice. 

7a 

7o  Copyright,  1995,  1998,  2002  Steven  B.  Leeb 
global  P  rs  rr  Xm  Xls  Xlr  we  J  B1  vds  vqs  vqr  vdr  T1 
7o  select  reference  frame  (see  ind.m) 
w  =  377. ; 

Lm  =  Xm/we;  7o  This  is  the  magnetizing  inductance 
Lis  =  Xls/we;  7o  Stator  leakage 

Llr  =  Xlr/we;  7o  Rotor  leakage 

Las  =  Lis  +  Lm; 

Lar  =  Llr  +  Lm; 

.lamqs  =  y(:  ,1)  ; 
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lamds  =  y( : ,2) ; 
lamqr  =  y(: ,3) ; 
lamdr  =  y(: ,4) ; 
wr  =  y(:,5)j 
th  =  w*t;  %y  C : , 6) ; 

D  =  Lm*Lm  -  Las*Lar ; 
idr  =  (Lm*lamds  -  Las*lamdr) /D ; 

iqr  =  (Lm*lamqs  -  Las*lamqr)/D; 

iqs  =  (Lm*lamqr  -  Lar*lamqs)/D; 

ids  =  (Lm*lamdr  -  Lar*lamds)/D; 

T  =  (3/2)*(P/2)*(lamqr.*idr  -  lamdr. *iqr) ; 

ias  =  cos(th) .*ids  -  sin(th) .*iqs; 

ibs  =  cos(th  -  2.*pi/3.) .*ids  -  sinCth  -  2.*pi/3.) .*iqs; 

ics  =  cos(th  +  2.*pi/3.) ,*ids  -  sinCth  +  2.*pi/3.) .*iqs; 

vas  =  cos(th).*vds  -  sin(th) . *vqs; 

vbs  =  cosCth  -  2 . *pi/3 . ) . *vds  -  sinCth  -  2 . *pi/3 . ) . *vqs ; 

vcs  =  cosCth  +  2 . *pi/3 . ) . *vds  -  sinCth  +  2.*pi/3.) .*vqs; 

m  =  [ids  iqs  idr  iqr] ; 

m2  =  [T  ias  ibs  ics  vas  vbs  vcs] ; 

E.4  PrepJFFT 

%  This  file  will  run  an  offline  version  of  prep  on  data  that  has  already  been  sampled 
%  to  have  128  pts/period 
%  A  is  current  data 
%  B  is  voltage  data 

function  [vHarm,  iHarm]  =  prep_fft(A,  B,  rows) 

i  -  i;  ; 

j  =  1; 

nBasis  *  128; 
shift  =  64; 

harm  =  [1357];  %  The  harmonics  we  want  to  look  at 
nHarm  =  4; 

in  =  zeros (l^nBasis) ; 
out  =  zeros (1, nBasis) ; 
iHarm  -  zeros (1? rows *nHarm*2) ; 
vHarm  =  zeros (1, rows*2) ; 

for  i=l:l:rows 

for  j=l:l:nBasis 
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in(j)  =  B(j  +  (i  -  1)  *  shift);  %  load  up  the  vector  to  be  FFTd 

end 

out  =  fft(in,  nBasis) ; 
vHarm(i*2-l)  =  real (out (2) ) ; 
vHarm(i*2)  =  imag(out(2) ) ; 
for  j=l:l:nBasis 

in(j)  =  A(j+(i-l)*shift) ; 

end 

out  =  fft(in, nBasis) ; 
for  j=l:l:nHarm 

iHarm((i-l)*2*nHarm+((j*2)-l))  =  real (out (harm(j )+l) ) ; 

-  iHarm((i-l)*2*nHarm  +  2*j)  =  imag(out(harm(j)+l)) ; 
end 

end 


E.5  Prep  _Rotate 

function  E  =prep_rotate (vHarm, iHarm,rows) 

pqbax  =  zeros (1,2); 
pqtmp  =  zeros (1,2); 

i=l; 

j=i; 

nrot=l ; 

nHarm  =4; 

harm  =  [1  3  5  7]; 

E  =  zeros(rows,2*nHarm) ; 

for  i=l:l:rows; 

pqbar(l)  =  vHarm(i*2  -  1); 
pqbar(2)  =  vHarm (i*2)  ; 

mag  =  sqrt (pqbar (1)  *  pqbar(l)  +  pqbar(2)  *  pqbar(2)); 

if  (mag  ~=  0) 

pqbar (1)  =  pqbar (l)/mag; 
pqbar (2)  =  pqbar (2) /mag; 

end 

nrot  =1; 

pqtmp  ( 1 )  =  pqbar  (1)  ; 
pqtmp  (2)  =  pqbar  (2); 


for  j=l:l:nHarm 
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a  =  iHarm((i-l)*2*nHarm+((j*2)-l)) ; 

b  =  iHarm((i-l)*2*nHarm  +  2*j);  • 

pqtmp  =  prep_multiply(pqtmp,  pqbar,  nrot,  harm(j)) ; 

E(i,(2*j)-1)  =  a*pqtmp(l)  +  b*pqtmp(2); 

E(i,2*j)  =  -(b*pqtmp(l)  -  a*pqtmp(2)) ; 

end 

end 

E.6  Matlabprep  / 

%  E  =  matlabprep(A,B,rows) 

% 

%  Given  current  and  voltage  sampled  with  128  points  per  period,  matlabprep  will 
%  output  arbitrarily  scaled  spectral  envelopes. 

% 

%  A  is  the  vector  of  current  values 
%  B  is  the  vector  of  voltage  values 
%  rows  is  the  number  of  cycles  of  the  60Hz  waveform 

%  '  1  i  ^  . 

%  This  function,  and  all  of  the  functions  which  it  invokes,  are 
%  based  on  routines  found  in  softprep.c  by  Steve  Shaw. 

% 

%  Robert  Cox,  2003 

function  E=matlabprep(A,B,rows) 

[vHarm , iHarm] =pr ep_f f t (A , B , rows) ; 

E=pr ep^rotat e (vHarm , iHarm , rows ) ; 

E.7  Indsim 

indparam; 

t=0: 1/(128*60): 8; 

[t ,y]=ode45(@ind,t, [0  000  0]); 

wr=y(: ,5) ;  ^ 

[m,m2]=convind(t,y) ;  ; 

Ias=m2(:,2); 

%Vas=m2(: ,5); 

Vcb=m2( : ,7)-m2( : ,6) ; 

rows=480 ;  \ 

time=(l : 481) /60 ; 

E=ones(481,8) ; 

E(l, :)=0; 

E(2:481, :)=matlabprep(Ias,Vcb,rows) ; 
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P=E( : ,2) ; 
plot (time, P) 


